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Abstract 

Limit Theorems for Internal Aggregation Models 

by 

Lionel Timothy Levine 
Doctor of Philosophy in Mathematics 

University of California, Berkeley 

Professor Yuval Peres, Chair 

We study the scaling limits of three different aggregation models on Z*^: internal DLA, in 
which particles perform random walks until reaching an unoccupied site; the rotor-router 
model, in which particles perform deterministic analogues of random walks; and the divisible 
sandpile, in which each site distributes its excess mass equally among its neighbors. As the 
lattice spacing tends to zero, all three models are found to have the same scaling limit, 
which we describe as the solution to a certain PDE free boundary problem in R'^. In 
particular, internal DLA has a deterministic scaling limit. We find that the scaling limits 
are quadrature domains, which have arisen independently in many fields such as potential 
theory and fiuid dynamics. Our results apply both to the case of multiple point sources and 
to the Diaconis- Fulton smash sum of domains. 

In the special case when all particles start at a single site, we show that the scaling 
limit is a Euclidean ball in M.'^ and give quantitative bounds on the rate of convergence to a 
ball. For the divisible sandpile, the error in the radius is bounded by a constant independent 
of the total starting mass. For the rotor-router model in Z*^, the inner error grows at most 
logarithmically in the radius r, while the outer error is at most order r^~^/'^logr. We also 
improve on the previously best known bounds of Le Borgne and Rossin in and Fey and 
Redig in higher dimensions for the shape of the classical abelian sandpile model. 

Lastly, we study the sandpile group of a regular tree whose leaves are collapsed 
to a single sink vertex, and determine the decomposition of the full sandpile group as a 
product of cyclic groups. For the regular ternary tree of height n, for example, the sandpile 
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group is isomorphic to (^3)^" (Zy)^" © ... Z2n-i_i © Z2n-i. We use this result to 
prove that rotor-router aggregation on the regular tree yields a perfect ball. 
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Chapter 1 



Introduction 



1.1 Three Models with the Same ScaUng Limit 

Given finite sets A,Bc Z'^, Diaconis and Fulton [15] defined the smash sum A(BB 
as a certain random set whose cardinality is the sum of the cardinalities of A and B. Write 
An B = {xi, . . . , Xk}- To construct the smash sum, begin with the union Co = AU B and 
for each j = 1 , . . . , A: let 

Cj = Cj-iU{yj} 

where yj is the endpoint of a simple random walk started at xj and stopped on exiting Cj-i. 
Then define A(B B = Ck- The key observation of [15] is that the law A® B does not 
depend on the ordering of the points Xj. The sum of two squares in J? overlapping in a 
smaller square is pictured in Figure p3| 



In Theorem 1.1.3 below, we prove that as the lattice spacing goes to zero, the 
smash sum A® B has a deterministic scaling limit in W^. Before stating our main results, 
we describe some related models and describe our technique for identifying their common 
scaling limit, which comes from the theory of free boundary problems in PDE. 

The Diaconis-Fulton smash sum generalizes the model of internal diffusion-limited 
aggregation ( "internal DL A" ) studied in |29j , and in fact was part of the original motivation 
for that paper. In classical internal DLA, we start with n particles at the origin o G and 
let each perform a simple random walk until it reaches an unoccupied site. The resulting 
random set of n occupied sites in Z*^ can be described as the n-fold smash sum of {o} 
with itself. We will use the term internal DLA to refer to particles which perform simple 



Figure 1.1: Smash sum of two squares in 1? overlapping in a smaller square, for internal 
DLA (top left), the rotor-router model (top right), and the divisible sandpile. 

random walks in until reaching an unoccupied site, starting from an arbitrary initial 
configuration. In this broader sense of the term, both the Diaconis-Fulton sum and the 
model studied in [29J are particular cases of internal DLA. 

In defining the smash sum A©S, various alternatives to random walk are possible. 
Rotor-router walk is a deterministic analogue of random walk, first studied by Priezzhev et 
al. [38] under the name "Eulerian walkers." At each site in I? is a rotor pointing north, 
south, east or west. A particle performs a nearest-neighbor walk on the lattice according 
to the following rule: during each time step, the rotor at the particle's current location is 
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rotated clockwise by 90 degrees, and the particle takes a step in the direction of the newly 
rotated rotor. In higher dimensions, the model can be defined analogously by repeatedly 
cycling the rotors through an ordering of the 2d cardinal directions in Z"^. The sum of two 
squares in I? using rotor-router walk is pictured in Figure all rotors began pointing 
west. The shading in the figure indicates the final rotor directions, with four different shades 
corresponding to the four possible directions. 

The divisible sandpile model uses continuous amounts of mass in place of discrete 
particles. A lattice site is full if it has mass at least 1. Any full site can topple by keeping 
mass 1 for itself and distributing the excess mass equally among its neighbors. At each 
time step, we choose a full site and topple it. As time goes to infinity, provided each full 
site is eventually toppled, the mass approaches a limiting distribution in which each site 
has mass < 1. Note that individual topplings do not commute. However, the divisible 
sandpile is "abelian" in the sense that any sequence of topplings produces the same limiting 
mass distribution; this is proved in Lemma 2.2.1, Figure [LT] shows the limiting domain of 
occupied sites resulting from starting mass 1 on each of two squares in Z^, and mass 2 on 
the smaller square where they intersect. 

Figure [LT] raises a few natural questions: as the underlying lattice spacing becomes 
finer and finer, will the smash sum A (B B tend to some limiting shape in M'^, and if so, 
what is this shape? Will it be the same limiting shape for all three models? To see how we 
might identify the limiting shape, consider the divisible sandpile odometer function 

u(x) = total mass emitted from x. (1-1) 

Since each neighbor y ~ x emits an equal amount of mass to each of its 2d neighbors, the 
total mass received by x from its neighbors is Yly^x'^iv)^ hence 

Au{x) = iy{x) - a{x) (1.2) 

where (t{x) and i^^x) are the initial and final amounts of mass at x, respectively. Here A is 
the discrete Laplacian in Z'^, defined by 

^u{x) = ^^u{y) -u{x). (1.3) 



Equation (1.2) suggests the following approach to finding the limiting shape. We 
first construct a function on Z'^ whose Laplacian is o" — 1; an example is the function 

7(x) = -\x\^ - ^ gi{x,y)a{y) (1.4) 
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where in dimension d > 3 the Green's function gi{x,y) is the expected number of times 
a simple random walk started at x visits y (in dimension d = 2 we use the recurrent 
potential kernel in place of the Green's function). Here \x\ denotes the Euclidean norm 



(xf + ... + x2)i/2. By (1.2), since u < 1 the sum u + 7 is a superharmonic function on Z ; 
that is, A(u + 7) < 0. Moreover if / > 7 is any superharmonic function lying above 7, 
then / — 7 — li is superharmonic on the domain D = {x ^ U^\v{x) = 1} of fully occupied 
sites, and nonnegative outside D, hence nonnegative everywhere. Thus we have proved the 
following lemma. 

Lemma 1.1.1. Let a he a nonnegative function on Z'^ with finite support. Then the odome- 



ter function (1.1) for the divisible sandpile started with mass (y{x) at each site x is given 
by 

u = s — 7 



where 7 is given by (I.4), o,nd 

s{x) = inf{/(x)|/ is superharmonic on Wf" and / > 7} 
is the least superharmonic majorant ofj. 



Lemma 1.1.1 allows us to formulate the problem in a way which translates naturally 
to the continuum. Given a function a on M"^ representing the initial mass density, by analogy 



with (1.4) we define the obstacle 

7(x)=-|x|^-/ g{x,y)a{y)dy 

where g{x,y) is the Green's function on M'^ proportional to \x — y\'^~'^ in dimensions d > 3 
and to — log |x — y| in dimension two. We then let 

s(x) = inf{/(x)|/ is continuous, superharmonic and / > 7}. 

The odometer function for a is then given by u = s — 7, and the final domain of occupied 
sites is given by 

D = {x G M'^|s(x) > 7(x)}. (1.5) 

This domain D is called the noncoincidence set for the obstacle problem with obstacle 7; 
for an in-depth discussion of the obstacle problem, see [T9j. 



Figure 1.2: The obstacles 7 corresponding to starting mass 1 on each of two overlapping 
disks (top) and mass 100 on each of two nonoverlapping disks. 



If A, B are bounded open sets in M*^, we define the smash sum of A and B as 

A®B = A\JB[JD (1.6) 



where D is given by (1.5) with o" = 1^ + 1^. In the two-dimensional setting, an alternative 
definition of the smash sum in terms of quadrature identities is mentioned in |23] . 

In this thesis we prove, among other things, that if any of our three aggregation 
models is run on finer and finer lattices with initial mass densities converging in an appro- 
priate sense to a, the resulting domains of occupied sites will converge in an appropriate 



sense to the domain D given by (1.5). We will always work in dimension d > 2; for a 
discussion of the rotor-router model in one dimension, see [32] . 

Let us define the appropriate notion of convergence of domains, which amounts 
essentially to convergence in the Hausdorff metric. Fix a sequence 5n i representing the 
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lattice spacing. Given domains C Sn'^'^ and D C M*^, write An D if for any e > 



(1.7) 



for all sufficiently large n. Here 



A = {x £ D\B{x,€) C L>} 



and 

= {xe R'^\B{x,e)nD ^ 0} 

are the inner and outer e- neighborhoods of D. For x G Sn'^'^ we write x'-' = [x + ^ , a; — ^] 
For t gM write [t] for the closest integer to t. 

Throughout this thesis, to avoid trivialities we work in dimension d >2. Our first 
main result is the following. 

Theorem 1.1.2. Let 17 C M'^ 6e a bounded open set, and let a : M*^ Z>o be a bounded 
function which is continuous almost everywhere, satisfying {o" > 1} = 0. Let Dn,Rn,Ln 
be the domains of occupied sites formed from the divisible sandpile, rotor-router model, and 
internal DLA, respectively, in the lattice 6n'^'^ started from source density 



an{x) 



a{y)dy 



Then as n ^ oo 



Dn,Rn -^Dun; 
and if 5n < 1/ log n, then with probability one 



In ^ Dun 



where D is given by (1.5), and the convergence is in the sense of (1.1). 



Remark. When forming the rotor-router domains Rn, the initial rotors in each lattice 5n'^'^ 
may be chosen arbitrarily. 



We prove a somewhat more general form of Theorem 1.1.2 which allows for some 
fiexibility in how the discrete density an is constructed from a. In particular, taking a = 
1^ + 1_B obtain the following theorem, which explains the similarity of the three smash 
sums pictured in Figure |1.1[ 
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Theorem 1.1.3. Let A,B dW^ he hounded open sets whose houndaries have measure zero. 
Let Dn,Rn,Ln he the smash sum of Andnl^'^ and i3n5„Z'^, formed using divisihle sandpile, 
rotor-router and internal DLA dynamics, respectively. Then as n ^ oo 

Dn,Rn ^ A®B- 

and if 5n < 1/ log n, then with prohahility one 

Ir,.^ A®B 



where A(B B is given hy (1.6), and the convergence is in the sense of (1.1). 



For the divisible sandpile, Theorem 1.1.2 can be generalized by dropping the re- 



quirement that a be integer valued; see Theorem 3.2.7 for the precise statement. Taking 
a real-valued is more problematic in the case of the rotor-router model and internal DLA, 
since these models work with discrete particles. Still, one might wonder if, for example, 
given a domain A C W^, starting each even site in An 5„,Z'^ with one particle and each odd 
site with two particles, the resulting domains Rn,In would converge to the noncoincidence 
set D for density a = |1a- This is in fact the case: if a-a is a density on Sn'^'^, as long 
as a certain "smoothing" of o"„ converges to a, the rotor-router and internal DLA domains 



started from source density an will converge to D. See Theorems 3.3.7 and 3.4.1 for the 
precise statements. 



1.2 Single Point Sources 



One interesting case not covered by Theorems 1.1.2 and 1.1.3 is the case of point 
sources. Lawler, Bramson and Griffeath [29] showed that the scaling limit of internal DLA in 
Z'^ with a single point source of particles is a Euclidean ball. In chapter [2] we prove analogous 
results for rotor-router aggregation and the divisible sandpile, and give quantitative bounds 
on the rate of convergence to a ball. Let An be the domain of n sites in Z'^ formed from 
rotor-router aggregation starting from a point source of n particles at the origin. Thus An 
is defined inductively by the rule 



Ar. 



An^l U {Xn} 



where Xn is the endpoint of a rotor-router walk started at the origin in Z'^ and stopped on 
first exiting An-i- For example, in Z^, if all rotors initially point north, the sequence will 



Figure 1.3: Rotor-router aggregate of one million particles in 1? . Each site is colored 
according to the direction of its rotor. 

begin = {(0,0)}, A2 = {(0, 0), (1, 0)}, A3 = {(0, 0), (1, 0), (0, -1)}. The region ^1,000,000 
is pictured in Figure |1.3[ 

Jim Propp observed from simulations in two dimensions that the regions An are 
extraordinarily close to circular, and asked why this was so |26| I39j. Despite the impressive 
empirical evidence for circularity, the best result known until now [S^ says only that if 
An is rescaled to have unit volume, the volume of the symmetric difference of An with a 
ball of unit volume tends to zero as a power of n, as n | 00. The main outline of the 
argument is summarized in [3l] . Fey and Redig ^18j also show that An contains a diamond. 
In particular, these results do not rule out the possibility of "holes" in An far from the 
boundary or of long tendrils extending far beyond the boundary of the ball, provided the 
volume of these features is negligible compared to n. 

Our main result on the shape of rotor-router aggregation with a single point source 
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is the following, which rules out the possibility of holes far from the boundary or of long 
tendrils. For r > let 

Br = {xeZ'^ : |x| < r}. 



Theorem 1.2.1. Let An be the region formed by rotor-router aggregation in Z'^ starting 
from n particles at the origin and any initial rotor state. There exist constants c, c' depending 
only on d, such that 

Br—clogr C An C B^^^^^,^ — i/dlogr) 

where r = {n/uJdY^'^ , and uJd is the volume of the unit ball in W^. 

We remark that the same result holds when the rotors evolve according to stacks 



of bounded discrepancy; see the remark following Lemma 2.4.1 



By way of comparison with Theorem 1.2.1 , if /„ is the internal DLA region formed 



from n particles started at the origin, the best known bounds |30] are (up to logarithmic 
factors) 

Bj._r^/3 C In C i?^_|_^l/3 

for all sufficiently large n, with probability one. 

Our next result treats the divisible sandpile with all mass initially concentrated at 
a point source. The resulting domain of fully occupied sites is extremely close to a ball; in 
fact, the error in the radius is bounded independent of the total mass. 

Theorem 1.2.2. For m > let C Ij'^ be the domain of fully occupied sites for the 
divisible sandpile formed from a pile of mass m at the origin. There exist constants c, c' 
depending only on d, such that 

— c C Dm C Bj.-\-c' 1 

where r = (m/ijJdY^'^ and uJd is the volume of the unit ball in M'^. 

The divisible sandpile is similar to the "oil game" studied by Van den Heuvel |39] . 
In the terminology of [18] , it also corresponds to the /i — > — oo limit of the classical abelian 
sandpile (defined below), that is, the abelian sandpile started from the initial condition in 
which every site has a very deep "hole." 
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Figure 1.4: Classical abelian sandpile aggregate of one million particles in 1? . Colors 
represent the number of grains at each site. 

In the classical abelian sandpile model , each site in Z"^ has an integer number of 
grains of sand; if a site has at least 2d grains, it topples, sending one grain to each neighbor. 
If n grains of sand are started at the origin in Z"^, write 5„ for the set of sites that are 
visited during the toppling process; in particular, although a site may be empty in the final 
state, we include it in Sn if it was occupied at any time during the evolution to the final 
state. 

Until now the best known constraints on the shape of Sn in two dimensions were 
due to Le Borgne and Rossin [31 , who proved that 

{x G |xi +X2 < - 1} C 5„ C {x € | xi,X2 < \fii/2}. 

Fey and Redig [TH] proved analogous bounds in higher dimensions, and extended these 
bounds to arbitrary values of the height parameter h. This parameter is discussed in 
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Figure 1.5: Known bounds on the shape of the classical abelian sandpile in 1?. The inner 
diamond and outer square are due to Le Borgne and Rossin [31j ; the inner and outer circles 
are those in Theorem 



1.2.3 



section 12.31 

The methods used to prove the near-perfect circularity of the divisible sandpile 



shape in Theorem 1.2.2 can be used to give constraints on the shape of the classical abelian 
sandpile, improving on the bounds of [18] and j3T] . 

Theorem 1.2.3. Let Sn be the set of sites that are visited by the classical abelian sandpile 
model in Z'^, starting from n particles at the origin. Write n = w^r"^. Then for any e > 
we have 



cir-C2 



where 



ci = {2d-l) 



-i/d 



id 



-i/d 



The constant C2 depends only on d, while depends only on d and e. 
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Note that Theorem 1.2.3 does not settle the question of the asymptotic shape of 



Sn, and Figure 1.4 indicates that the limiting shape in two dimensions may be a polygon 
rather than a disk. Even the existence of an asymptotic shape is not known, however. 



1.3 Multiple Point Sources 



Using our results for single point sources together with the construction of the 



smash sum (1.6), we can understand the limiting shape of our aggregation models started 



from multiple point sources. The answer turns out to be a smash sum of balls centered at 
the sources. For a; G M'* write x'' for the closest lattice point in Sn'^'^, breaking ties to the 



right. Our shape theorem for multiple point sources, which is deduced from Theorems 1.1.3 



1.2.1 and 1.2.2 along with the main result of [2.9j, is the following. 

Theorem 1.3.1. Fix xi, . . . ,Xk G M*^ and Ai, . . . , Afc > 0. Let Bi he the hall of volume Aj 
centered at Xj. Fix a sequence 6n I 0, and for x £ 6n1j'^ let 

k 
i=l 



an{x) 
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Let Dn,Rn,In be the domains of occupied sites in 6n'^'^ formed from the divisible sandpile, 
rotor-router model, and internal DLA, respectively, started from source density an- Then 
as n ^ CO 

Dn,Rn^ Bi® ...®Bk; (1.9) 
and if 5n < l/ra, then with probability one 



where © denotes the smash sum (1.6), and the convergence is in the sense of (1.1). 



Implicit in equation (1.9) is the associativity of the smash sum operation, which is 



not readily apparent from the definition ( 1.6 ). For a proof of associativity, see Lemma 3.5.1 
For related results in dimension two, see |4H Prop. 3.10] and [ 501 section 2.4]. 

We remark that a similar model of internal DLA with multiple point sources was 
studied by Gravner and Quastel |21 , who also obtained a variational solution. In their 
model, instead of starting with a fixed number of particles, each source Xj emits particles 
according to a Poisson process. The shape theorems of [21] concern convergence in the sense 



of volume, which is a weaker form of convergence than (1.7). 



1.4 Quadrature Domains 

By analogy with the discrete case, we would expect that volumes add under the 
smash sum operation; that is, 

C{A®B) = C{A) + C{B) 
where C denotes Lebesgue measure in W^. Although this additivity is not immediately 



apparent from the definition (1.6), it holds for all bounded open A,B<zW^ provided their 



boundaries have measure zero; see Corollary [3.1.14 



We can derive a more general class of identities known as quadrature identities 
involving integrals of harmonic functions over A ® B. Let us first consider the discrete 
case. If /i is a superharmonic function on Z'^, and o" is a mass configuration for the divisible 
sandpile (so each site x G Z'^ has mass o"(x)), the sum X^^g^d h{x)a{x) can only decrease 
when we perform a toppling. Thus 

^ h{x)u{x) < ^ h{x)(T{x), (1.10) 



14 



where v is the final mass configuration. We therefore expect the domain D given by (1.5) 
to satisfy the quadrature inequality 



D 



h{x)dx < I h{x)a{x)dx 



(1.11) 



D 



for all integrable superharmonic functions h on D. For a proof under suitable smoothness 



assumptions on a and /i, see Proposition 3.1.11 see also 

A domain D C M'^ satisfying an inequality of the form ( 1.11 ) is called a quadrature 
domain for a. Such domains are widely studied in potential theory and have a variety of 
appUcations in fluid dynamics HO] . For more on quadrature domains and their connection 



with the obstacle problem, see [H [HI [Ml [23 HH [12]- Equation (1.10) can be regarded as 
a discrete analogue of a quadrature inequality; in this sense our aggregation models, in 
particular the divisible sandpile, produce discrete analogues of quadrature domains. 

In Proposition |3.5.5[ we show that the smash sum of balls Bi® . . .® Bk arising in 



Theorem 1 1 . 3 . 1 1 obeys the classical quadrature identity 

k 

h{x)dx = Aj/i(a 



(1.12) 



i=l 



'Bi®...®Bk 

for all harmonic functions h on Bi (B ■ ■ ■ (B B^. This can be regarded as a generalization 
of the classical mean value property of harmonic functions, which corresponds to the case 
k = 1. Using results of Gustafsson and Sakai [42 on quadrature domains in the plane, 



we can deduce the following theorem, which is proved in section 3.5 

Theorem 1.4.1. Let Bi, . . . , B^ be disks in with distinct centers. The boundary of the 
smash sum Bi (B ■ ■ ■ (B Bk lies on an algebraic curve of degree 2k. More precisely, there is 
a polynomial P £ M[xi, X2] of the form 

P{xi,X2) = [xf + + lower order terms 

and there is a finite set of points E C M^, possibly empty, such that 

d{Bi ®...(BBk) = {{xi,X2) G R^\P{xi,X2) =0}-E. 

For example, if Bi and B2 are disks of equal radius r > 1 centered at (1,0) and 
(—1,0), then d{Bi © B2) is given by the quartic curve HI] 

(1.13) 



(xl + xl)^ - 2r^ {xl + xl) - 2{xl - xj) = 0. 



This curve describes the shape of the rotor-router model with two point sources pictured in 
Figure [1.6 
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1.5 Aggregation on Trees 

Let T be the infinite d-regular tree. To define rotor-router walk on a tree, for each 
vertex of T choose a cyclic ordering of its d neighbors. Each vertex is assigned a rotor 
which points to one of the neighboring vertices, and a particle walks by first rotating the 
rotor at each site it comes to, then stepping in the direction of the newly rotated rotor. 
Fix a vertex o £ T called the origin. Beginning with Ai = {o}, define the rotor-router 
aggregation cluster An inductively by 

An = An^l U {Xn}, n> 1 

where G T is the endpoint of a rotor-router walk started at o and stopped on first exiting 
An-i- We do not change the positions of the rotors when adding a new chip. Thus the 
sequence {An)n>i depends only on the choice of the initial rotor configuration. 

Our next result is the analogue of Theorem 1 1.2.1 on regular trees. Call a configu- 



ration of rotors acyclic if there are no directed cycles of rotors. On a tree, this is equivalent 
to forbidding directed cycles of length 2: for any pair of neighboring vertices t; ~ w, if the 
rotor at v points to w, then the rotor at w does not point to v. As the following result 
shows, provided we start with an acyclic configuration of rotors, the occupied cluster An is 
a perfect ball for suitable values of n. 

Theorem 1.5.1. Let T he the infinite d-regular tree, and let 

Br = {x £ T : |a;| < r} 

be the ball of radius r centered at the origin a G T, where \x\ is the length of the shortest 
path from a to x. Write 

br = #Br = l + d ^'^~^^^^~\ 

Let An be the region formed by rotor-router aggregation on the infinite d-regular tree, starting 
from n chips at a. If the initial rotor configuration is acyclic, then 

Ab^ = Br- 



The proof of Theorem 1.5.1 uses the sandpile group of a finite regular tree with 



the leaves collapsed to a single vertex. This is an abelian group defined for any graph G 



whose order is the number of spanning trees of G. In section 4.1 we recall the definition of 
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the sandpile group and prove the following decomposition theorem expressing the sandpile 
group of a finite regular tree as a product of cyclic groups. 

Theorem 1.5.2. Let Tn he the regular tree of degree d = a + 1 and height n, with leaves 
collapsed to a single sink vertex and an edge joining the root to the sink. Then writing Zp 
for the group (Z/pZ) © ... © {Ij/pli) with q summands, the sandpile group ofTn is given by 



n-l . 



For example, taking d = 3 we obtain that the sandpile group of the regular ternary 
tree of height n has the decomposition 



n — 4 



5P(r„)~(Z3)^ ©(Zr)^ ©...©Z2n-l_l©Z2n_l. 

Toumpakari [U] studied the sandpile group of the ball inside the infinite d- 
regular tree. Her setup differs slightly from ours in that there is no edge connecting the 
root to the sink. She found the rank, exponent, and order of SP{Bn) and conjectured a 
formula for the ranks of its Sylow p-subgroups. We use Theorem |1.5.2| to give a proof of 



her conjecture in section 4.1.3 



Chen and Schedler [TO] study the sandpile group of thick trees (i.e. trees with 
multiple edges) without collapsing the leaves to the sink. They obtain quite a different 
product formula in this setting. 



In section 4.2 we define the rotor-router group of a graph and show that it is 



isomorphic to the sandpile group. We then use this isomorphism to prove Theorem 1.5.1 

Much previous work on the rotor-router model has taken the form of comparing 
the behavior of rotor-router walk with the expected behavior of random walk. For example. 
Cooper and Spencer [12] show that for any configuration of chips on even lattice sites in Z'^, 
letting each chip perform rotor-router walk for n steps results in a configuration that differs 
by only constant error from the expected configuration had the chips performed independent 
random walks. We continue in this vein by investigating the recurrence and transience of 
rotor-router walk on trees. A walk which never returns to the origin visits each vertex 
only finitely many times, so the positions of the rotors after a walk has escaped to infinity 
are well-defined. We construct two "extremal" rotor configurations on the infinite ternary 
tree, one for which walks exactly alternate returning to the origin with escaping to infinity, 
and one for which every walk returns to the origin. The latter behavior is something of 
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a surprise: to our knowledge it represents the first example of rotor-router walk behaving 
fundamentally differently from the expected behavior of random walk. 

In between these two extreme variety of intermediate behaviors are possi- 

ble. We say that a binary word oi . . . a„ is an escape sequence for the infinite ternary tree if 
there exists an initial rotor configuration on the tree so that the A;-th chip escapes to infinity 
if and only if = 1. The following result characterizes all possible escape sequences on the 
ternary tree. 

Theorem 1.5.3. Let a = ai...an be a binary word. For j £ {1,2,3} write a^^^ = 
ajUj-^-^aj^Q . . .. Then a is an escape sequence for some rotor configuration on the infi- 
nite ternary tree if and only if for each j and k >2, every subword of a^^^ of length 2^ — 1 
contains at most 2^^^ ones. 



Theorem 1.5.3 is proved in section 4.2.3 by expressing the escape sequence corre- 
sponding to a rotor configuration on the full tree in terms of the escape sequences of the 
configurations on each of the principal subtrees. 
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Chapter 2 



Spherical Asymptotics for Point 
Sources 



This chapter is devoted to the proofs of Theorems 1.2.1 1.2.2 and 1.2.3 In sec 



tion 2.1, we derive the basic Green's function estimates that are used in the proofs. In 



section 2.2 we prove the abehan property and Theorem |1.2.2 for the divisible sandpile. In 



section 2.3 we adapt the methods used for the divisible sandpile to prove Theorem 1.2.3 for 



the classical abelian sandpile model. Section 2.4 is devoted to the proof of Theorem 1.2.1 
for the rotor-router model. 



2.1 Basic Estimate 



Write {Xk)k>o for simple random walk in Z"', and for d> 3 denote by 



giix,y)=E,#{k\Xk = y} 



(2.1) 



the expected number of visits to y by simple random walk started at x. This is the discrete 
harmonic Green's function in Z'^. For fixed x, the function gi{x,-) is harmonic except at 
X, where its discrete Laplacian is —1. Our notation gi is chosen to distinguish between the 
discrete Green's function in Z"' and its continuous counterpart g in M'^. For the definition 



of g, see section 3.1.2 Estimates relating the discrete and continuous Green's functions are 



discussed in section 13.1.51 

In dimension d = 2, simple random walk is recurrent, so the expectation on the 
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right side of (2.1 ) infinite. Here we define the potential kernel 

51 (x, y) = hm y) - x)) (2.2) 

n— »oo 

where 

gnx,y)=E,#{k<n\Xk = y}. 



The hmit defining gi in (2.2) is finite |28| and gi{x,-) is harmonic except at x, where 



its discrete Laplacian is —1. Note that (2.2) is the negative of the usual definition of the 



potential kernel; we have chosen this sign convention so that gi has the same Laplacian in 
dimension two as in higher dimensions. 

Fix a real number m > and consider the function on Z"^ 

7d(a;) = + 77151(0, x). (2.3) 

Let r be such that m = w^r'^, and let 

ld{x) = 7d(x) - 7d( [rj ei) (2.4) 

where ei is the first standard basis vector in Z^. The function plays a central role in our 
analysis. To see where it comes from, recall the divisible sandpile odometer function 

u{x) = total mass emitted from x. 

Let Djn C Z'^ be the domain of fully occupied sites for the divisible sandpile formed from a 
pile of mass m at the origin. For x £ D^, since each neighbor y of x emits an equal amount 
of mass to each of its 2d neighbors, we have 



^u{x) = ^^u{y) -u{a 



2d 

yr^x 

= mass received by x — mass emitted by x 
= 1 - m5ox- 

Moreover, u = on dDm, where for A C Ij'^ we write 

dA = {x £ A'^ \ x ^ y for some y € A} 

for the boundary of A. By construction, the function obeys the same Laplacian condition: 
A'jd = 1 — w-^o! and as we will see shortly, 7^ ~ on dBr- S incG WG expect the doina/in Drn 
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to be close to the ball Bj., we should expect that u k, ^d- In fact, we will first show that u 
is close to 7(i, and then use this to conclude that Dm is close to B^. 

We will use the following estimates for the Green's function |20| H5] : see also 
Theorems 1.5.4 and 1.6.2]. 



9i{x,y) 



-f log|x-y| +K + 0(|x-y| 2), d = 2 
ad\x-y\^-'^ + 0{\x-y\-'^), d>3. 



(2.5) 



Here 



where lu^ is the volume of the unit ball in M.'^, and k is a constant whose 



(d-2)a;d' 

value we will not need to know. Here and throughout this thesis, constants in error terms 
denoted O(-) depend only on d. 

We will need an estimate for the function jd near the boundary of the ball Bj.. We 



first consider dimension d = 2. From (2.5) we have 



where 



(xj = \x\ log \x\. 

TT 



(2.6) 



In the Taylor expansion of (j) about |x| = r 

(t>{x) = 0(r) - </.'(r)(r - \x\) + ^(/>"(t)(r - |:e|)^ 
the linear term vanishes, leaving 

72(x) = (1 + ^) {r - \x\f + 0{m\x\-'^) 
for some t between \x\ and r. 



(2.7) 



(2i 



In dimensions (i > 3, from (2.5) we have 

7d(x) = |xp + adm\x\^~'^ + 0{m\x\~'^). 



Setting (t){x) = \x\'^ + a^mlxp the linear term in the Taylor expansion (2.7) of </> about 
|x| = r again vanishes, yielding 

7^(x) = [l + {d- l){r/tY^ (r - \x\f + 0{m\x\-'^) 



for t between \x\ and r. Together with (2.8), this yields the following estimates in all 
dimensions d> 2. 
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Lemma 2.1.1. Let 7^ be given by (2.4)- For all x we have 



ld{x) > {r-\x\y + 



(2.9) 



Lemma 2.1.2. Let 7^ be given by (2.4)- Then uniformly in r, 



7d(x) = 0(1), x £ Br+i - Br-l. 



The following lemma is useful for x near the origin, where the error term in (2.9) 



blows up. 



Lemma 2.1.3. Let 7^ be given by (2. 4)- Then for sufficiently large r, we have 

Proof. Since 7d(x) — |xp is superharmonic, it attains its minimum in By./^ at a point z on 
the boundary. Thus for any x G B^j^ 

ld{x) - |xp > 7d(z) - 



hence by Lemma 2.1.1 



7d(2;) > (2r/3)2-(r/3)2 + 0(l) > -. □ 



Lemmas 2.1.1 and 2.1.3 together imply the following. 



Lemma 2.1.4. Let be given by (2. 4)- There is a constant a depending only on d, such 
that 'jd > —a everywhere. 



2.2 Divisible Sandpile 

2.2.1 Abelian Property 

In this section we prove the abelian property of the divisible sandpile mentioned 
in the introduction. We work in continuous time. Fix r > 0, and let T : [0, r] ^ Z'^ be a 
function having only finitely many discontinuities in the interval [0, t] for every t < t. The 
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value T(t) represents the site being toppled at time t. The odometer function at time t is 
given by 

ut{x) = C{T-\x)n[0,t]), 

where C denotes Lebesgue measure. We will say that T is a legal toppling function for an 
initial configuration if for every < t < r 

'^mt)) > 1 

where 

Mx) = Mx) + Aut{x) (2.10) 

is the amount of mass present at x at time t. If in addition z^,- < 1, we say that T is 
complete. The abelian property can now be stated as follows. 

Lemma 2.2.1. IfTi : [0,ri] andT2 : [0, T2\ — > Z'^ are complete legal toppling functions 

for an initial configuration uq, then for any x € Z"' 

C {T,-\x)) = C {T,\x)) . 

In particular, ri = T2 and the final configurations , identical. 

Proof. For i = 1,2 write 

ui{x) = C{T-\x)n[0,t]); 
i^iix) = iyo{x) + Auiix). 

Write ti* = ui_ and z^* = z^*.. Let t(0) = and let t{l) < t{2) < ... be the points of 
discontinuity for Ti. Let Xk be the value of Ti on the interval {t(k — l),t(k)). We will show 
by induction on k that 

u'^ixk) >ul^i.-^{xk). (2.11) 
Note that for any x ^ Xk, if u^^f^-^{x) > 0, then letting j < k he maximal such that xj = x, 



since Ti ^ x on the interval (t(j), t(k)), it follows from the inductive hypothesis (2.11 ) that 

u\x) = u\xj) > ul^^^ixj) = n,%)(x). (2.12) 
Since Ti is legal and T2 is complete, we have 
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hence 

Since Ti is constant on the interval (t(k — l),t(k)) we obtain 



By (2.12), each term in the sum on the right side is nonnegative, completing the inductive 
step. 



Since t{k) | ri as A; | oo, the right side of (2.12) converges to u^{x) as | oo 



hence > u^. After interchanging the roles of Ti and T2, the result follows. □ 

2.2.2 Proof of Theorem [1X2] 

Recall that a function s on Z'^ is superharmonic if As < 0. Given a function 7 on 
Z*^ the least superharmonic majorant of 7 is the function 

s{x) = inf{/(x) I / is superharmonic and / > 7}. 

Note that if / is superharmonic and / > 7 then 

Taking the infimum on the left side we obtain that s is superharmonic. 

Lemma 2.2.2. Let T : [0,t] ^ he a complete legal toppling function for the initial 
configuration vq, and let 

u{x) = £{T-^{x)) 

be the corresponding odometer function for the divisible sandpile. Then u = s + j, where 

7(x) = |x|^+ ^ gi{x,y)uo{y) 

and s is the least superharmonic majorant of —7. 



Proof. From (2.10) we have 

Since A7 = 1 — i/q, the difference u — 7 is superharmonic. As u is nonnegative, it follows 
that M — 7 > s. For the reverse inequality, note that s + 7 — u is superharmonic on the 
domain D = {x \ frix) = 1} of fully occupied sites and is nonnegative outside D, hence 
nonnegative inside D as well. □ 
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We now turn to the case of a point source mass m started at the origin: uq = m6o- 
More general starting densities are treated in Chapter |3j In the case of a point source of 
mass m, the natural question is to identify the shape of the resulting domain of fully 



occupied sites, i.e. sites x for which fi{x) = 1. According to Theorem 2.2.3, Dm is extremely 
close to a ball of volume m; in fact, the error in the radius is a constant independent of m. 
As before, for r > we write 

Br = {x£Z'^ : \x\ < r} 
for the lattice ball of radius r centered at the origin. 

Theorem 2.2.3. For m > let Dm C Z"^ be the domain of fully occupied sites for the 
divisible sandpile formed from a pile of size m at the origin. There exist constants c, c' 
depending only on d, such that 

Br—c C Dm C -Bf^c') 

where r = (m/u;^)^/"' and ujd is the volume of the unit ball in W^. 



The idea of the proof is to use Lemma 2.2.2 along with the basic estimates on 7 



Lemmas 2.1.1 and 2.1.2[ to obtain estimates on the odometer function 

u{x) = total mass emitted from x. 
We will need the following simple observation. 

Lemma 2.2.4. For every point x £ Dm — {0} there is a path a; = xo ~ a^i ~ • • . ~ a^fc = o 
in Dm with > u{xi) + 1. 

Proof. If Xi G Dm — {o}, let Xj+i be a neighbor of Xi maximizing Then Xj+i G Dm 

and 

= u{xi) + Au{xi) 
= u{xi) + 1, 

where in the last step we have used the fact that Xi € Dm- D 
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Proof of Theorem 2.2.3 We first treat the inner estimate. Let 7^ be given by (2.4). By 



Lemma 2.2.2 the function n — 7^ is superharmonic, so its minimum in the ball is attained 



on the boundary. Since u > 0, we have by Lemma 2.1.2 



u{x) - 'jdix) > 0(1), X G dBr. 



Hence by Lemma 2.1.1 



n(x) > (r -\x\f +0{r'^/\x\'^), 



X £ Br. 



(2.13) 



2.1.3 



we have 



It follows that there is a constant c, depending only on d, such that u{x) > whenever 
r/3 < |a;| < r — c. Thus Br-c — B^/s C D^. For x G B^./^, by Lemma 
u{x) > + 0(1) > 0, hence Bj./3 C Dm. 

we 



For the outer estimate, note that u — 7^ is harmonic on Dm. By Lemma 2.1.4 



have 7(i > — a everywhere, where a depends only on d. Since u vanishes on dDm it follows 



that li — 7(i < a on Dm. Now for any x G Dm with r — 1 < |rE| < r, we have by Lemma 2.1.2 

u{x) < 7d(x) + a < c' 



for a constant c' depending only on d. Lemma 2.2.4 now implies that Dm C Br+c'+i- D 



2.3 Classical Sandpile 

We consider a generalization of the classical abelian sandpile, proposed by Fey and 
Redig [18^ . Each site in Z"^ begins with a "hole" of depth H. Thus, each site absorbs the first 
H grains it receives, and thereafter functions normally, toppling once for each additional 
2d grains it receives. If H is negative, we can interpret this as saying that every site starts 
with h = —H grains of sand already present. Aggregation is only well-defined in the regime 
h < 2d — 2, since for h = 2d — 1 the addition of a single grain already causes every site in 
Z'^ to topple infinitely often. 

Let Sn,H be the set of sites that are visited if n particles start at the origin in X'^. 
Fey and Redig \18\ Theorem 4.7] prove that 

lim limsup — # {Sn,H A B^-i/d^) = 0, 

where n = uj^r'^, and A denotes symmetric difference. The following theorem strengthens 
this result. 
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Theorem 2.3.1. Fix an integer H > 2 — 2d. Let 5„ = Sn,H be the set of sites that are 
visited by the classical abelian sandpile model in , starting from n particles at the origin, 
if every lattice site begins with a hole of depth H. Write n = io^r'^. Then 

where 

ci = {2d - 1 + H)-^/'^ 

and C2 is a constant depending only on d. Moreover if H > \ — d, then for any e > we 
have 

where 

c[ = id-e + H)-^/^ 
and c'2 is independent of n but may depend on d, H and e. 

Note that the ratio ci/c'^ | 1 as | 00. Thus, the classical abelian sandpile 
run from an initial state in which each lattice site starts with a deep hole yields a shape 
very close to a ball. Intuitively, one can think of the classical sandpile with deep holes 



as approximating the divisible sandpile, whose limiting shape is a ball by Theorem 2.2.3 



Following this intuition, we can adapt the proof of Theorem 2.2.3 to prove Theorem 2.3.1[ 
just one additional averaging trick is needed, which we explain below. 
Consider the odometer function for the abelian sandpile 

u{x) = total number of grains emitted from x. 

Let Tn = {x\u[x) > 0} be the set of sites which topple at least once. Then 

r„ c s„ c T„ u ar„. 

In the final state, each site which has toppled retains between and 2d—l grains, in addition 
to the H that it absorbed. Hence 

H < Au{x) + nSox <2d-l + H, x eTj,. (2.14) 

We can improve the lower bound by averaging over a small box. For x G let 

Q{x,k) = {yeZ'^ : \\x - y\\oo < k} 
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be the box of side length 2k + 1 centered at x, and let 



Write 



= {x\Qix,k) CTn}. 



Le Borgne and Rossin [3T] observe that if T is a set of sites all of which topple, the number 
of grains remaining in T is at least the number of edges internal to T: indeed, for each 
internal edge, the endpoint that topples last sends the other a grain which never moves 
again. Since the box Q{x, k) has 2dk{2k + 1)°'"^ internal edges, we have 

2k 



Au^^\x) >{d + H) 



2k + r ^ " ■ 



(2.15) 



The following lemma is analogous to Lemma 2.2.4 



Lemma 2.3.2. For every point x G Tn adjacent to dTn there is a path x = ~ xi ~ . . . 
Xm = o in Tn with u{xi+i) > u{xi) + 1. 



Proof. By (2.14) we have 



^ u{y) > u{x^). 



Since u(xj_i) < u{xi), some term u{y) in the sum above must exceed u{xi). Let Xj+i = 
y. □ 



Proof of Theorem \2.3.1\ Let 

ldi.x) = {2d-l + H)\x\^ + ngi[o,x), 



and let 



Taking m = n/{2d — 1 + H) in Lemma 2.1.2 we have 



u{x) - Ux) > -Ux) = 0(1), X G dB, 



cir- 



(2.16) 



By (2.14), n — is superharmonic, so (2.16) holds in all of -Bcir- Hence by Lemma 2.1.1 



(2.17) 
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It follows that u is positive on Bc^r-C2 ~ -^cir/a for & suitable constant C2. For x S i^d-r/S) 
by Lemma 2.1.3 we have u{x) > cjr^/A + 0(1) > 0. Thus B^^r^c^ C r„ C Sn- 



For the outer estimate, let 

i^dix) = (d - e + H)\x\'^ + ngi{o,x). 
Choose k large enough so that {d + H) > d — e + H and define 

'>Pd{x) 



(2fc + i)-'^ ^ Vid(y). 

j/eQ(x,fc) 



Finally, let 



2.1.4 



ipdix) = il^dix) - ipd{[cir\ei). 

By (2.15), u^^") — tpd is subharmonic on Tn^\ Taking m = n/{d — e + H) in Lemma 
there is a constant a such that il)d > —a everywhere. Since u^^'^ < {2d + on 
OT^''^ it follows that u^'^) - < a + (2d + H)('^+^^^ on T^''^ Now for any x & Sn with 
c'j^r — 1 < |x| < c'^r we have by Lemma 



2.1.2 



(x) < Vd(2;) +a + {2d + Hf'^+^^^ < C2 



for a constant C2 depending only on d, if and e. Then u{x) < := (2A; + 1)'^C2. Lemma 
now implies that T„ C B^/^^+c'^: hence 



2.3.2 



Sn CTnU dTn C i?c'^^.+c^+l • □ 



We remark that the crude bound of {2d + used in the proof of the outer 

estimate can be improved to a bound of order k"^!!, and the final factor of {2k + l)'^ can 
be replaced by a constant factor independent of k and H, using the fact that a nonneg- 
ative function on Z'^ with bounded Laplacian cannot grow faster than quadratically; see 
Lemma I3.1.18[ 



2.4 Rotor-Router Model 

Given a function / on X'^, for a directed edge {x, y) write 

V/(x,2/) =/(y)-/(x). 
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Given a function s on directed edges in Z'^, write 

divs(x) = w;j^s{x,y). 



2d 

yr^x 

The discrete Laplacian of / is then given by 

A/(x) = divV/=^^/(y)-/(:E). 

2.4.1 Inner Estimate 

Fixing n > 1, consider the odometer function for rotor- router aggregation 

u{x) = total number exits from x by the first n particles. 

We learned the idea of using the odometer function to study the rotor-router shape from 
Matt Cook [n]. 

Lemma 2.4.1. For a directed edge {x,y) in , denote by K{x,y) the net number of cross- 
ings from X to y performed by the first n particles in rotor-router aggregation. Then 

Vu{x,y) = -2dK{x,y) + R{x,y) (2.18) 

for some edge function R which satisfies 

\R{x,y)\ <M-2 

for all edges (x, y). 

Remark. In the more general setting of rotor stacks of bounded discrepancy, the Ad — 2 will 
be replaced by a different constant here. 

Proof. Writing N(x, y) for the number of particles routed from x to y, we have 

u(x)-2d+l , u(x)+2d-l 

hence 

\Vu{x,y) -\-2dK{x,y)\ = \u{y) - u{x) + 2dN{x,y) -2dN{y,x)\ 
<M-2. □ 
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For the remainder of this section Co, Ci, . . . denote constants depending only on d. 
Lemma 2.4.2. Let Q d'L'^ — {o} be a finite set. Then 

^\y\^'^ < CoDiam(O). 

Proof. Let 

Sk = {y ^'L'^ : k<\y\ <k + 1}. 

Then 

^ |y|'-'<A:i-'#5fc<Co. 

Since J7 can intersect at most Diam(O) + 1 distinct sets Sk, the proof is complete. □ 

Lemma 2.4.3. Let G = Gb,. be the Green's function for simple random walk in Z*^ stopped 
on exiting Bf.. Let x G By with > r — p. Then 

<Ci/,. (2.19) 

k-y|<3p 

Proof. Let {Xt)t>Q denote simple random walk in Z"^, and let T be the first exit time from 
Bj.. For fixed y, the function 

A{x) = gi{x, y) - E,5i(Xt, y) (2.20) 
has Laplacian /S.A{x) = —6xy in Bj. and vanishes on dB^, hence A(x) = G{x,y). Let 



x,y Bj. and z ^ y. From (2.5 ) we have 



G2 

\9i{x,y) - gi{x,z)\ < _ , y.zi^x. 
\x y\ 



Using the triangle inequality together with (2.20), we obtain 



|G(x,y) -G{x,z)\ < \gi{x,y) -gi{x,z)\ + Ex\gi{XT,y) - gi{XT,z)\ 



x-yr ^ ^„ - ' \w -y\'^ ^' 



where Hx{w) = Fx{Xt = w). 

Write D = {y £ Br : \x - y\ < 3p}. Then 

^ ^ \G{x, y) - Gix, z)\ < Gsp + C2 H^iw) ^ ^ - y\'-''. (2.21) 

yeD^^y wedBr yen 

y^x ^7^^ 
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hi 



X 



Ho 



Ho 



Figure 2.1: Diagram for the Proof of Lemma 



2.4.4 



By Lemma 2.4.2 the inner sum on the right is at most CoDiam(D) = QCop, so the right 



side of (2.21 ) is bounded above by Cip for a suitable Ci. 



Finally, the terms in which y or z coincides with x make a negligible contribution 



to the sum in (2.19), since for y x Z 



\G{x,x) -G{x,y)\ < \gi{x,x) - gi{x,y)\ + E^\gi{XT,x) - gi{XT,y)\ < C4. □ 



Lemma 2.4.4. Let Hi,H2 be linear half-spaces in TL^ , not necessarily parallel to the coor- 
dinate axes. Let Ti be the first hitting time of Hi. If x ^ Hi U H2, then 

^ ' - 2 h2 V 2/127 

where hi is the distance from x to Hi. 

Proof. If one of Hi,H2 contains the other, the result is vacuous. Otherwise, let Hi be the 
half-space shifted parallel to Hf by distance 2/i2 in the direction of x, and let Tj be the 
first hitting time of Hi U Hi. Let {Xt)t>o denote simple random walk in Z'^, and write 
Mt for the (signed) distance from Xt to the hyperplane defining the boundary of Hi, with 
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Mq = hi. Then Mf is a martingale with bounded increments. Since E^Ti < oo, we obtain 
from optional stopping 



hence 



{X^^ G i^i) < 



/ii + l 
2/io 



Likewise, dM^^ — t is a martingale with bounded increments, giving 



-i 1 



< d(2/i2 + 1)^ Px G ^1 
<d(/li + l) (2/12 + 1) (1 + 



2/i, 



(2.22) 



(2.23) 



Let T = min(Ti,T2). Denoting by Dt the distance from Xt to the hyperplane 
defining the boundary of H2, the quantity 

N, = ^{D^ + {2h,-D,f)-t 

is a martingale. Writing p = ¥x{T = T2) we have 

= ENo = ENt > P^(2/i2)^ + {l-p)dhl - E^T 
> (1 +p)dhl-ExT 



hence by (2.23) 



E^T hi + 1 
P<-^<2^— 
a/i2 h2 



1 + 



2/22 



Finally by (2.221 



F(T,>n)<p + F(Xf^eH,)<l'^[l + ±-^\ □ 



Lemma 2.4.5. Let x € and let p = r + 1 — \x\. Let 

S*k = {yG Br : 2V <\x-y\< (2.24) 
Lei Tfc 6e i/ie /zrst hitting time of S^, and T the first exit time from Br- Then 

^An <T)< C22-\ 
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Figure 2.2: Diagram for the proof of Lemma 



2.4.5 



Proof. Let H be the outer half-space tangent to at the point z E dBr closest to x. Let 
Q be the cube of side length 2^ pl\fd centered at x. Then Q is disjoint from 5^, hence 

where Tqq and T/^ are the first hitting times of dQ and H. Let f^i, . . . , }i2d be the outer 



half-spaces defining the faces of Q, so that Q = n . . . H ^2^. By Lemma 2.4.4 we have 

2d 



1=1 



2d 



< 



5 ^ dist(3;,i?) + 1 



1 + 



1 



Since dist(j;, H) 
20(i'^/^(l -I- completes the proof. 



2^ dist{x, Hi) V 2dist(x,Fi)_ 
2; — -z| < /O and dist(x,ifj) = 2''^^p/Vd, and p > 1, taking C2 



□ 



Lemma 2.4.6. Let G = be the Green's function for random walk stopped on exiting 
Br- Let X & Br and let p = r — \x\. Then 

<C3/9log-. 
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Proof. Let St be given by (|2.24l) , and let 



W = {w £ d{S^ U dSk) : \w-x\<2^p}. 
Let T\Y be the first hitting time of W and T the first exit time from Br- For w £ W let 



Hx{w) = F^{Xr^yAT = w). 

Fixing y G 5^ and z ^ simple random walk started at x must hit W before hitting either 
y or z, hence 

\G{x,y)-Gix,z)\ < H,{w)\G{w,y)-G{w,z)\. 



For any w £ W we have \w\ > r — (2^ + l)p and |y — it;| < 3 ■ 2'^/?. Lemma 2.4.3 yields 



yest zr^y 



By Lemma 2.4.5 we have Ylw&w — so the above sum is at most 2C1C2P. 



Since the union of shells Sq , Si , . . . , S'^^^^^^^ ^ ^^-^ covers all of Br except for those points 
within distance p of x, and J2\x-y\<p'^z'^y \Gix,y) — G{x,z)\ < Cip by Lemma 
result follows. 



2.4.3 



the 

□ 



Proof of Theorem 1.2.1, Inner Estimate. Let k and R be defined as in Lemma [2.4. 1[ Since 
the net number of particles to enter a site x 7^ o is at most one, we have 2d div > — 1. 



Likewise 2d div k(o) = n — 1. Taking the divergence in (2.18), we obtain 



Au(x) < 1 + div R{x), X ^ 0; 
Au{o) = 1 - n + div R{o). 



(2.25) 
(2.26) 



Let T be the first exit time from Br, and define 



fix) = E^u(Xt) - E^T + n #{j < T\Xj = 0}. 

Then Af{x) = 1 for x G Br — {0} and A/(o) = 1 — n. Moreover / > on dBr- It follows 
from Lemma [2.1.2 with m = n that / > 7 — C4 on Br for a suitable constant C4. 
We have 

n(x) - Exu{Xt) = E E^. (^u(XfcAT) - u{^(k+i)AT)) ■ 

k>0 
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Each summand on the right side is zero on the event {T < k}, hence 

E^(^u{XkAT) - W.(-'^(fc+l)AT) I -^/cAt) = -Au{Xk)l{T>k}- 

Taking expectations and using ( |2.25 ) and (2.26), we obtain 

u{x) - EMXt) > 5^ Ex [l{T>k}{nl{x,=o} - 1 - div R{Xk)) 

= n #{k < T\Xk = o}- E,r - E, [l{r>fc}div R{Xk] 

k>0 



hence 



k>0 



^{T>k} ^ R{Xk,. 



(2.27) 



Since random walk exits with probabihty at least ^ every time it reaches a site adjacent 
to the boundary dBr, the expected time spent adjacent to the boundary before time T is 



at most 2d. Since \R\ < Ad, the terms in (2.27) with z £ dBr contribute at most IQd"^ to 
the sum. Thus 



k>0 



y,z£Br 



Sd^. 



For y G Br we have {X^ = y} n{T > k} = {X^^x = 2/}; hence 

u{x)- f{x)>-^Yl E ^'AXkAT = y)R{y,z)-8d^ 



(2.28) 



k>0 y,zeBr 
yr^z 



Write Pkiy) = ^xiXkAT = y)- Note that since V/ and k are antisymmetric, R is antisym- 
metric. Thus 



Y Pk{y)R{y,z) = - J2 Pk{z)R{y,z) 



y,z£Br 

y^z 



y,zeBr 
y^z 



E ''^'^-''^''^ R^y.z). 



y,z&Br 

y^z 



Summing over k and using the fact that \R\ < Ad, we conclude from (2.28) that 



i(x)>/(x)- Y |G(x,y)-G(x,z)|-8d2, 



y,zeBr 
yr^z 
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where G = is the Green's function for simple random walk stopped on exiting Br- By 
Lemma 12.4.61 we obtain 



u{x) > f{x) — C^{r — |x|) log — r — M' 

r — \x\ 



Using the fact that / > 7 — C4, we obtain from Lemma 2.1.1 



uix) > (r — Ixl)^ — Csfr — Ixl) log — '—j—r + O [ , ,, 

r — \x\ V ud" 



2.1.3 



The right side is positive provided r/3 < |x| < r — Cslogr. For x G -6^/3) by Lemma 
we have u{x) > — Csr log | > 0, hence -Br-Cslogr C A„. □ 

2.4.2 Outer Estimate 

The following result is due to Holroyd and Propp (unpublished); we include a 



proof for the sake of completeness. Notice that the bound in (2.29) does not depend on the 
number of particles. 

Proposition 2.4.7. LetT be a finite connected graph, and letY <Z Z he subsets of the vertex 
set ofT. Let s be a nonnegative integer-valued function on the vertices ofV. Let Hw{s,Y) 
be the expected number of particles stopping in Y if s{x) particles start at each vertex x 
and perform independent simple random walks stopped on first hitting Z. Let Hr{s,Y) be 
the number of particles stopping in Y if s{x) particles start at each vertex x and perform 
rotor-router walks stopped on first hitting Z . Let H{x) = Hw{lx, Y). Then 

\Hr{s,Y) - HUs,Y)\ < ^^|/7(u)-F(t;)| (2.29) 
independent of s and the initial positions of the rotors. 

Proof For each vertex u, arbitrarily choose a neighbor r]{u). Order the neighbors rj{u) = 
vi,V2, ■ ■ ■ ,Vd of u SO that the rotor at u points to fj+i immediately after pointing to Vi 
(indices mod d). We assign weight w{u,rj(u)) = to a rotor pointing from u to r](u), 
and weight w{u,Vi) = H{u) — H{vi) + uj{u,Vi^i) to a rotor pointing from u to Vi. These 
assignments are consistent since IL \s a harmonic function: ^i{H{u) — H{vi)) = 0. We 
also assign weight H{u) to a particle located at u. The sum of rotor and particle weights 
in any configuration is invariant under the operation of routing a particle and rotating the 
corresponding rotor. Initially, the sum of all particle weights is H^^s, Y). After all particles 
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Figure 2.3: Diagram for the proof of Lemma 



2.4.8 



have stopped, the sum of the particle weights is Hr{s,Y). Their difference is thus at most 
the change in rotor weights, which is bounded above by the sum in (2.29). □ 

For p G Z let 



Then 



{xeZ'^ : p< \x\ < p + 1}. 



Bp = {x£Z'^ : |x| < p} = cSo U . . . U Sp^i. 



(2.30) 



Note that for simple random walk started in Bp, the first exit time of Bp and first hitting 
time of Sp coincide. Our next result is a modification of Lemma 5(b) of [29]. 

Lemma 2.4.8. Fix p > 1 and y G Sp. For x & Bp let H{x) = Pi.(Xr = y), where T is the 
first hitting time of Sp. Then 

(2.31) 



\x-y\'^-^ 
for a constant J depending only on d. 

Proof. We induct on the distance |rE — y|, assuming the result holds for all x' with |x' — y| < 
^|rc— y|; the base case can be made trivial by choosing J sufficiently large. By Lemma 5(b) of 
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|29j . we can choose J large enough so that the result holds provided \y\ — \x\ > 2~'^~^\x — y\. 
Otherwise, let Hi be the outer half-space tangent to Sp at the point of Sp closest to x, and 
let H2 be the inner half-space tangent to the ball S of radius — y| about y, at the point 



of S closest to X. By Lemma 2.4.4 applied to these half-spaces, the probability that random 
walk started at x reaches S before hitting Sp is at most 2^^^. Writing T for the first hitting 



time of U 5p, we have 



H{x) < V = x')H{x') < 2^-'^ J ■ 



l-d 



/ ^ ^ xy^^'j' ^ }-^-^ V 7 — \ 2 

where we have used the inductive hypothesis to bound H{x'). □ 

The lazy random walk in Z'^ stays in place with probability 2, and moves to each 
of the 2d neighbors with probability We will need the following standard result, which 
can be derived e.g. from the estimates in [36', section 11.12; we include a proof for the sake 
of completeness. 

Lemma 2.4.9. Given u ^ v £ Z'^, lazy random walks started at u and v can be coupled 
with probability 1 — C/R before either reaches distance R from u, where C depends only on 
d. 

Proof. Let i be the coordinate such that Ui ^ Vi. To define a step of the coupling, choose 
one of the d coordinates uniformly at random. If the chosen coordinate is different from 
i, let the two walks take the same lazy step so that they still agree in this coordinate. If 
the chosen coordinate is i, let one walk take a step while the other stays in place. With 
probability ^ the walks will then be coupled. Otherwise, they are located at points u',v' 
with |n' — = 2. Moreover, P(|u — n'| > < ^ for a constant C depending only on d. 

From now on, whenever coordinate i is chosen, let the two walks take lazy steps in opposite 
directions. 

Let 

Rx = lx X,= ^ 

be the hyperplane bisecting the segment [n',t;']. Since the steps of one walk are reflections 
in Hi of the steps of the other, the walks couple when they hit Hi . Let Q be the cube of 
side length Rj\fd -|- 2 centered at n, and let H2 be a hyperplane defining one of the faces 



of Q. By Lemma 2.4.4 with hi = \ and /i2 = R/^^fd, the probability that one of the walks 
exits Q before the walks couple is at most 2d- 1^ (l + < 40d'^/'^{l + 2Vd)^ /R. □ 
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Lemma 2.4.10. With H defined as in Lemma 2.4.8, we have 



Y^\H{u)-H{v)\< J' log p 

u&Bp v^u 



for a constant J' depending only on d. 



Proof. Given u ^ Bp and u ~ u, by Lemma 2.4.9 lazy random walks started at u and 



V can be coupled with probability 1 — 2C /\u — y\ before either reaches distance |m — y\/2 



from u. If the walks reach this distance without coupling, by Lemma 2.4.8 each has still 
has probability at most J/\u — y\'^~^ of exiting Bp at y. By the strong Markov property it 
follows that 



\H{u)-H{v)\ < 

Summing in spherical shells about y, we obtain 

p 



2CJ 

\u - y\' 



t=i 



2C T 



We remark that Lemma 2.4.10 could also be inferred from Lemma 2.4.8 using [28l 



Thm. 1.7.1] in a ball of radius \u — y\/2 about u. 



Proof of Theorem 1.2.1, Outer Estimate. Fix integers p > r, h > 1. In the setting of 



Proposition 2.4.7 let G be the lattice ball Bp-^-h+l, and let Z = Sp^h- Fix y G Sp+h and 
let Y = {y}. For x £ Sp let s{x) be the number of particles stopped at x if all particles in 
rotor-router aggregation are stopped upon reaching Sp. Write 

H{x) = F,{Xt = y) 



where T is the first hitting time of Sp^h.. By Lemma 2.4.8 we have 



xeSp 



(2.32) 



where 

is the number of particles that ever visit the shell Sp in the course of rotor-router aggregation. 
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By Lemma 2.4.10 the sum in (2.29 ) is at most J' log /i, hence from Propositon 



2.4.7 



and (2.32) we have 



Hr{s.y) <^ + J'logh. 



Let p{0) = r, and define p{i) inductively by 

p{i + 1) = min [pii) + iV'(/f mm{p > pii)\Np < iVp(.)/2}} . 
Fixing h < p{i + 1) — we have 

2d-2 

h^-HoghK N^^^' logiV,(i) < iV^(,); 



(2.33) 



(2.34) 



so (2.33) with p = p{i) simplifies to 



Hr{s,y) < 



CN, 



Pii) 



(2.35) 



where C = J + J'. 

Since all particles that visit during rotor-router aggregation must pass 
through 5p(j), we have 

Np^i)+h< Yl Hr{s,y). (2.36) 

Let Mfc = #{An n Sk)- There are at most Mp(j)_(_^ nonzero terms in the sum on the right 
side of ( 2.36| ), and each term is bounded above by (2.35), hence 



> 



h 



d-l 



where the second inequality follows from iVp(j)-|-/i > Summing over h, we obtain 

(2.37) 



p(*+i)-i , 

E M,>—{p{i + i)-p{i)-ir 
p=p{i)+i 



The left side is at most -/Vp(j), hence 

p{i + 1) - p{i) < (2<iCiVp(,))i/^ < iv;/f 



provided iVp(,) > C := {2dCf'^-^. Thus the minimum in (^2.34^ is not attained by its first 
argument. It follows that Np(^^_^^l^ < hence Np(^^i^^^-^ < C for a sufficiently large 

constant a. 
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By the inner estimate, since the ball Br~c\ogr is entirely occupied, we have 
^ Mp < Wrf/ - uJdir - clog r)'^ 

p>r 

< cduj^r'^"^ log r. 



Write Xi = p{i + 1) — p{i) — 1; by (2.37) we have 



a log r 

^ xf < cduJdr'^-Hogr, 

By Jensen's inequality, subject to this constraint, '^Xi is maximized when all Xi are equal, 
in which case Xi < C"r^^^^'^ and 

p{a log r) = r + ^ < r + C"r^~^/'^ log r. (2.38) 

Since Np^^^i^gj.^ < C' we have -/Vp(aiogj,)_|_c/ = 0; that is, no particles reach the shell 



'Sp(aiogr)+C' ■ Taking c' = C" + C", we obtain from (2.38) 



log r) ■ 



□ 
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Chapter 3 



Scaling Limits for General Sources 



This chapter is devoted to proving Theorems 1.1.2 1.1.3 and 1.3.1 The proofs use 
many ideas from potential theory, and the relevant background is developed in section |3.1[ 



The proof of Theorem 1.1.2 is broken into three sections, one for each of the three aggrega- 
tion models. Of the three models the divisible sandpile is the most straightforward and is 
treated in section EiH The rotor-router model and internal DLA are treated in sections 13^ 



and 3.4 respectively. Finally, in section 3.5 we deduce Theorem 1.3.1 for multiple point 



sources from Theorem [1.1.3 along with our results for single point sources. 



3.1 Potential Theory Background 



In this section we review the basic properties of superharmonic potentials and of 
the least superharmonic majorant. For more background on potential theory in general, we 
refer the reader to [3 ; for the obstacle problem in particular, see [3 [12] . 



3.1.1 Least Superharmonic Majorant 

Since we will often be working with functions on R*^ which may not be twice 
differentiable, it is desirable to define superharmonicity without making reference to the 
Laplacian. Instead we use the mean value property. A function u on an open set C is 
superharmonic if it is lower-semicontinuous and for any ball B{x, r) C 0, 

^ ^ ' ■ (3.1) 



u{x) > Aru{x) := I u{y)dy. 

^dr"" JB{x,r) 
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Here Ud is the volume of the unit bah in M'^. We say that u is subharmonic if —u is 
super harmonic, and harmonic if it is both super- and subharmonic. 

The following properties of superharmonic functions are well known; for proofs, 
see e.g. |3], [16] or |35] . 

Lemma 3.1.1. Let u be a superharmonic function on an open set ncR'^. Then 
(i) u attains its minimum in 0, on the boundary, 
(a) If h is harmonic on 0, and h = u on dO,, then u > h. 
(Hi) If B[x,tq) C B{x,ri) C $7, then 

ArgU{x) > Ar^u{x). 

(iv) If u is twice differentiable on 17, then An < on 17. 

(v) If B C 0, is an open ball, and v is a function on 0, which is harmonic on B, continuous 
on B, and agrees with u on B'^, then v is superharmonic. 

Given a function 7 on M'^ which is bounded above, the least superharmonic majo- 
rant of 7 (also called the solution to the obstacle problem with obstacle 7) is the function 

s(x) = inf{/(3;)|/ is continuous, superharmonic and / > 7}. (3.2) 

Note that since 7 is bounded above, the infimum is taken over a nonempty set. 

Lemma 3.1.2. Let ^ be a uniformly continuous function which is bounded above, and let 



s be given by (3.2). Then 
(i) s is superharmonic. 
(a) s is continuous. 
(Hi) s is harmonic on the domain 

D = {xe R'^\s{x) > -fix)}. 
Proof, (i) Let / > 7 be continuous and superharmonic. Then f > s. By the mean value 



property (3.1), we have 

f>Arf> ArS. 
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Taking the infimum over / on the left side, we conclude that s > ArS. 
It remains to show that s is lower-semicontinuous. Let 

w(7,r)= sup h{x)-j{y)\. 

,\x-y\<r 



Since 

ArS > Ar^ > 7 — w(7, r) 
the function ArS + a;(7,r) is continuous, superharmonic, and lies above 7, so 

ArS < s < ArS + w(7, r). 

Since 7 is uniformly continuous, we have 01(7, r) | as r | 0, hence ArS s as r [ 0. 
Moreover if vq < ri, then by Lemma 3.1.1[ iii) 



ArgS = lim Ar^ArS > lim Ar^^ArS = Arj^S. 

Thus s is an increasing limit of continuous functions and hence lower-semicontinuous. 

(ii) Since s is defined as an infimum of continuous functions, it is also upper- 
semicontinuous. 

(iii) Given x £ D, write e = s(x) — 7(x). Choose 6 small enough so that for all 
y€B = B{x,6) 

|7(x) -7(y)| < I and \s{x) - s{y)\ < ^. 

Let / be the continuous function which is harmonic in B and agrees with s outside B. By 
Lemma 3.1.1 V), / is superharmonic. By Lemma 3.1.1 1), f attains its minimum in S at a 



point z G dB, hence iov y £ B 



f{y) > f{z) = s{z) > s{x) - I = 7(x) + I > 7(y). 



It follows that / > 7 everywhere, hence f > s. From Lemma 3.1.1[ ii) we conclude that 



f = s, and hence s is harmonic at x. □ 

3.1.2 Superharmonic Potentials 

Next we describe the particular class of obstacles which relate to the aggregation 
models we are studying. For a bounded measurable function a on M'^ with compact support, 
write 

Ga{x) = I g{x,y)a{y)dy, (3.3) 
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where 



flog|x-y|, d = 2; 



(3.4) 



= < 



a'd\x - y 



2-d 



d>3. 



Here 



2 



where ujd is the volume of the unit ball in W^. Note that (3.4) differs 



{d-2)uja 



by a factor of 2d from the standard harmonic potential in M'^; however, the normalization 
we have chosen is most convenient when working with the discrete Laplacian and random 
walk. 

The following result is standard; see jl6| Theorem 1.1.7.2]. 
Lemma 3.1.3. Let a be a measurable function on M'^ with compact support, 
(i) If a is bounded, then Ga is continuous, 
(a) If a is C^, then Ga is and 



mains true in the sense of distributions. For our applications, however, we will not need 
this fact, and the following lemma will suffice. 

Lemma 3.1.4. Let a be a bounded measurable function on M'^ with compact support. If 
a >0 on an open set i7 C M'^, then Ga is superharmonic on 17. 

Proof. Suppose B(x, r) C 11. Since for any fixed y the function /(x) = g{x, y) is superhar- 
monic in X, we have 



2da. 



(3.5) 



Regarding (ii), if we remove the smoothness assumption on a, equation (3.5 1 re- 




ArGa{x). □ 



By applying Lemma 3.1.4 both to a and to —a, we obtain that Ga is harmonic off the 
support of a. 
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Let B = B{o,r) be the ball of radius r centered at the origin in M'^. We compute 
in dimensions d > 3 



GlBix) 



d-2 I- 
{d-2y\x\ 



d-2 



\x\ < r 



\x\ > r. 



(3.6) 



Likewise in two dimensions 



G1b{x) 



r^(l — 21ogr) — |x| < r 



-2r^ log \x\ 



(3.7) 



\x\ > r. 



Fix a bounded nonnegative function a on with compact support, and let 

7(x) = — |x|^ — Ga{x). 



Let 



(3.8) 
(3.9) 



s{x) = inf{/(x)|/ is continuous, superharmonic and / > 7} 
be the least superharmonic majorant of 7, and let 

D = {xe R'^\s{x) > 7(x)} (3.10) 

be the noncoincidence set. 

Lemma 3.1.5. (i) 7(x) + |xp is subharmonic on M'^. 

(ii) If (T < M on an open set O C , then 7(x) — (M — l)|xp is superharmonic on Q. 
Proof, (i) By Lemma 3.1.4 since a is nonnegative, the function 7(x) + |xp = —Ga{x) is 



subharmonic on 



(ii) Let B = i?(o, R) be a ball containing the support of a. By (3.6) and (3.7), for 
X G B we have 



bp = CdR^ - G1b(x) 



where C2 = 1 — 2 log R and = ^"^^ d > 3. Hence for x G B we have 

7(x) - (M - l)|xp = -Gct(x) - M|x|2 

= G{M1b - (7){x) - CdMR^. 



Since a < MIb on 17, by Lemma 3.1.4 the function 7 — (M — l)|xp is superharmonic 
in B ni^. Since this holds for all sufficiently large R, it follows that 7 — (M — l)|xp is 
superharmonic on all of $7. □ 
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Lemma 3.1.6. Let u = s — ^ , where 7 and s are given by (3.8) and (3.9). Then 
(i) u{x) — |xp is superharmonic on M.'^. 

(a) If a < M on an open set 17 C W^, then u{x) + M|xp is subharmonic on il. 
Proof, (i) By Lemmas 3.1.2[ i) and 3.1.5[ i), the function 

u — = s — (7 + 

is the difference of a superharmonic and a subharmonic function, hence superharmonic on 

ia>d 



(ii) With Ar defined by (3.1), we have 

= I + 2x-y + \yf) dy 

^dr'^ JB(o,r) 



, , (dujdt'^'^)t'^ dt 
^dr'^ Jo 

dr"^ 



d + 2 



By Lemma 3.1.5 we have 



ArS + > ylr.7 + Ar\x\^ - |xp > 7. 



Since A^s is continuous and superharmonic, it fohows that ^4.^5 + > s, hence 

ArS + Ar\x\'^ = ArS + + -J-^ >S + 



Thus s + |xp is subharmonic on M"^, and hence by Lemma 3.1.5 li) the function 

u + M\x\^ = {s + \x\^) -{-i-{M - l)|xp) 

is subharmonic on VI. 



□ 



For A C M"', write A° for the interior of A and A for the closure of A. The 
boundary of A is dA = A — A° . 



Lemma 3.1.7. Let D be given by (3.10). Then {a > 1}° C D. 
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Proof. If o" > 1 in a ball B = B{x, r), by (3.6) and (3.7), for y & B we have 



7(y) = - Ga{y) = -c^ - G{a - lB){y). 



By Lemma 3.1.4 it follows that 7 is subharmonic in B. In particular, s > 7 in 5, so 
x(^D. □ 

The next lemma concerns the monotonicity of our model: by starting with more 
mass, we obtain a larger odometer and a larger noncoincidence set; see also 



Lemma 3.1.8. Let ai < a2 be functions on with compact support, and suppose that 
f^d (T2{x)dx < 00. Let 

7j(a;) = -|x|^ - Go-j(x), i = l,2. 
Let Si be the least superharmonic majorant of ji, let Ui = Si — 'ji, and let 

Di = {x\si{x) > 7i(x)}. 

Then ui < U2 and Di C D2. 
Proof. Let 

S = S2 + G(fT2 - <Tl). 

Then s is continuous and superharmonic, and since S2{x) > — |xp — Ga2{x) we have 

s{x) > — — Gai{x) 

hence s > si. Now 

U2 - Ul = S2 - Sl + G((T2 - (Ti) = S - Si > 0. 



Since Di is the support of Uj, the result follows. 



□ 



Our next lemma shows that we can restrict to a domain 17 C M"' when taking the 
least superharmonic majorant, provided that 0, contains the noncoincidence set. 



Lemma 3.1.9. Let 7, s, D be given by (3.8)-(3.10). Let Q (ZW^ be an open set with D (ZQ. 
Then 

s{x) = inf{/(x)|/ is superharmonic on fi, continuous, and f > 7}. (3.11) 
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Proof. Let / be any continuous function which is superharmonic on $7 and > 7. By 
Lemma 3.1.2[ iii), s is harmonic on D, so f — s is superharmonic on D and attains its 
minimum in D on the boundary. Hence f — s attains its minimum in ^7 at a point x where 
s{x) = 7(3;). Since / > 7 we conclude that f > s on 0, and hence everywhere. Thus s is at 



most the infimum in (3.11). Since the infimum in (3.11) is taken over a strictly larger set 
than that in (3.9), the reverse inequality is trivial. □ 



3.1.3 Boundary Regularity for the Obstacle Problem 



Next we turn to the regularity of the solution to the obstacle problem (3.9) and 



of the free boundary dD. There is a substantial literature on boundary regularity for the 
obstacle problem. In our setup, however, extra care must be taken near points where 
cr(x) = 1: at these points the obstacle (3.8) is harmonic, and the free boundary can be 



badly behaved. We show only the minimal amount of regularity required for the proofs of 
our main theorems. Much stronger regularity results are known in related settings; see, for 
example, [3 [S] . 

The following lemma shows that if the obstacle is sufficiently smooth, then the 
superharmonic majorant cannot separate too quickly from the obstacle near the boundary 
of the noncoincidence set. For the proof, we follow the sketch in Caffarelli [71 Theorem 2]. 



As usual, we write for the inner e-neighborhood of D, given by (1.8) 



Lemma 3.1.10. Let a be a function on M'^ with compact support. Let ^,s,D be given 



by (3.8)-(3.10), and write u = s — 7. Thenu is C^, and for y G dD^^ we have \Vu{y)\ < Cqe, 



for a constant Cq depending on a. 
Proof. Fix xq G dD, and define 

m-- 



+ (V7(xo),2; - xo). 

Let A be the maximum second 



3.1.3 



Since o" is C"*^, we have that 7 is C"^ by Lemma 
partial of 7 in the ball B = B{xQ,4:e). By the mean value theorem and Cauchy-Schwarz 
for x € B we have 



\L{x) - 7(x)| = |(V7(xo) - Vj{x^),x- xo)\ 
< AVd\xQ — x^:\\x — xo\ < ce^, 



(3.12) 



50 



where c = 16 A^/d. Hence 

s{x) > 7(x) > L{x) — ce^, x G B. 

Thus the function 

w = s — L + ce^ 

is nonnegative and superharmonic in B. Write w = wq + wi, where wq is harmonic and 
equal to w on dB. Then since s{xo) = 7(xo), we have 

woi^o) < wixo) = s{xo) - L{xo) + ce^ = ce^. 

By the Harnack inequahty, it follows that wq < c'e^ on the ball B' = B{xo, 2e), for a suitable 
constant c'. 

Since wi is nonnegative and vanishes on dB, it attains its maximum in i? at a 
point xi in the support of its Laplacian. Since Awi = As, by Lemma 3.1.2[ iii) we have 
s{xi) = 7(xi), hence 

'Wi{xi) < w{xi) = s{xi) — L{xi) + ce^ < 2ce^, 



where in the last step we have used (3.12). We conclude that < w < (2c + c')e^ on B' and 
hence |s — L| < (c + c')e'^ on B'. Thus on B' we have 



|n| = |s — 7I < |s — L| + I7 — L| < (2c + c )e 

In particular, u is differentiable at xq, and Vu(xo) = 0. Since s is harmonic in D and equal 
to 7 outside D, it follows that u is differentiable everywhere, and off dD. 

Given y £ dD^, let xq be the closest point in dD to y. Since B{y,e) C D, by 
Lemma 3.1.2[ iii) the function w is harmonic on B{y,e), so by the Cauchy estimate |3l 



Theorem 2.4] we have 



\Vs{y) - V7(xo)| = \Vw{y)\ < - sup w{z). (3.13) 

e z&B{y,e) 



Since B{y, e) C B' , the right side of (3.13) is at most C(2c + c')e, hence 

|Vu(y)| < \Vs{y) - V7(xo)| + |V7(y) - V7(xo)| < Coe 
where Co = C(c + c') + Thus u is on as well. □ 
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Next we show that mass is conserved in our modeh the amount of mass starting 
in D is Jj^ a{x)dx, while the amount of mass ending in D is C{D), the Lebesgue measure 
of D. Since no mass moves outside of we expect these to be equal. Although this seems 
intuitively obvious, the proof takes some work because we have no a priori control of the 
domain D. In particular, we first need to show that the boundary of D cannot have positive 
d-dimensional Lebesgue measure. 

Proposition 3.1.11. Let a be a function on with compact support, such that 



£(o--i(l)) = 0. Let D be given by (3J^. Then 

(i) C{dD) = 0. 

(ii) For any function h € C^{D) which is superharmonic on D, 



h{x)dx < / h{x)(j{x)dx. 
D Jd 

Note that by applying (ii) both to h and to —h, equality holds whenever h is 
harmonic on D. In particular, taking h = 1 yields the conservation of mass: Jj^cr(x)dx = 
C{D). 

The proof of part (i) follows Friedman [19', Ch. 2, Theorem 3.5]. It uses the 
Lebesgue density theorem, as stated in the next lemma. 

Lemma 3.1.12. Let A CM.'^ be a Lebesgue measurable set, and let 
A' = {xeA \ liminf e"'^£(B(x,e) n A'^) > O}. 

Then C{A') = 0. 

To prove the first part of Proposition |3.1.lT given a boundary point x £ dD, the 
idea is first to find a point y in the ball B{x, e) where u = s — 7 is relatively large, and then 



to argue using Lemma 3.1.10 that a ball B{y,ce) must be entirely contained in D. Taking 
A = dD in the Lebesgue density theorem, we obtain that x S A'. 



The proof of the second part of Proposition 3.1.11 uses Green's theorem in the 

form 

[ {uAh - Auh)dx = [ (u^ - ^h] dr. (3.14) 
Jd' JdD' \ dn dn J 

Here D' is the union of boxes x^ that are contained in D, and dx is the volume measure in 

D' , while dr is the {d — l)-dimensional surface measure on dD' . The partial derivatives on 

the right side of (3.14) are in the outward normal direction from dD'. 
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Proof of Propostion 3.1.11\ (i) Fix < A < 1. For x G dD with (t{x) < A, for small enough 

2 

f{y) = i{y) + 



e we have a < on B{x,e). By Lemma 3.1.5 ii) the function 



is superharmonic on B{x, e) Pi D, so by Lemma 3.1.2[ iii) the function 



w{y) = u{y) - - 

= s{y) - fiy) + (1 - \){x,y) - ^^kP 

is subharmonic on B(x, e) H D. Since w{x) = 0, its maximum is not attained on dD, so it 
must be attained on dB{x, e), so there is a point y with |x — y| = e and 

1 - A 



u{y) > 



By Lemma 3.1.10 we have | Vti| < (1 + c)Coe in the ball -B(y, ce). Taking c : 



l-A 
4Co ' 



we obtain for z G -B(y, ce) 

m(z) = n(y) + (Vu(y*), z - y) 

> ^^e^ - c(l + c)Coe2 > 0. 

Thus for any x & dD r\ {a < \} 

£{B{x, (1 + c)e) n {dDf) > ujdicef. 
By the Lebesgue density theorem it follows that 

C{dDr\{a < A}) = 0. 



(3.15) 



By Lemma 3.1.7 we have o" < 1 on dD. Taking A | 1 in (3.15), we obtain 



C{dD) < C{dD n {cj < 1}) + C{a-'{1)) = 



(ii) Fix 5 > and let 

D' = IJ X 
where = x + [— f, f]"^. By Lemmas 



□ 



3.1.2 iii) and 



3.1.3 Ii), in D' we have 



Aii = As - A7 = 2d(l - a). 
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Now by Green's theorem (3.14), since u is nonnegative and h is superharmonic, 
/ (1 — ajhdx = — - / Ann ax < — - / u- 7^/1 dr 



(3.16) 



where n denotes the unit outward normal vector to dD'. By Lemma 3.1.10 the integral on 
the right side is bounded by 



dD' 



u 



dh 
d~ 



n 



+ 



du , 
d: 



-h 



n 



dr < CoVdi6^\\Vh\\^ + 6\\h\\^)Cd-i{dD'), (3.17) 



where Cd-i denotes {d — l)-dimensional Lebesgue measure. Let 

s= u -° 

Since C{dD) = 0, given e > we can choose 5 small enough so that C{S) < e. Since 
dD' C dS, we have 



Cd-i[dD') < Cd^i{dS) < 



2dC{S) 



Since D is open, C{D') ] C{D) as (5 J, 0. Taking 5 < e smaller if necessary so that 
C{D') > C{D) - e, we obtain from ( |3l6| ) and ( [sTtI ) 



/ (1 - a)hdx < (M + l)e + CoVd{\\h\\oo + \\Vh\\^)e 
Jd 

where M is the maximum of |fT|. Since this holds for any e > 0, we conclude that 
h{x)dx < h{x)a{x)dx. □ 



We will need a version of Proposition 3.1.11 which does not assume that a is or 
even continuous. We can replace the assumption and the condition that £(cr^^(l)) = 
by the following condition. 



For all X G M*^ either cr{x) < A or ct{x) > 1 



for a constant A < 1. Then we have the following result. 



(3.18) 



Proposition 3.1.13. Let a be a bounded function on M'^ with compact support. Let D be 



given by (3.10), and let D = D U {a > 1}°. If a is continuous almost everywhere and 
satisfies (S.IS ), then 



(i) C{dD) = 0. 
(a) C.(yD^ = f^a{x)dx. 
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(in) C[D) = J^a{x)dx. 

In particular, taking a = ly^ + 1^, we have 

A® B = AU BU D = D 



where © denotes the smash sum (1.6). From (iii) we obtain the volume additivity of the 
smash sum. 

Corollary 3.1.14. Let A,B C be bounded open sets whose boundaries have measure 
zero. Then 

C{A®B) = C{A)+C{B). 



Proof of Proposition 3.1.13\ (i) Fix e > 0. Since a is continuous almost everywhere, there 
exist functions (Tq < a < ai with f^d{(^i — o'o)dx < e. Scaling by a factor of 1 + 5 
for sufficiently small 6, we can ensure that C{a^^{l)) = 0, so that (Tq and ai satisfy the 



hypotheses of Proposition 3.1.11 For i = 0, 1 let 

7i(x) = - Gai{x), 

and let Sj be the least superharmonic major ant of 7j. Choose a with A < a < 1 such that 
C{a^^{a)) = for i = 0, 1, and write 

A = {si > 7i}; 

Si = D,U{ai > a}. 



By (3.18) we have {(Tq ^ ck} C {a > 1}", hence by Lemma 3.1.8 



So = DqU {(70 > a} C D U {a > 1}° = D 



C D C DiU{a>l} = Si. 



(3.19) 
(3.20) 



For z = 0, 1 write 

1, X G A 

cri(x), X^Di. 
By Proposition 3.1.11| [ii) with = 1, we have 

/ a\ = £(A) + / = / 



(3.21) 
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For < ao < a, we have 



C{Si - So) <C{Sin {a°o < ao})+C{{ao < a°Q < a}). 
Since > a on Si, the first term is bounded by 

£{Si n {ctq < ao}) < - 



(3.22) 



< 



a — ao 



a — ao a — ao 

where the equahty in the middle step fohows from ( |3.21 ). 

By Proposition 3.1.11[ i) we have C{dDi) = 0, hence C{dSi) = 0. Taking oq 



a - we obtain from (3.19), (3.20), and (3.22) 



C{dD) <C{Si-So) =C{Si-So) < + C{{a - < ao < a}). 

Since this holds for any e > 0, the result follows, 
(ii) Write 

1, xe D 



a°{x) = { 



ct(x), X ^ D. 



From (3.21) we have 



Jm.'i JR<* 



fji. 



The left and right side differ by at most e. Since /jjd co ^ Im.<^^ — /igdO'i) and e > i 
arbitrary, it follows that f^a o = jj^^ cr° , hence 



IS 



[ a = [ a° - [ a = C{D). 
Jd Jw Jd'' 



(3.23) 



(iii) By Lemma 3.1.7 if cr{x) > 1 and x ^ D, then a is discontinuous at x. Thus 



(7 = 1 almost everywhere on D — D, and we obtain from (3.23) 



!a = C{D)+ f a = C{D). □ 

JD Jd-D 



Our next lemma describes the domain resulting from starting mass m > 1 on a 
bah in M^. Not surprisingly, the result is another ball, concentric with the original, of m 
times the volume. In particular, if m is an integer, the m-fold smash sum of a ball with 
itself is again a ball. 
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Lemma 3.1.15. Fix m > 1, and let D be given by (3.10) with a = ml^^gj.). Then 
D = B{o,m^/'^r). 

Proof. Since 7(3;) = —\x\'^ — Ga{x) is spherically symmetric and the least superharmonic ma- 
jorant commutes with rotations, D is a ball centered at the origin. By Proposition 3.1.13[ ii) 



we have C{D) = mC{B{o,r)). □ 
Next we show that the noncoincidence set is bounded; see also \4:2\ Cor. 7.2]. 

Lemma 3.1.16. Let a be a function on Mf^ with compact support, satisfying < o" < M . 
Let D be given by j^S.lClj^ . Then D is bounded. 



Proof. Let B = B{o, r) be a ball containing the support of a. By Lemmas 3.1.8 and 3.1.15 



the difference s — 7 is supported in B{o, M^^'^r). □ 

3.1.4 Convergence of Obstacles, Majorants and Domains 

In broad outline, many of our arguments have the following basic structure: 

convergence of densities =^ convergence of obstacles 

=^ convergence of majorants 
=^ convergence of domains. 

An appropriate convergence of starting densities is built into the hypotheses of the theorems. 
From these hypotheses we use Green's function estimates to deduce the relevant conver- 
gence of obstacles. Next, as we have already seen in Lemmas 3.1.5 and 3.1.6[ properties 



of the obstacle can often be parlayed into corresponding properties of the least superhar- 
monic majorant. Finally, deducing convergence of domains (i.e., noncoincidence sets) from 
the convergence of the majorants often requires additional regularity assumptions. The 
following lemma illustrates this basic three-step approach. 

Lemma 3.1.17. Let a and an, n = 1,2, . . . be densities on M*^ satisfying 

0<cr,an< MIb 

for a constant M and ball B = B{o,R). Suppose that a is continuous except on a set of 
Lebesgue measure zero, and that 

an{x) ^ a{x) (3.24) 
as n ^ 00, for all continuity points x of a. Then 
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(i) Gan Ga uniformly on compact subsets ofM.'^. 

(ii) s„ — > s uniformly on compact subsets ofM.'^, where s,s„ are the least superharmonic 
majorants of the functions 7 = — |xp — Ga, and 7„ = —\x\'^ — Gan, respectively. 

(Hi) For any e > we have C ^(n) for o-H sufficiently large n, where D,D(^n) 
noncoincidence sets {s > 7} and {s„ > 7„}, respectively. 



Proof, (i) Since a, an are supported on i?, we have 



\Gan{x) - Ga{x)\ < / \g{x,y)\\an{y) - a{y)\dy. 

IB 



Fix < e < 1. Since < cr, o"„, < M, we have by (3.6) and (3.7) 



B{x,e) 

Now let Bpi = B - B(x, e) and 



\g{x,y)\\an{y) - a{y)\dy <3Me\ 



(3.25) 



(3.26) 



An = {ye Bo : |cj„(y) - a{y)\ < e for all n > N}. 



Then by (3.24) 



\J An^Bo- DC{a) 

N>1 

where DC{a) is the set of discontinuities of a. As C{DC{a)) = and the sets A^ are 
monotone increasing in N ^ we can choose large enough so that 



C{Am) > C{Bo) 



\og{R/e) 



(3.27) 



Given n > N, for y € Bq we have 



R 

\g{x,y)\ < a^e^-^log- < 

e L{Bo - An) 



Splitting the right side of (3.25) into separate integrals over B{x,e), An, and Bq — An, we 
obtain 

\Ganix)-Ga{x)\<6Me^ + e j \g{x,y)\dy + 2Made 

< (3i?^logmax(i2, |x|) + (2arf + 3)M)e. 
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(ii) By Lemmas 3.1.8 and 3.1.15 the noncoincidence sets D(^n) ^-re contained in 



the ball Bi = M^/'^B. Given e > and a compact set K containing choose large 
enough so that |7„ — 7I < e on X for all n > N. Then 

s + e>7 + e>7n 



on K, so the function = max(s + e,jn) is superharmonic on K. By Lemma 3.1.9 
have fn>Sn, and hence s + e > Sn on K . Likewise 

Sn + e>7n + e>7 



we 



on K, so the function /„ = max(s„ + e,-y) is superharmonic on K. By Lemma 3.1.9 we 
have fn^s and hence s„ + e > s on X. Thus s„ — > s uniformly on K. 

(iii) Let /3 > be the minimum value of s — 7 on D,., and choose large enough 
so that I7 — 7„|, \s — Sn\ < /3/2 on D for all n > N. Then 

Sn-7n>s--f-(3>0 



on hence C -D(„) for all n > N. 



□ 



According to the following lemma, a nonnegative function on Z'^ with bounded 
Laplacian can grow at most quadratically. We will use this fact repeatedly. 

Lemma 3.1.18. Fix < /? < 1. There is a constant cp such that any nonnegative function 
f on Tj'^ with f{o) = and |A/| < X in B{o,R) satisfies 

/(x) < C/3A|x|^ xeB{o,PR). 

Proof. Consider the functions 



/±(y) = /(y)±A|yp. 



(3.28) 



By hypothesis, /+ is subharmonic and /„ is superharmonic on B(o, R). Given x G B{o, f3R), 
let r = min(i?, 2\x\), and let h± be the harmonic function on the ball B = B(o, r — 1) which 
agrees with f± on dB. Then 



f+<h+<h^+ 2Ar2 < /- + 2Xr'^. 
By the Harnack inequality |28l Theorem 1.7.2] there is a constant C/3 such that 

h+{x) < Cf3h+{o). 



(3.29) 
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Since /-(o) = /(o) = we have from (3.29) 



hence 



with C/3 = 8c/3 — 1. 



h+{o) < /_(o) + 2Ar2 < 8A|x|2. 



f{x) < h-\-{x) — A|x|^ < c/3A|x| 



□ 



The fohowing continuous version of Lemma 3.1.18 is proved in the same way, using 
the continuous Harnack inequality in place of the discrete one, and replacing |yp by \y\'^/2d 



in (3.28). 



Lemma 3.1.19. Fix < /5 < 1. There is a constant c'^ such that the following holds. Let 
f he a nonnegative function on M'^ with f{o) = 0, and let 

|2 



/±(x) = /(x)±A 



\x\ 
~2d 



If f+ is suhharmonic and /_ is superharmonic in B{o,R), then 

f{x) <Cf^X\x\'^, xeB{o,(3R). 
In the following lemma, denotes the inner e-neighborhood of A, as defined by 



(1.8) 



Lemma 3.1.20. Let A,BcM.'^ be bounded open sets. For any e > there exists rj > 
with 

{AUB)^ C A^UB^. 

Proof. Let K be the closure of {A U B)^, and let Yn = U -Bi/„,. Since K is contained in 
\jYn = AL) B, the sets K CiYn form an open cover of K, which has a finite subcover, i.e. 
K cYn for some n. □ 



3.1.5 Discrete Potential Theory 

Fix a sequence (5n i with 6i — 1. In this section we relate discrete superharmonic 
potentials in the lattice 6n'^'^ to their continuous couterparts in M*^. If j4 is a domain in 
6n'^'^, write 



A° = A + 



2 ' 2 
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points 


x'' = 






.° = .+ [-%,%]^ 


sets 


A--- = 


A n 5„Z'^ 




functions 


r = 









Table 3.1: Notation for transitioning between Euclidean space and the lattice. 



for the corresponding union of cubes in M*^. If ^ is a domain in M"^, write A'' = An 5n'^'^- 
Given x G M*^, write 

for the closest lattice point to x, breaking ties to the right. For a function / on Sn^'^, write 



X- = {x-—,x+ 2 



f°ix) = fix----) 

for the corresponding step function on W^. Likewise, for a function / on M'^, write /•' 
flSnZd-- These notations are summarized in Table 1. 

We define the discrete Laplacian of a function / on Sn^'^ to be 



A/(x) = 5-2(^^/(y)-/(x)) 



According to the following lemma, if / is sufficiently smooth, then its discrete Laplacian on 
dn'^'^ approximates its Laplacian on W^. 

Lemma 3.1.21. /// has continuous third derivative in a 6n-neighborhood of x E Sn'^'^, and 
A is the maximum pure third partial of f in this neighborhood, then 

\Af{x)-2dAf--{x)\<^Adn. 
Proof. By Taylor's theorem with remainder 

fix + SnCi) - 2f{x) + f{x-Snei) = -Q^^n + 6^(^ + *^')^^ ~ Q^^^ ~ ^^^^^^ 
for some < < 1. Summing over i = 1, . . . ,d and dividing by 5"^ gives the result. □ 
In three and higher dimensions, for x, y G 5^^'^ write 

gn{x, y) = Sl-'^giiS-^x, S'^y), d>3, (3.30) 
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where 

is the Green's function for simple random walk on T!^ . The scaling in (3.30) is chosen so 
that Axgn{x,y) = —6~'^l^rc=yy. In two dimensions, write 



9n{x, y) = -a{5j^ ^x, (5„ ^y) + - log 6n, d = 2, 

where 

a(x, y) = lim ( < m\Xk = x} - E^.#{fc < m\Xk = y}) 

m^oo 

is the recurrent potential kernel on Z^. 
Lemma 3.1.22. In all dimensions d>2, 

9n{x, y) = g{x, y) + 0{6l\x - y[~'^) 



(3.31) 



where g is given by (3.4)- 



Proof. In dimensions d > 3 we have from (3.30) and the standard estimate for the discrete 



Green's function (2.5) 



gn{x,y) = 5^ \ad 



\x - y\ 

5n 



2-d 



+ 



\x - y\ 

(5n 



ad\x-y\''-'' + 0{6l\x-y\~''). 



Likewise, in dimension two, using the standard estimate for the potential kernel (2.5) in 



(3.31 ) gives 



2 |x — y| 2 

gn{x, y) = - - log — h - log 5n + O 

n On TT 



\x - y\ 



. □ 



Our next result adapts Lemma 3.1.17[ i) to the discrete setting. We list here our 
standing assumptions on the starting densities. Let cr be a function on M'^ with compact 
support, such that 

< o- < M (3.32) 
for some absolute constant M. Suppose that a satisfies 



C{DC{a)) = 



(3.33) 
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where DC (a) denotes the set of points in R'^ where a is discontinuous. 
For n = 1, 2, . . . let be a function on dnl^'^ satisfying 



0< a„ < M 



and 



Finally, suppose that 

There is a ball B CM-^ containing the supports of a and for all n. 



(3.34) 
(3.35) 

(3.36) 



Although for maximum generality we have chosen to state our hypotheses on a and o"„ 
separately, we remark that the above hypotheses on ct„ are satisfied in the particular case 
when an is given by averaging a over a small box: 



f7„(x) = / a{y)dy 



In parallel to (3.3), for x E 6n'^'^ write 



(3.37) 



(3.38) 



where gn is given by (3.30) and (3.31). 
Lemma 3.1.23. If a, cr„ satisfy 13.3^^ - ^3.3^ , then 

(GnO-n)^ Ga 

uniformly on compact subsets o/M'^. 

Proof. Let K C M."^ he compact. By the triangle inequality 

|(G„a„)° - Gal < \{Gn(Tn)° - Ga°\ + \Ga° - Ga\. 

By Lemma 3.1.17[ i) the second term on the right side is < e/2 on for sufficiently large 
n. The first term is at most 



\{Gnan)°{x)-Ga°{x)\< ^ a„(y) 



^ian{x-\y) - I g{x,z)dz 



- ^ X] / \9nix--,y) - g{x,z)\dz 
\gn{x-\z--) - g{x,z)\dz. 



M 



(3.39) 
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By Lemma 3.1.22 we have 



\gnix-'- , z'-'-) - g{x, z)\ < \gn{x-'- , z'-'-) - g{x'-'- , z-'-)\ + \g{x-'- , z'-'-) - g{x, z) 

< C5n\x-Z\^-'^ 



for a constant C depending only on d. Integrating (3.39) in spherical shells about x, we 
obtain 

\{Gnanf - G<J^\ < CduJdMRSn, 

where R is the radius of B. Taking n large enough so that the right side is < e/2, the proof 
is complete. □ 



Our next result adapts Lemma 3.1.9 to the discrete setting. Let an be a function 
on with finite support, and let 7^ be the function on defined by 



^n{x) = -\x\'^ - Gn(Tn{x). 



Let 



Sn — inf{/(x)|/ is superharmonic on 5n TL^ and / > 7n} 
be the discrete least superharmonic majorant of 7n, and let 

Dn = {x^ 5n \Sn{x) > 7„(x)}. 



(3.40) 



(3.41) 



(3.42) 



Lemma 3.1.24. Letjn,Sn,Dn be given by m^-m^. Ifnc 6nZ'^ satisfies L>„ C n, 
then 

Sn{x) = inf{/(j;)|/ is superharmonic on 17 and / > 7„}. (3.43) 

Proof. Let / be any function which is superharmonic on 0, and > 7^. Since s„ is harmonic 
on Dn, so f — Sn is superharmonic on Dn and attains its minimum in Dn U dDn on the 
boundary. Hence / — s„ attains its minimum in 17 U dO, at a point x where Sn{x) = ^n{x). 
Since f > jn conclude that f > Sn on and hence everywhere. Thus s„ is at most the 



infimum in (3.43). Since the infimum in (3.43) is taken over a strictly larger set than that 
in (3.41), the reverse inequality is trivial. □ 
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3.2 Divisible Sandpile 



3.2.1 Convergence of Odometers 

By the odometer function for the divisible sandpile on 6n Z"^ with source density 
an, we will mean the function 

Un{x) = 6n ■ total mass emitted from x if each site y starts with mass (Jniv)- 

Theorem 3.2.1. Let Un he the odometer function for the divisible sandpile on b^^L'^ with 



source density an- If a, an satisfy (3.32)-(3.36}, then 



u. 



□ 



uniformly, 



where 



Ga{x), 



Ga is given by (3.3), and s is the least superharmonic majorant 0/7. 



Lemma 3.2.2. Let Dn be the set of fully occupied sites for the divisible sandpile in 6n'^'^ 
started from source density an- There is a ball J7 C M'^ with 



\jD°UDcn. 



n>l 



Proof. Let An be the set of fully occupied sites for the divisible sandpile in J^Z'^ started 



from source density t(x) = MIx^b, where B is given by (3.36). From the abelian property, 
Lemma [2.2.1[ we have D„ C An- By the inner bound of Theorem 1.2.2 if we start with mass 
m = 26^'^C{B) at the origin in {5^Z^, the resulting set of fully occupied sites contains B"'^ 
by the abelian property it follows that if we start with mass Mm at the origin in br^^ , the 



resulting set J7„ of fully occupied sites contains An- By the outer bound of Theorem 1.2.2 



is contained in a ball O of volume 3Mi2(i?). By Lemma 3.1.16 we can enlarge 17 if 
necessary to contain D. □ 

For X E (JnZ-^ write 



ln{x) 



\x^ — Gnan{x), 



where G„ is defined by (3.38). Denote by s„ the least superharmonic majorant of 7„ in the 
lattice 6nZ'^- 
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Lemma 3.2.3. Let Q be as in Lemma 3.2.2. There is a constant M' independent of n, 
such that |7„| < M' in 



Proof. By Lemma 3.1.22 we have for x S O'' 



\Gnan{,x)\ < (5^ ^ \gn{x,y)\an{y) 
<2M6i \9ix,y)\ 

yeB-- 

< CMR^ log R 

for a constant C depending only on d, where R is the radius of O. It follows that |7„| < 
{CM + I) R^ log RirvVi--. □ 

Lemma 3.2.4. Fix x G W^, and for y G ^nZ*^ let 

a{y) = 5ign[x-\y)- g{x,z)dz. (3.44) 

JyO 

There are constants Ci , C2 depending only on d, such that 

(i) |a(y)| <Ci5i+^|x-y|i-^. 

(ii) Ifyi ~ y2, then \a{yi) - a{y2)\ < C2Sl+'^\x - yi]'"^. 

if z G then 



Proof, (i) By Lemma 



3.1.22 



\9n{x-'-,y) -g{x,z)\ < \gn{x-\y) - g{x-\y)\ + \g{x-\y) - g{x,z) 
< C6n\x-y\'~'' 

for a constant C depending only on d. Integrating over y^ gives the result. 



(ii) Let p = y2 — yi- By Lemma 3.1.22 we have 



< 



y? 



y? 



\a{yi) -a(y2)| < I ^ \9n{x--,yi) -gn{x"\y2) - g{x, z) + g{x, z + p)\dz 
\9{x''',yi) - 9{x"",y2) - g{x, z) + g{x, z + p)\dz + 
+ 0{5i+^\x-yi\-''). 
Writing q = x'' — x — yi + z, the quantity inside the integral can be expressed as 

\f{o)-f[p)-f[q)+f[p + q)l 
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where 

f{w) = g{w,x'-- - yi). 
Since \p\, |g|, |p + ^| < 5n{'s/~d +1), we have by Taylor's theorem with remainder 



- f{p) - f{q) + f{p + q)\< M5l 

( 



dxidxj 



< < 



6£{d-l){d-2)adSl\x-yi\~'^, d > 3 
{24/TT)5l\x-yi\-^, d = 2. □ 



Lemma 3.2.5. Let J7 be an open ball as in Lemma 3.2.2. and let J7i be a ball with Q (Z^i. 
Let 

(j)n = -A{Snln:^). 

Then 

uniformly on Jl. 

Proof. Let I'nix) be the amount of mass present at x in the final state of the divisible 



sandpile on 6n'^ ■ By Lemma 2.2.2 we have = n„ + 7n, hence 



As„ = {l/n - (Tn) + {(Tn - 1) = fn - 1- 



In particular, |As„| < 1. 
Since 



and Sn = Gn4>n in we have for x E O 



= - ^ (t)n{y)a{y) 



where a{y) is given by (3.44). Hence 



G4>n{x) - Sn{x) = ^ Asn{y)a{y) + ^ {sn{y)a{z) - Sn{z)a{y)). (3.45) 



yen^ 



yr^z 
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Let R,Ri be the radii of By Lemma 3.2.4[ i), summing in spherical shells 



about X, the first sum in (3.45) is bounded in absolute value by 



< 



yen^ 



2Ri 



duJdr 



d-l 



)dr = CidbJdR\5n- 



To bound the second sum in (3.45), note that Sn = 7n outside Jl, so 



\sn{y)a{z) - Sn{z)a{y)\ < \jn{y)\\a{y) - a{z)\ + \a{y)\\jn{y) - 7n{z)\- (3.46) 
By Lemmas 3.2.3 and 3.2.4[ ii), the first term is bounded by 

hn{y)\Hy) - «(^)l < CsM'^^+'^lx - y\-^. 
Fix e > 0, and let 0,2 be a ball with 0,i C 0,2- Since 7 is uniformly continuous on 0,2, and 



7„ — > 7 uniformly on O2 by Lemma 3.1.23 for sufficiently large n we have 



hn{y) - 7n{z)\ < |7„(y) - 7(y)| + \-f{y) - 7(^)1 + h{z) - 7„(z)| < e. 
Thus by Lemma 3.2.4[ i) the second term in (3.46) is bounded by 

Hy)\\ln{y) - ln{z)\ < C,6i+''\x - yji-'^e. 
Since x £ 0, and y is adjacent to dOi, we have |x — ?/| > Ri — R — 5n, so the second term 



in (3.45) is bounded in absolute value by 

26~^#da-i (c2M'6l+'^{Ri - Ry^ + Ci6i+'^{Ri - R)^''^e) < C^e 



for sufficiently large n. 

Lemma 3.2.6. — > s uniformly on compact subsets ofM.'^. 



□ 



Proof. By Lemma 3.2.2 there is a ball 0, containing D and Z?^ for all n. Outside 0, we have 

□ □ 



uniformly on compact sets by Lemma 3.1.23 



To show convergence in Q, write 
s{x) = 



s{y)\-'v ( My, 



(3.47) 



68 



where rj is the standard smooth monifier 

Cei/(N'-i), |x| < 1 
0, |x| > 1 

normahzed so that Jj^a 7]dx = 1 (see 1171 Appendix C.4]). Then s is smooth and superhar- 



monic. Fix e > 0. By Lemma 3.1.2 'ii) and compactness, s is uniformly continuous on J7, so 



taking A sufficiently small in (3.47) we have |s — s| < e in $7. Let be the maximum third 



partial of s in $7. By Lemma 3.1.21 the function 



qn{x) = s'-'-{x) - l:Ae6n\x\'^ 
6 



is superharmonic in fi''. By Lemma 3.1.23 we have 7^ — > 7 uniformly in fl. Taking N large 



enough so that gAe(5„|xp < e in ^7 and |7„ — 7''| < e in i} ' for all n > N, we obtain 

Qn > s' — e > s" — 2e > 7 ' — 2e > 7„ — 3e 

in rj''. In particular, the function /„ — max(^j2 ~t~ 36, 7n) is superharmonic in 17". By 
Lemma 3.1.24| it follows that fn > Sn, hence 



Sn < + 3e < s" + 3e < s" + 4e 

in 0,''. By the uniform continuity of s on 0,, taking larger if necessary we have |s — < e 
in Q, and hence < s + 5e in $7 for all n > N. 
For the reverse inequality, let 

= -A(s„ln.:) 



where J7i is an open ball containing $7. By Lemma 3.2.5 we have 



< e 



and hence 



Gcb°>l°-e>^- 2e. 



on 17 for sufficiently large n. Since </>^ is nonnegative on 17, by Lemma 
G(j)^ is superharmonic on 17, so by (3.48) the function 



3.1.4 



(3.48) 
the function 



V'n = max(G(/)° + 2e, 7) 
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is superharmonic on 0, for sufficiently large n. By Lemma 3.1.9 it follows that ipn ^ s, hence 



> G(j)^ — e > s — 3e 



on n for sufficiently large n. 



□ 



Proof of Theorem \3.2.1\ Let $7 be as in Lemma [3. 2. 2[ By Lemmas |3 . 1 . 2 3| and |3 . 2 . 6] we have 

In 



^ 7 and ^ s uniformly on 0,. By Lemma 



2.2.2 



we have Un = Sn — Jn- Since 



tt^ = = s — 7 off O, we conclude that — > s — 7 uniformly. 



□ 



3.2.2 Convergence of Domains 



□ 



In addition to the conditions (3.32 )-(3.36 ) assumed in the previous section, we 



assume in this section that the source density a satisfies 

For all X G M*^ either cr(x) < A or (j{x) > 1 
for a constant A < 1. We also assume that 



(3.49) 



(3.50) 



{a > 1} = {o- > l}''. 
Moreover, we assume that for any e > there exists A^(e) such that 

If X G {cj > then fT„(x) > 1 for all n > iV(e); 

and 

If X ^ {cj > ly, then cj„(x) < A for ah n > N{e). 

As before, we have chosen to state the hypotheses on a and Cn separately for maximum 
generality, but all hypotheses on (T„ are satisfied in the particular case when cr„ is given by 



(3.51) 
(3.52) 



averaging o" in a small box (3.37). 
We set 

7(x) = — |xp — Ga{x) 
with s the least superharmonic majorant of 7 and 

D = {xe R'^\s{x) > 7(x)}. 

We also write 

5 = L>U{xGM^|cj(x) > 1}°. 

For a domain A C M'^, denote by and A"" its inner and outer open e-neighborhoods, 
respectively. 
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Theorem 3.2.7. Let a and Un satisfy {3.3^ - \3.3^ and \340i) - \3.5^ . For n > 1 let Dn 

be the domain of fully occupied sites for the divisible sandpile in b^^'^ started from source 
density On- For any e > we have for large enough n 



(3.53) 



According to the following lemma, near any occupied site x € Dn lying outside W", 
we can find a site y where the odometer Un is relatively large. This is a discrete version of a 
standard argument for the obstacle problem; see for example Friedman [19], Ch. 2 Lemma 
3.1. 

Lemma 3.2.8. Fix e > and x G Dn with x ^ D''. If n is sufficiently large, there is a 
point y G Sn"^"^ with \x — y\ < | + 5n and 



Un{y) > Unix) + 



1 - A 



Proof By (3.50), we have {cr > 1} c D. Thus if x ^ D^, the ball B = B{x, §)•• is disjoint 
from {cr > l}*^/^. In particular, if n > A^(§), then by (3.52) we have (T„ < A on B. Thus 
the function 

w{y) = Un{y) - (1 - \)\x - 

is subharmonic on BODn, so it attains its maximum on the boundary. Since w{x) > 0, the 
maximum cannot be attained on dDn, where Un vanishes; so it is attained at some point 
y € dB, and 

n„(y) >u;(y) + (l-A)(|)'. 
Since w{y) > w{x) = Un{x), the proof is complete. □ 



Proof of Theorem 3.2.7. Fix e > 0. By Lemma 3.1.20 we have 

D^ C Dr,U {a > l}r, 

for some r/ > 0. Let u = s — -f. Since the closure of Z?^ is compact and contained in D, we 
have u > on for some > 0. By Theorem 3.2.1 for sufficiently large n we have 
Un > u — ^rriri > on D^j, hence C Dn- Likewise, by (3.51) we have {a > 1}':^ C Dn for 
large enough n. Thus D^ C Dn- 

For the other inclusion, fix x € (5„Z'^ with x ^ D"". Since u vanishes on the ball 
B = B{x, |), by Theorem 3.2.1 we have n„ < ^-^e^ on B ' U dB " for all sufficiently large 
n. By Lemma 3.2.8 we conclude that x ^ Dn, and hence Dn C D'' ''. □ 
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3.3 Rotor-Router Model 



In trying to adapt the proofs of Theorems 3.2.1 and 3.2.7 to the rotor-router 
model, we are faced with two main problems. The first is to define an appropriate notion 
of convergence of integer- valued densities (Jn on to a real- valued density cr on M . The 
requirement that (T„ take only integer values is of course imposed on us by the nature of 
the rotor-router model itself, since unlike the divisible sandpile, the rotor-router works with 
discrete, indivisible particles. The second problem is to find an appropriate analogue of 
Lemma |2. 2. 21 for the rotor-router model. 

Although these two problems may seem unrelated, the single technique of smooth- 
ing neatly takes care of them both. To illustrate the basic idea, suppose we are given a 
domain j4 C M"^, and let (t„ be the function on Sfi^j taking the value 1 on odd lattice 
points in A\ and the value 2 on even lattice points in A'\ while vanishing outside A'. We 
would like to find a sense in which cr„ converges to the real- valued density a = |1a- One 
approach is to average Un in a box whose side length L„ goes to zero more slowly than the 
lattice spacing: j while Ln/^n t oo as n | co. The resulting "smoothed" version of an 
converges to a pointwise away from the boundary of A. 

By smoothing the odometer function in the same way, we can obtain an approx- 



imate analogue of Lemma 2.2.2 for the rotor-router model. Rather than average in a box 
as described above, however, it is technically more convenient to average according to the 
distribution of a lazy random walk run for a fixed number a(n) of steps. Denote by (Xjfc),fc>o 
the lazy random walk in (J^Z'^ which stays in place with probability ^ and moves to each of 
the 2d neighbors with probability Given a function / on 5n'^'^, define its k-smoothing 

=E(Xfc|Xo = x). (3.54) 

From the Markov property we have SkS^ = Sk+i- Also, the discrete Laplacian can be 
written as 

A = 25-2(^1 -5o). (3.55) 

In particular, lS.Sk = SkA. 

3.3.1 Convergence of Odometers 

For n = 1,2,... let an be an integer- valued function on 6n'^'^ satisfying < cr„ < 
M. We assume as usual that there is a ball B CM.'^ containing the support of an for all n. 
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Let (T be a function on M'^ supported in B satisfying (3.32) and (3.33). In place of condition 



(3.35) we assume that there exist integers a{n) j oo with (5„a(n) [ such that 

('S'a(„)0-„)°(x) a{x), X ^ DC{a). 



(3.56) 



By the odometer function for rotor-router aggregation starting from source density 
an, we win mean the function 

Un{x) = 5^ ■ number of particles emitted from x 

if (Jn{y) particles start at each site y. 

Theorem 3.3.1. Let Un he the odometer function for rotor-router aggregation on 5n'^'^ 



starting from source density an- If a, an satisfy (3.32)-(3.34), (3.36) and (3.56), then 
li^ — > s — 7 uniformly, where 



7(x) 



\x\ 



g{x,y)a{y)dy 



f{v) - fix) 



and s is the least superharmonic majorant of ^. 

Given a function / on Sfi'^ , for an cdgG (x, write 

V/(x,y) = 

Given a function k on edg es in Syi*^^ f write 

divK;(x) = ^^^i^ix,y)- 

y 

The discrete Laplacian on 5„,Z'^ is then given by 
Af{x) = divVf = 6~^ 



The following "rescaled" version of Lemma 2.4.1 is proved in the same way. 

Lemma 3.3.2. For an edge {x,y) in b^'^ , denote by K{x,y) the net number of crossings 
from X to y performed by particles during a sequence of rotor-router moves. Let 

u{x) = 6n ■ number of particles emitted from x during this sequence. 

Then 

Vu{x, y) = 5n{-2dK{x, y) + p{x, y)). (3.57) 
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for some edge function p which satisfies 

\p{x,y)\ <M-2 
independent of x, y, n and the chosen sequence of moves. 

Proof. Writing N{x^ y) for the number of particles routed from x to y, for any y,z ^ x we 
have 

\N{x,y)-N{x,z)\ < 1. 
Since u{x) = 6'^ J2yr^x ^{^^ y)i obtain 

6~^u{x) -2d + l< 2dN{x, y) < 5~'^u{x) + 2d-l 

hence 

\Vuix,y)+2d6nK{x,y)\ = 6n\S~'^u{y) - 6~'^u{x) + 2dN{x,y) - 2dN{y,x)\ 

< (4d - 2)6n. □ 

Recall that the divisible sandpile odometer function has Laplacian 1 — (T„ inside 
the set Dn of fully occupied sites. The next lemma shows that the same is approximately 
true of the smoothed rotor-router odometer function. Denote by Rn the occupied shape for 
rotor-router aggregation on (5„Z°' starting from source density C7„. 

Lemma 3.3.3. |A5fcU„(x) - F^{Xk £ Rn) + 5fcCT„(x)| < Co/VkTl. 



Proof. Let k and p be defined as in Lemma 3.3.2 Since each site x starts with an{x) 
particles and ends with either one particle or none accordingly as x S Rn or x ^ Rn, we 
have 

2d6ndiv n = an - Ir„. 



Taking the divergence in (3.57) we obtain 

Aun = 1r„ - (Jn + 5n div p. (3.58) 
Using the fact that Sk and A commute, it follows that 

ASfcn„(x) = P^(Xfc G Rn) - E^.c7„(Xfc) + 5n div p{Xk). (3.59) 
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Since Vn and k are antisymmetric, p is antisymmetric by (3.57). Thus the last term in 



(3.59) can be written 

5n div p{Xk) = ^Y1 ^-^^^ = y^^^y^ ' 

^Y,^,{Xk = z)p{y, 



yr^z 



\- ^(Px.(Xfc = y)- P,(Xfc = z))p{y, , 



4d 



(3.60) 



yr^z 



where the sums are taken over ah pairs of neighboring sites y,z ^n^'^ ■ 

We can couple lazy random walk started at x with lazy random walk X'^ 
started at a uniform neighbor of x so that the probability of not coupling in k steps is at 
most C I^Jk + 1, where C is a constant depending only on d [36 . Since the total variation 
distance between X]. and X'^ is at most the probability of not coupling, we obtain 

C 



^|P, {Xk = y)-^Mk = z)\ < 



y^z 



Using the fact that p{y,z) is uniformly bounded in (3.60), taking Cq = ACd completes the 
proof. □ 

The next lemma shows that smoothing the odometer function does not introduce 
much extra error. 

Lemma 3.3.4. \SkUn - Un\ < 5l{\Mk + CQ^/k + l). 



Proof. From (3.55) we have 



k-l 



\SkUn Un\ ^ ^ ^ I Sj-\-lUn SjUr, 
j=0 



6\ 



j=0 



But by Lemma 3.3.3 



\ASiUn\ <M + 



Co 
VJTT' 



Summing over j yields the result. 



□ 
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Let 



ln{x) = - G„5„(„)Crn(x) 



and let Sn be the least superharmonic majorant of 7n. By Lemma 2.2.2 the difference Sn—^n 
is the odometer function for the divisible sandpile on starting from the smoothed 

source density a„ = Sa[n)^n- Note that by Lemma 3.2.2 there is a ball 17 C M*^ containing 
the supports of s — 7 and of Sn — 'Jn for all n. The next lemma compares the smoothed 



rotor-router odometer for the source density cr„ with the divisible sandpile odometer for the 
smoothed density 

Lemma 3.3.5. Let 17 C M'^ be a ball centered at the origin containing the supports of s — ^ 
and of Sn — 7n, for all n. Then 



on all of bn%^ , where r is the radius ofQ,. 
Proof. Since 



7n - Cor^ainyy^ 



(3.61) 



A7n = -1 + Sa(n)(^n, 



by Lemma 3.3.3 the function 

f{x) = Sa(n)Un{x) + 7n(x) + Coa(n)~^/^(r^ - 
is superharmonic on br^'^ ■ Since / > 7n on Jl'', the function il)n = max(/, 7^) is superhar- 



monic on Jl'', hence V'n ^ Sn by Lemma 3.1.24 Thus f > Sn on Q'', so (3.61) holds on Jl" 
and hence everywhere. □ 

Lemma 3.3.6. Let Rn be the set of occupied sites for rotor-router aggregation in b^^'^ 
started from source density an- There is a ball ncR'^ with \JR° cn. 

Proof. By assumption there is a ball B C M.'^ containing the support of (T„ for all n. Let 
An be the set of occupied sites for rotor-router aggregation in bn'^'^ started from source 
density t{x) = Ml^g^- From the abelian property we have i?„ C An- By the inner 

if we start with [25~'^£(i?)J particles at the origin in bn'^'^, the 



bound of Theorem 



1.2.1 



resulting set of occupied sites contains B ' ; by the abelian property it follows that if we 
start with M [2b-'^C{B)\ particles at the origin, the resulting set 0„ of fully occupied sites 



contains An- By the outer bound of Theorem 1.2.1 r2„ is contained in a ball 17 of volume 
3MC{B). □ 
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Proof of Theorem 3.3.1 By Lemma 3.2.2 there is a ball 17 containing the support of s — 7 



and of Sn — 7n for all n. By Lemma 3.3.6 we can enlarge if necessary to contain the 



support of 5c,(„)Mn for all n. By Lemma 3.3.3, the function 

(l){x) = Sa(n)Un{x) - Sn{x) + 7n(x) + Coa(n)~^/^|x|^. 

is subharmonic on the set 

Rn = {x G Rn : y £ Rn whenever ||x — y||i < a(n)} 



since (-'^^a(n) ^ Rn) = 1 for X G From Lemma 3.3.3 we have |AS'„(„)n„| < M + Co, 
so for X ^ Rn we have by Lemma 



3.1.18 



Sa{n)Un{x) < c{M + Co)a(n)^5^. 
By the maximum principle in Rn 

(/.(x) < c{M + Co)a{nf5l + Coa(n)-i/V2, x G Rn- 
From Lemma 13.3.51 we obtain 

- Cor2a(n)-i/2 < S^(n)Un -Sn + ln< c{M + CQ)a{nf5l + Cor'a(n)-i/2 (3.62) 
on all of 5n'^'^, where r is the radius of 



By Lemmas 



3.1.23 



and 



3.2.6 



we have 7 and — > s uniformly on $7. Since 



a(n) t 00 and (5„Q(n) | 0, we conclude from (3.62) that (5'Q,(„)ii„)^ — > s — 7 uniformly on 
rj. Since both 5'q,(„)U„ and s — 7 vanish outside this convergence is uniform on M*^. By 
Lemma 13.3.41 we have 

|'S'a(}i)'Wn ~ "Un] ^ 

uniformly on 6n'^'^, and hence v!^ — > s — 7 uniformly on M*^. □ 

3.3.2 Convergence of Domains 

In addition to the assumptions of the previous section, in this section we require 



that 



We also assume that 



For ah X G M'^ either cr(x) > 1 or a{x) = 0. 



{a > 1} = {o- > 1}^ 



(3.63) 
(3.64) 
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Moreover, we assume that for any e > there exists A^(e) such that 



and 



If X G {cr > then fT„(x) > 1 for all n > N{e); 



If X ^ {cj > ly, then o-„(x) = for ah n > N{e). 



(3.65) 



(3.66) 



Theorem 3.3.7. Let a and satisfy (3.32)-(3.34), (3.36), (3.56) and (3.63)-(3.66). For 



n > 1 let Rn = {un > 0} be the domain of sites in 611^ that emit a particle in the votov- 
router model started from source density an. For any e> Q we have for all sufficiently large 
n 

D'i cRnC D''\ (3.67) 

where 

5 = {s > 7} U{o- > 1}". 
Lemma 3.3.8. Fix e > and n > N{e/2). Given x ^ D'^ and 5n < p < e/2, let 



Mpix) = #B{x,p)nR„ 



If Un < S^m on B{x,p), then 



Tfl 



Proof. Since at least Mp{x) particles must enter the ball B{x,p), we have 

yedB{x,p-5n) 



There are Mp[x) — Mp-s„{x) terms in the sum on the left side, and each term is at most 

m{J\fp{x) - Mp-s^x)) > Mp{x). □ 



5^m, hence 



The following lemma can be seen as a weak analogue of Lemma 3.2.8 for the 



divisible sandpile. In Lemma 3.3.11 below, we obtain a more exact analogue under slightly 
stronger hypotheses. 
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Lemma 3.3.9. For any e > 0, if n is sufficiently large and x G i?„ with x D"", then there 
is a point y £ Rn with \x — y\ < e/2 and 

, N y ^ 

''"^^^ - 4 log {2Me/2dnr)' 
Proof. Take n > N{e/2) large enough so that 4(5n < e. Let m be the maximum value of 



on B{x,e). Since x G i?„ we have Mo{x) = 1. Iteratively applying Lemma 3.3.8 

obtain 



we 



m 



m — 1 



J\fo{x) > exp 



Since 



we conclude that 



M,/2{x) <#B{x,e/2) <2iOd 



25n) 



< 



log {2ua{e/5n)''). 



Solving for m yields the result. 



□ 



The following lemma shows that far away from D, the rotor-router odometer grows 
at most quadratically as we move away from the boundary of 

Lemma 3.3.10. For any e > and k >1, if n is sufficiently large and z U D^, then 

Un{x) < 4cdk^S'^, X G B{z, kSn)" ■, 



where c = 0^/2 ^-^ the constant in Lemma 3.1.18 



Proof. Taking the divergence in (3.57), since — 2d5„, div ^(a;) is the net number of particles 
entering x, we obtain 

Au„ = 1r„ - cFn + 5n div R. (3.68) 



By (3.66), if n > N{e/2) is large enough so that k6n < e/4, we have cr„ = on the ball 



B = B{z, 2kdn)"- Since \R\ <Ad — 2, we have |Aun| < 4(i on B, so by Lemma 3.1.18 

Un{x) < Acd\x — z\^ 

for X G 5(2;, k5n)"- 



□ 



The following lemma is analogous to Lemma 3.2.8 for the divisible sandpile. 
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Lemma 3.3.11. Fix e > and k > 4Co, where Cq is the constant in Lemma 3.3.3 There 



is a constant Ci such that if x £ Rn, x ^ satisfies 

SkUn{x) > Cik'^51 

and n is sufficiently large, there exists y G ^^Z"^ with \x — y\ < f + (5n and 



Proof. Let 



SkUniy) > SkUn{x) + 



R^^ = {x S Rn : y £ Rn whenever ||x — y||i < k6n}. 



By Lemma 3.3.3, we have 

ASkUn{y)>l-Skan{y)-Co/Vk, y & R^n'^ ■ 

Note that CQ/\/~k < ^. Take n > N^e/A) large enough so that k5n < e/4; then SkCTn 
vanishes on the bah B = B{x, e/2)'' by ( 3.66| ). Thus the function 



f{y) = SkUniy) -hx -y 



(k) 

is subharmonic on i?„ := Rn Ci B, so it attains its maximum on the boundary. 



If y G dR'n\ there is a point z i?„ with |y — z| < k6n. By Lemma 



3.3.3 



we have 



|Au„| < M + Co, hence by Lemma 3.1.18 it follows that 



Un < 4c{M + Co)k'^6l 

in the ball B{z, 2k6n). Taking Ci = 4c(M + Cq) we find that 

SkUniy) < Cik'^6l < SkUn{x). 

Thus / cannot attain its maximum in Rn on dRn , and hence must attain its maximum at 
a point y G dB. Since f{y) > f{x) we conclude that 

1, 



SkUniy) = f{y) + ok - y 
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Proof of Theorem 3.3.7 Fix e > 0. By Lemma 3.1.20 we have 



D, C Dr,\J {a > l}r, 

for some rj > 0. Since the closure of D^j is compact and contained in D, we have u > on 
for some > 0. By Theorem 3.3.1 for sufficiently large n we have Un> u — \rnr^ > 



on D':I, hence D':: C Rn- Likewise, by (3.51) we have {a > 1}^ C Rn for large enough n. 



Thus D'^ C Rn. 

For the other inclusion, fix x G (^nZ*^ with x ^ W". Since u vanishes on the 
ball B = B(x,e), by Theorem 



3.3.1 



we have u„ < ^ on B' for all sufficiently large n. 
Let k = 4Co, and take n large enough so that k6n < e/4. Then SkUn < ^ on the 
ball i?(x, 3e/4)''. By Lemma 3.3.11 it follows that S^Un < Cik^S"^ on the smaller ball 



B' = B{x,e/^Y'' . In particular, for y £ B' we have 



SkUniu) 



< 



C2k''/^+^6l 



MXk = y) 

For sufficiently large n the right side is at most e(5„/41og (2iOci{^/26n)'^), so we conclude from 



Lemma 3.3.9 that x 



□ 



3.4 Internal DLA 

Our hypotheses for convergence of internal DLA domains are the same as those for 
the rotor-router model: cr is a bounded, nonnegative, compactly supported function on M'^ 
that is continuous almost everywhere, and {an}n>i is a sequence of uniformly bounded func- 
tions on Sfi'^ with uniformly bounded supports, whose "smoothings" converge pointwise 
to a at all continuity points of a. 



Theorem 3.4.1. Let a and satisfy (3.32)-(3.34), (3.36), (3.56) and (3.63)-(3.66). For 



n> 1 let In be the random domain of occupied sites for internal DLA in 5„Z'^ started from 
source density an. For all e > we have with probability one 



D'^ C In C D'^" for all sufficiently large n. 



(3.69) 



where 



D = {s>j}\j{a> 1}°. 
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3.4.1 Inner Estimate 

Fix n > 1, and label the particles in b^r{L'^ by the integers 1, . . . , m„, where = 
XlxG^nZ'* cr„(x). Let Xi be the starting location of the particle labeled so that 

= x} = cr„(x). 

For each i = 1, . . . ,m„ let (X^)t>o be a simple random walk in such that Xq = x^, 
with X* and independent for i ^ j. 

For z E Jn^*^ and e > 0, consider the stopping times 

Ti=ini{t>{)\Xl = zY 
Ti=ini{t>^\Xi^D--]; 

v^ = \ni{t>f)\Xl<^{Xl^yr}^]. (3.70) 

The stopping time z^* is defined inductively in i with z^^ = 0. We think of building up the 
internal DLA cluster one site at a time, by letting the particle labeled i walk until it exits 
the set of sites already occupied by particles with smaller labels. Thus is the number of 
steps taken by the particle labeled i, and is the location where it stops. 
Fix z G and consider the random variables 

m„ 

1=1 

m„ 

Le = l{;y!<T|<T|}- 
i=l 

These sums can be interpreted in terms of the following two-stage procedure. During the 
first stage, we form the occupied internal DLA cluster by adding on one site at a time, as 
described above. In the second stage, we let each particle continue walking from where it 
stopped until it exits D];. Then Ai^ counts the total number of particles that visit z during 
both stages, while counts the number of particles that visit z during the second stage. 
In particular, if < Me, then z was visited during the first stage and hence belongs to 
the occupied cluster. 

The sum Le is difficult to estimate directly because the indicator random variables 
in the sum are not independent. Following [29], we can bound Le by a sum of independent 
indicators as follows. For each site y G D'^, let {Y^^)t>o be a simple random walk in (J^Z*^ 
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such that Yq = y, with and Y^ independent for x ^ y. Let 



where 



r|=inf{t>0|y,^ = z}; 
r*' = inf {t >^\Yy ^D'-]. 

Thus counts the number of walks Y^ that hit z before exiting D'^. Since the sites X* 
are distinct, we can couple the walks {^^} and {X*} so that < L^. 
Define 

where 

ff„,,(y,z)=E#{t<r,^|y/ = z} 
is the Green's function for simple random walk in ^nZ*^ stopped on exiting D'-^. Then 



1=1 



QnAz.z) Y^{an{y)-\)^{Tl<Tt) 



(3.71) 



where in the last step we have used the identity 



Thus fn,e solves the Dirichlet problem 

Afn,e = d-Hl-an), 
fn,e = 0, 



on D: 



on (9L); 



(3.72) 



(3.73) 



Note that the divisible sandpile odometer function u„ for the source density a„ 



solves exactly the same Dirichlet problem, with the domain D'^ in (3.73) replaced by the 



domain Dn of fully occupied sites. Our strategy will be first to use Theorem 3.2.7 to argue 
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that since — > D, the solutions to the Dirichlet problems in Dn and D]; should be close; 
next, by Theorem 3.2.1 since u„ ^ u it follows that the functions fn,e and u are close. Since 
u is strictly positive in D^, we obtain in this way a lower bound on fn^e, and hence a lower 
bound on — L^). Finally, using large deviations for sums of independent indicators, 

we conclude that with high probability < Me, and hence with high probability every 



point z S belongs to the occupied cluster /„. The core of the argument is Lemma 3.4.3 
below, which gives the desired lower bound on fne- 



In order to apply Theorems 3.2.1 and 3.2.7 we must have discrete densities cr„ 
which converge pointwise to a at all continuity points of a. Recall, however, that in order 
to allow a to assume non-integer values, we have chosen not to assume that (T„ converges 
to fj, but rather only that the smoothed density Sa(n)(^n converges to cr; see (3.56). For 
this reason, we will need to run the divisible sandpile on (JnZ'^ starting from the smoothed 
density rather than from cj,„, and we will use the following smoothed version of equation 



(3.73). 



on D '). 



(3.74) 



Here /c > 1 is arbitrary, and e' = e + k5n- 



The following lemma is proved in the same way as Lemma 3.3.4 



Lemma 3.4.2. Let e' = e + k5n- Then for z E D'^, we have 

\SkfnA^) - fnA^)\ < \Mk. 



Proof. From (|3.55|) and (|3.74|) we have for z G D'^, 

\Sj+lfn,e{z) - Sjfn,e{z)\ 



\Skfn,e{^) - fnA^)\ ^ 



El 

j=0 



f^\AS,fnAz)\ 



3=0 



< -Mk. □ 
- 2 



Lemma 3.4.3. Fix e > 0, and let P > be the minimum value of u = s — 'j on D^. There 
exists < r] < e such that for all sufficiently large n 

fn,viz) > \l35-\ z G D--. 
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Proof. Since u is uniformly continuous on D, we can choose r] > small enough so that 
u < (3/8 outside Z)2r?- Let ry' = r/ + (5„a(n). Since 5na{n) | 0, we have u < (3/8 outside Z)^/ if 
n is sufficiently large. Since u > (3 on D^, for large enough n we have dD^i <f_ and hence 
r( < e. Let Un be the odometer function for the divisible sandpile on Sji'^ started from 
source density Sa(n)(^n- We have — > n uniformly by Theorem 



3.2.1 



so for n sufficiently 



large we have \un — m''| < (3/8, hence u„ < /3/4 on dD';^,. 



Let Dn be the domain of fully occupied sites for the divisible sandpile in (5„Z° 



started from source density Sa(n)^n- By Theorem 3.2.7 we have D]^ C Dn for all sufficiently 
large n. Now (3.73) implies that Sa(^n)fn,ri — S~'^Un is harmonic on Z)^,. By Lemmas 
and 



3.4.2 



3.1.18 



for z E dD';^, we have 



\Sa{n)fn,-n{z)\ < \fn,ri{z)\ + ^Ma{n) 

< cMa{nf + ^Ma(n). 

Since 6na{n) [ 0, for sufficiently large n we have 

(3 

Sa(n)fn,r,{z) > "^^^ ^ ^ ^^V' 

Since Sa[n)fn,ri — ^n'^i^n attains its minimum in D'-^, on the boundary, for z G D]^, we have 

~ if) 

7^ 
16 



'S'a(n)/n,ry(^;) > [Un{z) 



> S' 



u{z) 



Since u > /? in D];, taking n large enough so that a{n)5l < (3/8M, we conclude from 
Lemma 13.4.21 that 

9/3 ^ 



1652 



on Dy. 

Lemma 3.4.4. We have 
Moreover 



\Ma{n) > \(35-\ 



□ 



9n,ri{z, z) 



ETW^ < M- 

9n,ri\Z, Z) 
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Proof. Let 



r| = inf{t > oiy^^ = z} 



be the first hitting time of z for the walk Y^. By (3.72 ) and the symmetry of gn,r), we have 



EL, = ^ F{Ty < ry) 

,r,{V,z) 



(3.75) 



yeD' 



. 9n,ri{z, z) 



= 7 ^ 9n,rj{z,y). 

9n,rj{Z,Z) ^^^^ 

The smii in the last line is KzTqd-^^. 

To prove the inequality for note that KAi^] is bounded above by M times 

the sum on the right side of (3.75 ). □ 

The next lemma, using a martingale argument to compute the expected time for 
simple random walk to exit a ball, is well known. 

Lemma 3.4.5. Fix r > and let B = B{o,rY'- C (5„Z'^, and let T he the first hitting time 
ofdB. Then 

r 



+ 



Proof. Since — t is a martingale with bounded increments, and EqT < oo, by 

optional stopping we have 



EoT = 6-^Eo\Xt\^ = 5-\r + 0{6n)f. □ 



The next lemma, which bounds the expected number of times simple random walk 
returns to the origin before reaching distance r, is also well known. 

Lemma 3.4.6. Fix r > and let B = B{o,ry'' C Sn'^'^, and let Gb be the Green's function 
for simple random walk stopped on exiting B. If n is sufficiently large, then 

r 

Gb{o,o) < log— . 

On 
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Proof. In dimension d > 3 the result is trivial since simple random walk on 5„Z'^ is transient. 
In dimension two, consider the function 

f{x) = Gsio, x) - gn{o, x) 



where gn is the rescaled potential kernel defined in (3.31). Since / is harmonic in B it 



attains its maximum on the boundary, hence by Lemma 3.1.22 we have for x ^ B 



f{x)<-logr + o(^]. 



□ 



Since gn{o, o) = f log^m the result follows on setting x = a. 

We will use the following large deviation bound; for a proof, see [2] Cor. A. 14]. 



Lemma 3.4.7. If N is a sum of finitely many independent indicator random variables, 
then for all X > 

F{\N - EN\ > XEN) < 26-"^^^ 
where cx > is a constant depending only on X. 



Let 



in = {Xi^W' < f'} Cin 



where is given by (3.70), and 



= inf {t > 0\Xi ^ D--}. 



The inner estimate of Theorem 3.4.1 follows immediately from the lemma below. Although 



for the inner estimate it suffices to prove Lemma 3.4.8 with /„ in place of In, we will make 



use of the stronger statement with In in the proof of the outer estimate in the next section. 
Lemma 3.4.8. For any e > 0, 

P(D(;' C In for all but finitely many n) = 1. 

Proof. For z E let £z{n) be the event that z ^ In- By Borel-Cantelli it suffices to show 
that 
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By Lemma 



3.1.20 



since D = D U {a > 1}° we have 

C D,, U{a> 1}, 



for some e' > 0. By (3.65), for n > iV(e') the terms in (3.76) with z G {cr > Ij^' vanish, so 
it suffices to show 

n>l zeD'-'; 



By Lemma 3.4.3 there exists < r] < e' such that 



(3.78) 



for all sufficiently large n, where /3 > is the minimum value of u on D^/. Fixing z G D;'/, 



since Lrj < Lrj we have 



in)) < = L^) 

< ^{M^ < L^) 

< ^{Mri <a)+ F{Lri > a) 



(3.79) 



for a real number a to be chosen below. By Lemma 
independent indicators, we have 



3.4.7 



since and My are sums of 



(L^ > (1 + A) E Ly) < 26-''^^^'^ 
{My < (1 - A) E My) < 2e-"^^-^" 



(3.80) 



where cx depends only on A, chosen below. Now since z G D^/ and D is bounded by 
Lemma 3.1.16[ we have 

B{z,e -7]) CDyC B{z,R) 



hence by Lemma 3.4.5 



(3.81) 



By Lemma 3.4.4 it follows that EMy < MR^ /6lgn,y{z, z). Taking 



a = EL„ + 



4^Sl9n,r,{z,z) 



in (3.79), and letting A = P/4MR^ in (3.80), we have 



(3.82) 
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hence 



a > (1 + A)EL^. 



Moreover, from (3.78) we have 



hence by (3.82) 



a < E,M 



< (1 - A) EA^^^. 



Thus we obtain from (3.79) and (3.80) 



(3.83) 



By Lemmas 3.4.4 and 3.4.6 along with (3.81) 



EL. 



> 



9n,ri{z,z) 

1 fe' -T] 



2\ 5n ) \og{R/5n) ■ 



Using (3.83), the sum in (3.77) is thus bounded by 



^ ^ P(f.(n)) < S-'u^.R' ■ 4exp (- 

n>lz&D''; n>l ^ 



cx{e'-rjf 
2Sllog{R/6r, 



< oo. □ 



3.4.2 Outer Estimate 

For X € write 

Q{x,h) = {yeZ'^ : ||x-y||oo < h} 

for the cube of side length 2/i + 1 centered at x. According to the next lemma, if we start 
a simple random walk at distance h from a hyperplane H C Z*^, then the walk is farily 
"spread out" by the time it hits H, in the sense that its chance of first hitting H at any 
particular point z has the same order of magnitude as z ranges over a (d — l)-dimensional 
cube of side length order h. 
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Lemma 3.4.9. Fix y with yi = h, and let T he the first hitting time of the hyperplane 
H = {x ^ U^\x\ = 0}. Let F = H n Q{y, h). For any z ^ F we have 

¥y{XT = z)> ah^-'^ 
for a constant a depending only on d. 
Proof. For fixed w £ H, the function 

fix) = F,{Xt = w) 

is harmonic in the ball B = B{y, h). By the Harnack inequality |28| Theorem 1.7.2] we have 

/(x)>c/(y), x£B{y,h/2) 

for a constant c depending only on d. By translation invariance, it follows that for w' G 
H n -B(w, /i/2), letting x = y + w — w' we have 

¥y{XT = w') = fix) > cfiy) = cFy {Xt = w). 

Iterating, we obtain for any w' £ H f] Q{w, h) 

Fy {Xt = w') > Fy {Xt = w). 



(3.84) 



Let T' be the first exit time of the cube Q{y, h — 1). Since F is a boundary face of this cube, 
we have {Xt' G F} C {Xt € F}. Let zq be the closest point to y in H. Taking w' = zq, 



we have w' £ H D Q{w, h) whenever w £ F. Summing (3.84) over w £ F, we obtain 

1 

2d' 



{2h + l)^-ic~^ Fy {Xt = zo) > Fy{XT e F) > F{Xt' G F) 



Now for any z £ F, taking w = zq and w' = z in (3.84), we conclude that 
Fy{XT = z)> c^'^ Fy {Xt = zq) > —{2h + 1)^-"^. □ 



2d 



Lemma 3.4.10. Let h,p be positive integers. Let N be the number of particles that ever 
visit the cube Q{o, p) during the internal DLA process, if one particle starts at each site 
y £TL^ — Q{o, p + h), and k additional particles start at sites yi,. . . ,yk ^ Q{o,p + h). If 
k< Ih'^, then 

N < Binom(A;,p), 
where p < 1 is a constant depending only on d. 
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Proof. Let Fj be a half-space defining a face of the cube Q = Q{o, p), such that dist{y j, Fj) > 
h. Let Zj G Fj be the closest point to yj in Fj, and let Bj = Q{zj, h/2) n Fj. Let be the 
random set of sites where the particles starting at yi, . . . ,yj stop. Since #Bj > ^h'^ > 2k, 
there is a hyperplane Hj parallel to Fj and intersecting Bj, such that 

#Bj n Hj n < ^h'^'\ (3.85) 

Denote by Aj the event that the particle starting at yj ever visits Q. On this event, the 
particle must pass through Hj at a site which was already occupied by an earlier particle. 
Thus if {Xt}t>o is the random walk performed by the particle starting at yj, then 

A, c {Xt G Zj.i}, 



where T is the first hitting time of Hj. By Lemma 3.4.9 every site z ^ BjT] Hj satisfies 



¥{Xt = z)> ah^-'^. 



From ( |3.85[ ), since #Bj n Hj > we obtain F{Xt ^ Zj^i) > a/2, hence 

F{Aj\J^j.i) < p 

where p = 1 — a/2, and J^i is the o"- algebra generated by the walks performed by the particles 
starting at yi, . . . ,yi. Thus we can couple the indicators with i.i.d. indicators Ij > 
of mean p to obtain 

k k 
iV = ^ 1^^. < ^ Ij = Bmom{k, p) . □ 
j=i i=i 

The next lemma shows that if few enough particles start outside a cube Q(o, 3p), 
it is highly unlikely that any of them will reach the smaller cube Q{o,p). 

Lemma 3.4.11. Let p, k be positive integers with 

k<\{l-p'n'p' 



where p < 1 is the constant in Lemma 3.4-10. Let N be the number of particles that ever 



visit the cube Q{o, p) during the internal DLA process, if one particle starts at each site 
y — Q{o, 3/9), and k additional particles start at sites yi, . . . ,yk ^ Q{o, 3/3). Then 

F{N > 0) < coe-^i'' 

where cq,ci > are constants depending only on d. 
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Proof. Let Nj be the number of particles that ever visit the cube Qj = Q{o, pj), where 
Let kj = jpf^k, and let Aj be the event that Nj < kj. Taking h = pj — pj^i in Lemma 



3.4.10 



since 



we obtain 



Hence 



Nj+ilA, < Bmom{Nj,p). 



;+i\Aj) > P(Binom(A;j,p) < kj+i) 
> 1 - 26-"^^ 



(3.86) 



where in the second line we have used Lemma 3.4.7 with A = — P- 
Now let 



log k — log p 



log(l/p) 

so that p^^'^p < kj < p. On the event Aj, at most p particles visit the cube Q{o,2p). Since 
the first particle to visit each cube Q{o, 2p — i) stops there, at most p — i particles visit 



Q{o,2p — i). Taking i = p we obtain ¥{N = 0) > ¥{Aj). From (3.86) we conclude that 

P(Ar = 0) > F{Ai) F {A2\Ai) {Aj\Aj^i) > 1 - 2je~^^. 
The right side is at least 1 — cqc^^^'^ for suitable constants cq, ci. 



□ 



Proof of Theorem 3.4-1 The inner estimate is immediate from Lemma 3.4.8, For the outer 
estimate, let 

Nn = #{l < i < mn\v' > f'} = m„ - 
be the number of particles that leave D ' before aggreg ating to the cluster. Let Kq = 



3.4.10 



I { ^^^ — ) where p is the constant in Lemma 
the event 

{Nn<KQ{e/6nfioTa\\n>no]. 



For e > and uq >\, consider 



By (3.63) and (3.64), the closure of D contains the support of a, so by Proposition 
(5^m„ = 5^ a-aix) ^ / a{x)dx = C{b). 



3.1.13 



(3.87) 
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Moreover, by Proposition 3.1.13[ i), for sufficiently small rj we have 



C{D-D2rj) < -i^oe • 
Taking n large enough so that 5„ < rj, we obtain 

Si#D-:^ = C{bf) > C{D2r,) > C{D) - ^-Koe^ 
Thus for sufficiently large n, on the event D':^ C In we have 



n 



dx~d 



<mn-5-^C{D) + -Koe%~ 



where in the last line we have used (3.87). From Lemma 3.4.8 we obtain 



{Fno) > ^{D'^ C in for all n > no) T 1 



(3.88) 



as no t oo. 

By compactness, we can find finitely many cubes Qi, . . . , Qm of side length p = 
2e/3\/d centered at points in d{D^), with d{D^) C \JQi. Taking k = [Ko{e/5n)'^\ in 
Lemma 



3.4.11 



since 3Qi is disjoint from D, we obtain for n > no 



Summing over i yields 

F{{D' n /„ / 0} n Fno) < come-'^P. 
By Borel-Cantelli, if G is the event that /„ ^ D''''' for infinitely many n, then P(F„,gnG) = 0. 



From (3.88) we conclude that P(G) = 0. 



□ 



3.5 Multiple Point Sources 



This section is devoted to proving Theorems 1.3.1 and 1.4.1 
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3.5.1 Associativity and Hausdorff Continuity of the Smash Sum 



In this section we establish two basic properties of the smash sum (1.6) that are 



needed to prove Theorems 1.3.1 and 1.4.1 In Lemma 3.5.1 we show that the smash sum is 



associative, and in Lemma 3.5.3 we show that it is continuous in the Hausdorff metric. 

Lemma 3.5.1. Let A,B,C C M'^ be bounded open sets whose boundaries have measure 
zero. Then 

{A®B)®C = A®{B®C) 

= AUBUCUD 



where D is given by (3.10) with a = 1a + + Re- 
proof. Let 



7(x) 



G(U + 1b)(x) 



and 



7(x) = -\x\'^ -G{Ia®b + Ic){x). 



Let ti = s — 7 and u = s — 7, where s, s are the least superharmonic majorants of 7, 7. Then 

{A®B)®C = {A®B)\JC\J{u>{)} 

= ^ U S U {n > 0} U C U > 0} 

= ylU5uCu{M + u>0}. (3.89) 

Let Vn be the final mass density for the divisible sandpile in S'li^j started from source density 



1a-' + 1_B'-) ^-iid let Un be the corresponding odometer function. By Theorem 3.2.1 we have 
ti„ — > li as n — > cxD. Moreover, by Theorem 3.2.7 we have 



(3.90) 



for all X ^ d{A © B). Let Un be the odometer function for the divisible sandpile on 5n'^'^ 
started from source density Vn + Ic--- By Proposition 3.1.13[ i) the right side of (3.90) 
is continuous almost everywhere, so by Theorem 3.2.1 we have n„ ^ n as n 



00. On 



the other hand, by the abelian property. Lemma 2.2.1 the sum n„ + Un is the odometer 
function for the divisible sandpile on JnZ'^ started from source density Ia- + 1b- + lc--i so 

-u := s — 7, where 

|2 



by Theorem 3.2.1 

7(x) = - 



\x\ 



G{1a + Ib + 1c){x) 
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and s is the least superharmonic majorant of 7. In particular, 



U + U = lim {Un + Un) = "U, 



so the right side of (3.89 ) is equal io AU B U C U D. 



□ 



The following lemma shows that the smash sum of two sets with a small intersection 
cannot extend very far beyond their union. As usual, C denotes Lebesgue measure in W^, 
and A"" denotes the outer e-neighborhood of a set A C W^. 

Lemma 3.5.2. Let A,Bc M'^ be bounded open sets whose boundaries have measure zero, 
and let p = {C{A n B))^/'^. There is a constant c, independent of A and B, such that 

A®B C [AuByP. 



Proof. Let 7,5 be given by (3.8) and (3.9) with a = 1a + 1b, and write u = s — 7. Fix 

X 



£{A®B)-{A\JB) and let r = dist(x, AUS). Let B = B{x,r/2). By Lemma |3T2| iii) 



s is harmonic m. A (B B — {A[J B) , so by Lemma 3.1.4 the function 



w{y) = u{y) -\x- 

= s(y) + + G{Ia + lB){y) -\x- y\^ 

is harmonic on the intersection {A®B)r\ B; hence it attains its maximum on the boundary. 
Since w{x) > the maximum cannot be attained on d{A (B B), so it is attained at some 
point y G dB, and 

(3.91) 



u{y) > w{y) + ^ > ^ 



If z is any point outside A (B B, then by Lemma 3.1.6 and Lemma 3.1.19 with 
A = Ad, there is a constant c' > 1 such that u < c'h^ on B{z, h) for all h. Taking 



h = r/2vc', we conclude from (3.91) that B{y,h) C A(B B. Since B{y,h) is disjoint from 



Au B, we have by Corollary 3.1.14 



C{A UB)+ ujdh'^ < C{A ®B)= C{A) + C{B), 
hence ujdh'^ ^ H B). Taking c = 2\/c' /ij^]J'^ yields r = cuh/'^h < cp. 



□ 



The following lemma, together with Lemma 3.1.8 shows that the smash sum 
operation is continuous in the Hausdorff metric. 
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Lemma 3.5.3. Let A,Bc M"^ be bounded open sets whose boundaries have measure zero. 
For any e > there exists rj > such that 



{A e B), C ® C ® B"^ C {A® B) 



(3.92) 



Proof. By Lemma 3.1.20, since A (B B = AU B D D, we have 

{AeB)e C Ae' U B^> U A' 



for some e' > 0. By Lemma 3.1.17lii) with a = 1a + and (7„ = + l-Bi/„i for 



3.5.1 



sufficiently small r] we have D^' C © -6,^. This proves the first inclusion in (3.92). 
The second inclusion is immediate from Lemma 13.1.81 

For the final inclusion, write A' = A^ - A and B' = B'' - B. Since C{A') j 
C{dA) = and C{B') | C{dB) = as | 0, for small enough r] we have by Lemmas 
and [3X2 

AiQB"^ = A®B®A' ®B' 

C {{AeB)uA''UB''Y~'^ 

c{AeBY. □ 



3.5.2 Smash Sums of Balls 



In this section we deduce Theorem 1.3.1 from our other results. The following 
result on internal DLA with a single point source is a restatement of the main result of ^29j . 



Theorem 3.5.4. Fix A > 0, and let In be the random set of occupied sites for internal 
DLA in 6'ii^^ ^ starting with m — ['^^n '^J particles at the origin. If Sn ^ 1 /ri for all then 
for any e > 0, we have with probability one 

B,; C In C B""" for all sufficiently large n, 

where B is the ball of volume A centered at the origin in M"'. 



Proof of Theorem L3.1. For i = 1, . . . , A; let D^, Rn, In be the domain of occupied sites in 
starting from a single point source of |_(5~'^AjJ particles at x'-^. By Theorems 



1.2.1 



and 1.2.2 for any > we have 

{BiY^ C D'n, Rn C {BiY'-'- for all i and ah sufficiently large n. (3.93) 
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Moreover if 5^ ^ l/''^, then by Theorem 3.5.4 we have with probabihty one 



(-Bj);^ C In C [BiY" for all i and all sufficiently large n. 



(3.94) 



Next we argue that the domains Dn-,Rn-,In can be understood as smash sums of 
D^^ Rl^, In as i ranges over the integers 1, . . . ,k. By the abelian property fl5|, the domain 
/„ is the Diaconis-Fulton smash sum in 5n'^'^ of the domains Likewise, if ri is an 
arbitrary rotor configuration on 6n'^'^, let be the smash sum of and R^ formed using 
rotor-router dynamics with initial rotor configuration ri, and let r2 be the resulting final 
rotor configuration. For i > 3 define inductively as the smash sum of S!^^ and R^ 
formed using rotor-router dynamics with initial rotor configuration rj_i, and let rj be the 
resulting final rotor configuration. Then i?„ = S^- Finally, by Lemma 



2.2.1 



the domain 



Dn contains the smash sum of domains formed using divisible sandpile dynamics, and 
Dr, is contained in the smash sum of the domains Dl, U dDl, . 



Fixing e > 0, by Theorem 1.1.3 and Lemma 3.1.8 it follows from (3.93) that for 
all sufficiently large n 

^•/2 CL>„,i?„,/„ (3 95) 

where A is the smash sum of the balls {Bi)^, and A is the smash sum of the balls {Bi)"^ . By 



Lemma 3.5.3 we can take r] sufficiently small so that 



A/2 CACAC Z?^/2 



where D = Bi 



Bk- Together with (3.95), this completes the proof. 



□ 



3.5.3 Algebraic Boundary 



In this section we prove Theorem 1.4.1 showing that the boundary of a smash sum 
of disks in M? is an algebraic curve. The first step is to establish the classical quadrature 



identity (1.12), which holds in any dimension. We will prove the following quadrature 
inequality generalizing (1.12) for superharmonic functions on a smash sum of balls in M'^. 



Proposition 3.5.5. Fix xi, . . . ,Xk G M'^ and Ai, . . . , > 0. Let Bi be the ball of volume 
Xi centered at Xj. Then 



j 

JBi 



3iw...te-Bfc 

for all integrable superharmonic functions h on Bi® . 



h{x)dx < \ih{xi) 



Bh. 
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We could deduce Proposition 3.5.5 from Proposition 3.1.11 by integrating smooth 
approximations to h against smooth densities 0"„ converging to the sum of the indicators of 
the Bi. However, it more convenient to use the following result of M. Sakai, which is proved 
by a similar type of approximation argument. 

Theorem 3.5.6. ;42, Thm. 7.5] Let ^ CW^ be a bounded open set, and let a be a bounded 
function on M"^ supported on i7 satisfying a > 1 on 17. Let 7, s,D be given by (3.8)-(3.10). 



h{x)dx < I h{x)a{x)dx 
D J D 

for all integrable superharmonic functions h on D. 

We remark that the form of the obstacle in |42J is superficially different from ours: 
taking w = 1 and uo = a \n section 7 of |l2], the obstacle is given by 



i)(x) = Gl 



B 



Ga 



for a large ball B, rather than our choice (3.8) of 

7(3;) = —\x\'^ — Ga. 



Note, however that — ip is constant on B by (3.6) and (3.7). Thus if B is sufficiently large 



the two obstacle problems have the same noncoincidence set D by Lemmas [3.1.9 and 3.1.16 



To prove Proposition 3.5.5 using Theorem 3.5.6 we must produce an appropriate 
density a on M'^ strictly exceeding 1 on its support. We will take a to be twice the sum of 
indicators of balls of half the volume of the B^. 



Proof of Proposition 3.5.5 Let B'- be the ball of volume Aj/2 centered at Xi, and consider 
the sum of indicators 



a = 2j]l^,. 



i=l 



Let 7,8,1? be given by (3.8)-(3.10). By Theorem 3.5.6 and the mean value proprerty for 



superharmonic functions (3.1), we have 



/I k n k 

/ h{x)dx<2j2 h{x)dx<Y^ 
•^^ i=i -^^'i i=i 



XihiyXi 



for all integrable superharmonic functions h on D. 



98 



It remains to show that D = Bi 



Bk- By Lemma 3.1.7 we have B[ d D for 



all i, hence by Lemma 3.5.1 



D = B[®B[ 



By Lemma 3.1.15 we have Bi = B'- (B B'-, completing the proof. 



□ 



Theorem 1.4.1 now follows from Proposition 3.5.5 by a result of Gustafsson 



sec. 6]; see also |23| Lemma 1.1(a)]. The methods of [i22j rely heavily on complex analysis, 
which is why they apply only in dimension two. Indeed, in higher dimensions it is not 
known whether classical quadrature domains have boundaries given by algebraic surfaces 
Ch. 2]. 
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Chapter 4 



Aggregation on Trees 



4.1 The Sandpile Group of a Tree 

In this section we describe a short exact sequence relating the sandpile group of 
a tree to those of its principal subtrees. In the case of a regular tree this sequence splits, 
enabling us to compute the full decomposition of the sandpile group into cyclic subgroups. 
This resolves in the affirmative a conjecture of E. Toumpakari concerning the ranks of the 
Sylow p-subgroups. 

We begin with a simple combinatorial problem. Let Tn be the d-regular tree of 
height n. Collapse all the leaves of T„ to a single vertex s, the sink, and add an edge 
connecting the root to the sink. 

Lemma 4.1.1. Let tn be the number of oriented spanning trees of rooted at the sink. 
Then 

tn = ttl{dtn-l-{d-l)ttl). 

Proof. If the edge (r, s) from the root to the sink is included in the spanning tree, then each 
of the principal branches of T„ may be assigned an oriented spanning tree independently, 
so there are t'^\ such spanning trees. On the other hand, if (r, s) is not included in the 
spanning tree, there is a directed path r — > xi — > . . . — > Xn-i — > s in the spanning tree 
from the root to the sink. In this case, every principal branch except the one rooted at 
xi may be assigned an oriented spanning tree independently; within the branch rooted at 
xi, every subbranch except the one rooted at X2 may be assigned an oriented spanning 



tree independently; and so on (see Figure 4.1). Since there are {d — 1)" possible paths 
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Figure 4.1: The two cases in the proof of Lemma 



4.1.1 



n-l 



.d-2 
''fc ■ 



k=l 



j;i — > . . . — > Xn-i, we conclude that 

tn = tt\ + {d-ir-'Yltt 

Substituting n — 1 for n we find that 



n-2 



k=l 



hence from (4.1 ) 



tn = ttlitn-1 + {d- l)(t„_i - ^2))- □ 



(4.1) 



From Lemma 4.1.1 one can readily show by induction that 

n-2 

t^ = {l + a + ... + a'^-i) J](l + a + . . . + a'^)"" 

k=l 



where a = d — 1. A variant of this formula was found by Toumpakari [57] , who gives an 
algebraic proof. 

For any graph G there is an abelian group, the sandpile group, whose order is the 
number of oriented spanning trees of G rooted at a fixed vertex; its definition and properties 
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are reviewed in section 4.1.1 A product formula such as the one above immediately raises 



the question of an analogous factorization of the sandpile group. Our main goal in this 



section is to prove Theorem 1.5.2, which establishes precisely such a factorization. In 



section 4.2 we will apply this result to study the rotor-router model on regular trees. 



The next three sections are organized as follows. In section 4.1.1 we review the 
definition and basic properties of the sandpile group, and characterize the recurrent states 
on a tree explicitly in terms of what we call critical vertices. We prove a general result, 



Theorem [4.1.5 , relating the sandpile group of an arbitrary tree to the sandpile groups of its 



principal branches (i.e. the subtrees rooted at the children of the root). This result takes 



the form of an isomorphism between two quotients. In section 4.1.2| we take advantage of 
the symmetry of regular trees to define a projection map back onto the subgroup being quo- 
tiented. This allows us to express SP(Tn) as the direct sum of a cyclic group and a quotient 



of the direct sum of d — 1 copies of SP{Tn-i), which enables us to prove Theorem 1.5.2 by 



induction. Finally, in section 4.1.3 we prove Toumpakari's conjecture [57 



4.1.1 General Trees 

Let T be a finite rooted tree. Collapse all the leaves to a single vertex s, the sink, 
and add an edge connecting the root to the sink. If the vertices of T are xi, . . . ,Xn = s, the 
sandpile group of T is defined by 

SP{T) = iTj^ 

where A C is the lattice 

Here A^;- is the vector in Z" taking value 1 at each neighbor of Xj, value — deg(x) at x, and 
value elsewhere. The sandpile group of a graph was originally defined in [5l EJ [Ti] . 

A nonnegative vector u G may be thought of as a chip configuration with Ui 
chips at vertex Xj. A vertex x 7^ s is unstable if u{x) > deg(x). An unstable vertex may 
topple, sending one chip to each neighbor. Note that the operation of toppling the vertex x 
corresponds to adding the vector A^. to u. We say that u is stable if no vertex is unstable. 
Given chip configurations u and v, we define u + v as the stable configuration resulting from 
starting with n(x) + v{x) chips at each vertex x and toppling any unstable vertices; the 
order in which topplings are performed does not affect the final configuration, as observed 
by Dhar |14] and (in a more general setting) Diaconis and Fulton |15] . 
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A chip configuration u is recurrent [TF il u + v = u for a nonnegative configu- 
ration V. In [13j it is proved that every equivalence class mod A has a unique recurrent 
representative. Thus the sandpile group SP{T) may be thought of as the set of recurrent 
configurations under the operation of addition followed by toppling. Note that if f is a 
nonnegative configuration, its recurrent representative is given by 

V := V + e 

where e is the identity element of SP{T) (the recurrent representative of 0); indeed, i) is 
recurrent since e is recurrent, and v = v (mod A) since e G A. Note that if u is a recurrent 
configuration and v is a nonnegative configuration, then 

u + v = u + {v + e) = {u + e)+v = u + v. (4.2) 



The following result is a simple variant of the "burning algorithm" [M]; see also 
da Cor. 2.6]. 

Lemma 4.1.2. Let (3{x) he the number of edges from x to the sink. A chip configuration 
on T is recurrent if and only if adding (3{x) chips at each vertex x causes every vertex to 
topple exactly once. 

Proof. Note that 

/3 = deg(s)s + A, = deg(s)s - ^ A^. (4.3) 
If every vertex topples exactly once, then 

u + /3 = n + /? + ^A^ = u, 

so u is recurrent. Conversely, suppose u is recurrent. Since /3 G A we have /? = e, hence 



from (4.2) 



u + l3 = u + e = u. 



By (4.3), since {/S.x}x^s are linearly independent, every vertex topples exactly once. □ 



We first characterize the recurrent configurations of a tree explicitly. The charac- 
terization uses the following inductive definition. Denote by C{x) the set of children of a 
vertex x G T. 
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Figure 4.2: A recurrent configuration on the ternary tree of height 5. Critical vertices 
are circled; if any of the circled vertices had fewer chips, the configuration would not be 
recurrent. 



Definition. A vertex x G T is critical for a chip configuration u x ^ s and 

u{x) < #{y G C{x) I y is critical}. (4.4) 
Proposition 4.1.3. A configuration u G SP(T) is recurrent if and only if equality holds 



in (4-4) for every critical vertex x. 
Proof. If X is critical, then 

u{x) + #{y G C(x) I y is not critical} < deg(x) — 1. 



(4.5) 



Thus after chips are added as prescribed by Lemma 4.1.2 inducting upward in decreasing 
distance to the root, if x 7^ r is critical, its parent must topple before it does. In particular. 



if strict inequality holds in (4.4), and hence in (4.5), for some vertex x, that vertex will 



never topple, so u is not recurrent. 



Conversely, suppose equality holds in (4.4), hence in (4.5), for every critical x 



Begin toppling vertices in order of decreasing distance from the root. Note that a non- 
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critical vertex x satisfies 

u{x) + £ C{x) I y is not critical} > deg(x), (4-6) 



Inducting upward, every non-critical vertex topples once. Hence by equality in ( |4.5| ), once 
all vertices other than the root are stable, every critical vertex x has either toppled (if its 
parent toppled) or is left with exactly deg(x) — 1 chips (if its parent did not topple). In 
particular, the root now topples, as it was given an extra chip in the beginning. Now if 
a; is a critical vertex that has not yet toppled, its parent is also such a vertex. Inducting 
downward from the root, since all of these vertices are primed with deg(x) — 1 chips, they 
each topple once, and u is recurrent. □ 

Write Ti, . . . ,Tk for the principal branches of T (i.e. the subtrees rooted at the 
children of the root). If Ui is a chip configuration on Ti, and a is an integer, we will use the 

notation | | for the configuration on T which has a chips at the root and coin- 

Ul,...,Uk 



cides with m on Ti . The following result is an immediate consequence of Proposition |4.1.3 

( " 

Lemma 4.1.4. Let u = I 

Y Ml, . . . ,nfc 

(i) If u is recurrent, each Ui is recurrent. 

(a) If ui, . . . ,Uk are recurrent and a = k, then u is recurrent. 

Write 6x for a single chip at a vertex x, and denote hy x = e + dx the recurrent 



form of 6x. Note that by (4.2), if u is recurrent then 

u + X = u + 6x. (4-7) 
Theorem 4.1.5. Let Ti, . . . ,Tk be the principal branches ofT. Then 

k 

SP(T)/{r)^^SP{Ti)/{{h,...,h)) 

i=l 

where r, ri are the roots ofT, Ti respectively. 
Proof. Define 4) : SP{T) 0^^^ SP{Ti) by 

^ {ui,...,Uk). (4.8) 
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Lemma 4.1.4[ i) ensures this map is well-defined. Note that if 



+ r, 



Ul,...,Uk 



Vl,...,Vk 



by (4.7) either b < k and Uj = Vi for all i; oi b = k and the root topples, in which case Ui 



Vi+Ti. Thus 4> descends to amap of quotients : SP{T)/{f) ©j=i SP{Ti) / {{fi, . . . ,ffc)). 

By adding two configurations without allowing the root to topple, the configura- 



tions on each branch add independently, hence by (4.7) 




-|- dr 



VI,..., Vk 



Ul + Vl, . . . ,Uk + Vk 



for some nonnegative integers c, d. Thus (p is a group homomorphism. Moreover, 
surjective by Lemma 4.1.4[ ii). Finally, if 



IS 



{ui,...,Uk) = {vi,...,Vk) +c{ri,---,rk), 



then by (4.7), allowing the root to topple exactly c times, we obtain 




+ {c{k + l) + d)f. 



Vl,...,Vk 



for a suitable integer d. Thus (j) is injective. 



□ 



4.1.2 Regular Trees 



In this section we show that for regular trees. Theorem 4.1.5 can be strengthened 
to express SP{T) as a direct sum. 

Let T„ be the regular tree of degree d and height n, with leaves collapsed to the 
sink vertex and an edge added from the root to the sink as in section |4.1.1[ The chip 
configurations which are constant on the levels of T„ form a subgroup of SP{Tn). If each 
vertex at height k has chips, we can represent the configuration as a vector (oi, . . . , a„_i). 
If such a recurrent configuration is zero on a level, all vertices above that level are critical, so 



by Proposition 4.1.3 they must have d — 1 chips each. The recurrent configurations constant 
on levels are thus in bijection with integer vectors (oi, . . . , a„_i) with < a, < d — 1 subject 
to the constraint that if = then ai = . . . = aj_i = d — 1. 
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f 2f 3f 4f 5r 6r 7r 



3r lOr llf 12r 13r 14f 15r 



2012012220 1 2 1 2 

0111222201 1 1 2 2 2 

2222222011 1 1 1 1 1 

Figure 4.3: Multiples of the root f in the ternary tree of height 4. Each column vector 
represents a chip configuration which is constant on levels of the tree. 



The following lemma uses the lexicographic order given by a < 6 if for some k 
we have On-i = bn-i, ■ ■ ■ ,ak+i = bk+i and < bk- In the cyclic lexicographic order on 
recurrent vectors we have also (d — 1, . . . , d — 1) < (d — 1, . . . , d — 1, 0). 

Lemma 4.1.6. Ifu,v are recurrent configurations on that are constant on levels, write 
u V if V follows u in the cyclic lexicographic order on the set of recurrent vectors. Then 
for every integer k >0, we have 

kr {k + l)f. 



Figure 4.3 demonstrates the lemma for a ternary tree of height 4. 



Proof. By (4.7) we have 



{k + l)f = kr + Sr- 

Thus if A;r = (ai, . . . , a„_i) with ai < d—1, then (k + l)r = (ai + 1, 02, 
Otherwise, if not all ai equal d — 1, let j > 1 be such that ai = ... 
Uj < d — 1. Adding a chip at the root initiates the toppling cascade 



as desired. 
■ d—1 and 



/ d \ 

d-1 
d-1 

d-1 
d-1 
aj 



( 





d 

d-1 

d-1 
d-1 
aj 



\ 



fd-l\ 

d 

d-1 
d-1 



( d-l\ 

d-1 
d-1 



d 

\^ a-n-i J 



( d-1 \ 

d-1 
d-1 

d-1 


+ 1 
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as desired. If all = d — 1 the cascade will travel all the way down, ending m {d — 1, . . . ,d — 
1,0) as desired. □ 

Proposition 4.1.7. Let Tn be the regular tree of degree d and height n, and let R(Tn) be 
the subgroup of SP(Tn) generated by f. Then R{Tn) consists of all recurrent configurations 
that are constant on levels, and its order is 

*R{Tn) = (4.9) 

Proof. Since the toppling rule is symmetric, all configurations in R{Tn) are constant on 
levels. The number of such recurrent configurations is the number of vectors of the form 
(d — 1, . . . , d — 1, 0, flj, . . . , a„_i), with G [d — 1], which is 



Moreover, by Lemma 4.1.6 any such vector can be expressed as a multiple of f, so R{T.n) 



contains all the recurrent configurations that are constant on levels. □ 

Index the vertices of the d-regular tree of height n by words of length < n — 2 in 
the alphabet {1, . . . ,d — 1}. Let be the automorphism of the tree given by 

ai{wi . . . Wk) = wi . . . {wi + 1) . . . Wk 

with the sum taken mod d — 1; k < i then ai{w) = w. Given a map a : [n — 2] — > [d — 1] 
let (Jq be the composition 13?=!^ ^t^^^ • 

If a is an automorphism of the form (Jq, write au for the chip configuration au{x) = 
u{ax). Writing u (B v for addition in the sandpile group and u + v for the ordinary vector 
sum, we have 

m 

u®v = u + v + ^ A^^. 

j=i 

Since aA^ = A^-x we obtain 

m 

a{u ev) = o-(m) + a{v) + ^ A^^^ . 

i=i 

The configuration on the right side is stable, recurrent, and = (t{u) + cr{v) (mod A), so it 
is equal to cj(u) © cr(v). Thus a is an automorphism of the sandpile group. 
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Figure 4.4: A non-regular tree for which Proposition 



4.1.8 



fails. 



Proposition 4.1.8. Let Tn he the regular tree of degree d and height n, and let R{Tn) = (r) 
he the suhgroup of SP(Tn) generated hy the root. Then 

with d — 1 summands of SP{Tn-i) on the right side. 
Proof Define p : SP{Tn) SP{Tn) by 



p{u) = {d- 1)^ ^ aaU. 

a:[n~2]-*[d-l] 



(4.10) 



By construction p{u) is constant on levels, so the image of p lies in R{Tn) by Proposi- 
tion 4.1.7 Given u £ R(Tn), since u is constant on levels we have ^"^(n) = u for all a. Since 
there are {d — 1)"~^ terms in the sum (4.10), we obtain 



p{u) = {d- l)"n = u 



where the second inequality follows from (4.9). Thus R(Tn) is a summand of SP(T„), and 
the result follows from Theorem I4.1.5[ □ 



Proposition 4.1.8 fails for general trees. For exai nple , if T is the tree consisting of a 

2 



root with 2 children each of which have 3 children (Figure 



4.4|), then r = | | has order 

3,3 
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, 2 , 

10 and the element x = \ \ satisfies 4x = f, so x has order 40. The total number of 

0,3 

recurrent configurations is 4 • 4 • 3 - 8 = 40, so SP{T) ~ Z/40Z, and R{T) ~ Z/IOZ is not 
a summand. 

Write Zp as a shorthand for (Z/pZ) © ... © (Z/pZ) with g summands. 

Theorem 4.1.9. Let Tn he the regular tree of degree d = a + 1 and height n, with leaves 
collapsed to the sink vertex and an edge joining the root to the sink. Then 

Proof. Induct on n. For the base case n = 2 we have 

SP{T2) = R{T2) ^'Ld = Zl+a. 



Write g„ = l + ffl + -- - + a" • By Proposition 4.1.7 the root subgroup R{Tn) is cyclic of 



order qn. By Proposition 4.1.8 and the inductive hypothesis, it follows that 



SP{Tn ^^^^ 



7a— 1 



"("^-^^ © ... ©zf '(""^^ □ 



-'qn ^ "-'q„-l ^ qn-2 VLV . . . VLV 



4.1.3 Proof of Toumpakari's Conjecture 

As before write a = d — 1 and 

qn = l + a + ... + a 



n-l 



If p is a prime not dividing d{d — 1), let tp be the least positive n for which p\qn- Then 

p if a = 1 (mod p) 

ordp(a), else. 

Here ordp(a) is the least positive k for which p\a'' — 1. Note that p\qn if and only if tp\n. 
The following result was conjectured by E. Toumpakari in [U] (where the factor of d — 2 
was left out, presumably an oversight). 
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Theorem 4.1.10. Let Bn he the hall of radius n in the d-regular tree, with each leaf con- 
nected hy d — 1 edges to the sink, hut with no edge connecting the root to the sink. Let p he 
a prime not dividing d{d — 1), and let Sp{n) he the Sylow-p suhgroup of the sandpile group 
SP{Bn). Then 



rank(S'p(n)) 



d{d-2)Y^ m<n {d-l) 

m,=n (mod tp) 



if n ^ — 1 (mod tp); 



d{d-2)J2 m<n (d-l)™ + d-l, if n = -l(mod tp). 

m=n (mod tp) 



Proof. By Theorem 4.1.5 we have 



SP{Bn)/{r) 



SP{Tr 



n+1, 



SP{T, 



n+l) 



{R{Tn+l), • • • , R{Tn+l)) 



4.1.7 



we have R{Tn+i] 



with d summands. By Proposition 

sp{Bn)/{r) ~ z;;^^^ ezl^'-i^^^+i) ezl'^"!)''(''+^) 



q^j^-^, SO from Theorem 4.1.9 



(4.11) 



By Proposition 7.2 of [47 , the root subgroup (f) of SP{Bn) has order d{d — 1)". Thus for 
p not dividing d{d — 1) the Sylow p-subgroup of SP{Bn) is the same as that of the quotient 



SP{Bn) /{r). Each summand Z^^. in (4.11) contributes 1 to the rank of Sp{n) if tp\k and 
otherwise. If n ^ — 1 (mod tp), the total rank is therefore 

rank{Sp{n))= ^ (a - l)a"-*^(a + 1) 

2<k<n 
tp I k 

= d{d-2) (^-ir- 

0<m<n-2 
m=n (mod tp) 



In the case that n = — 1 (mod tp), the first summand Zg,^^-^ in (4.11) contributes an addi- 
tional rank a = d — 1 to Sp{n). □ 



4.2 The Rotor-Router Model on Trees 



This section is devoted to proving Theorems 1.5.1 and 1.5.3 



4.2.1 The Rotor-Router Group 

In this section we define the rotor-router group of a graph and show it is isomorphic 



to the sandpile group. This isomorphism, Theorem 4.2.5 is mentioned in the physics 
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literature; see |38l [37] . To our knowledge the details of the proof are not written down 
anywhere. While our main focus is on the tree, the isomorphism is just as easily proved for 
general graphs, and it seems to us worthwhile to record the general proof here. 

Let G be a strongly connected finite directed graph without loops. Fix a sink 
vertex s in G, and write Rec{G) for the set of oriented spanning trees of G rooted at the 
sink. Given a configuration of rotors T, write ex{T) for the configuration resulting from 
starting a chip at x and letting it walk according to the rotor-router rule until it reaches 
the sink. (Note that if the chip visits a vertex infinitely often, it visits all of its neighbors 
infinitely often; since G is strongly connected, the chip eventually reaches the sink.) We 
view T as a subgraph of G in which every vertex except the sink has out-degree one. Note 
that such a subgraph is an oriented spanning tree rooted at the sink if and only if it has no 
oriented cycles. 

Lemma 4.2.1. //T G Rec{G), then e^iT) G Rec{G). 

Proof. Let Y be any collection of vertices of G. If the chip started at x reaches the sink 
without ever visiting Y, then the rotors at vertices in Y point the same way in ex(T) as 
they do in T, so they do not form an oriented cycle. If the chip does visit Y, let y G y be 
the last vertex it visits. Then either y = s, or the rotor at y points to a vertex not in Y; in 
either case, the rotors at vertices in Y do not form an oriented cycle. □ 

We will need slightly more refined information about the intermediate states that 
occur before the chip falls into the sink. These states may contain oriented cycles, but only 
of a very restricted form. For a vertex x we write CyCx{G) for the set of rotor configurations 
U such that 

(i) U contains an oriented cycle; and 

(ii) If the rotor U{x) is deleted, the resulting configuration contains no oriented cycles. 

Lemma 4.2.2. Starting from a rotor configuration Tq G Rec{G) with a chip at xq, let 
and Xj. he the rotor configuration and chip location after k steps. Then 

(i) Ifn i Rec{G), then n e CyCx,{G). 

(ii) IfTk G Rec{G), then Xk ^ {xq, . . . ,Xk-i}. 
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Proof, (i) It suffices to show that any oriented cycle in contains Xk- Let Y be any set 
of vertices of G not containing Xk- If Y is disjoint from {xq, . . . , x^.i}, then the rotors at 
vertices in Y point the same way in as they do in Tq, so they do not form an oriented 
cycle. Otherwise, let y G y be the vertex visited latest before time k. The rotor at y points 
to a vertex not in Y, so the rotors at vertices in Y do not form an oriented cycle. 

(ii) Suppose Xk G {xq, . . . ,Xk-i}- Let yo = ^k, and for i = 0,1,... let y^+i = 
TkiVi)- Then the last exit from x^ before time k was to yi, and by induction if yi, . . . , yj_i 
are disjoint from x^, then yj_i was visited before time k, and the last exit from yj_i before 
time k was to y^. It follows that y^ = Xk for some i, and hence T^. contains an oriented 
cycle. □ 

Lemma 4.2.3. IfTi,T2 € Rec{G) and e^iTi) = e^(r2), then Ti = T2. 

Proof. We will show that T can be recovered from e^(T) by reversing one rotor step at a 
time. Given rotor configurations U, U' and vertices y, y', we say that ([/', y') is a predecessor 
of (U, y) if a chip at y' with rotors configured according to U' would move to y in a single step 
with resulting rotors configured according to U. For each neighbor z ^ y with U{z) = y 
there is a unique predecessor of the form ([/', z), which we will denote Pz{U,y). 

Suppose {U,y) is an intermediate state in the evolution from T to ex(T). li U ^ 
Rec{G), then by case (i) of Lemma 4.2.2 there is a cycle of rotors U{y) = yi,U{yi) = 
y2, . . . , U{yn) = y. If U{z) = y and z / y„, then z is not in this cycle, so the predecessor 
Pz(U, y) has a cycle disjoint from its chip location. Thus Pz{U, y) does not belong to Rec{G) 



or to GyCz{G), so by Lemma 4.2.2 it cannot be an intermediate state in the evolution from 
T to ex(T). The state immediately preceding {U,y) in the evolution from T to ex(T) must 
therefore be Py^{U,y). 



Now suppose U £ Rec{G). By case (ii) of Lemma 4.2.2 U is the rotor configuration 
when y is ffist visited. If y = x, then U = T. Otherwise, let x = xq — > xi — > . . . ^ = s 
be the path in U from x to the sink. Then the last exit from x before visiting y was to xi. 
By induction, if xi, . . . ,Xj-i are different from y, then Xj_i was visited before y and the 
last exit from xj-i before visiting y was to Xj. It follows that xj = y for some j, and the 
state immediately preceding {U,y) must be Pxj-iiU,y). □ 

Thus for any vertex x of G, the operation Cx of adding a chip at x and routing it 
to the sink acts invertibly on the set of states Rec{G) whose rotors form oriented spanning 
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trees rooted at the sink. It is for this reason that we call these states recurrent. We 
define the rotor-router group RR{G) as the subgroup of the permutation group of Rec{G) 
generated by {e^lzeG- Note that if there are two (indistinguishable) chips on G and each 
takes a single step according to the rotor-router rule, the resulting rotor configuration does 
not depend on the order of the two steps. Thus the operators commute, and the group 
RR{G) is abelian; for a general discussion of this property, which is shared by a number of 
models including the abelian sandpile and the rotor-router, see |15] . 

Lemma 4.2.4. RR{G) acts transitively on Rec{G). 

Proof. Given Ti,T2 G Rec{G), let u{x) be the number of rotor turns needed to get from 
Ti(x) to T2(x). Let v{x) be the number of chips ending up at x if u{y) chips start at each 
vertex y, with rotors starting in configuration Ti, and each chip takes a single step. After 
each chip has taken a step, the rotors are in configuration T2, hence 

u{x)ex j Ti = ^ v{x)ex j T2. 
^xeV{G) / \x€V{G) / 

Letting g = YyxeViG)i'"'i^) ~ 'vi^))(^x we obtain T2 = gTi. □ 

Theorem 4.2.5. Let G be a strongly connected finite directed graph without loops, let 
RR{G) be its rotor-router group, and SP{G) its sandpile group. Then RR{G) ~ SP{G). 

Proof. Let V be the vertex set of G. Recall [l3] that the sandpile group can be expressed 
as the quotient 

SPiG)^Z^/{s,Ax)x^v 



where s E y is the sink and 



Ax = Y2 y ~ outdeg(x)x. 



Define (j) : ^ RR{G) by 



(piu) = ^ u{x)ex. 

Since adding outdeg(x) chips at x and letting each chip take one step results in one chip at 
each neighbor y ^ x with the rotors unchanged, we have 

outdeg(j:)e^ = Cy. 

y*-x 



114 



Thus (t){^x) = and (j) descends to a map SP{G) RR{G). This map is trivially surjective; 
to show it is injective, by Lemma 4.2.4| we have 



#RR{G) > #Rec{G) = #SP{G), 
where the equality on the right is the matrix-tree theorem |46] . □ 

4.2.2 Aggregation on Regular Trees 

Let T be the infinite d-regular tree. Fix an origin vertex o in T. In rotor-router 
aggregation, we grow a cluster of points in T by repeatedly starting chips at the origin and 
letting them walk until they exit the cluster. Beginning with Ai = {a}, define the cluster 
An inductively by 

An = An^l U {Xn}, 71 > 1. 

where Xn G T is the endpoint of a rotor-router walk started at o and stopped on first exiting 
An-i- We do not change the positions of the rotors when adding a new chip. 

In this section we use the group isomorphism proved in the last section to show 



that An is a perfect ball for suitable values of n (Theorem 1.5.1 ). The proof makes crucial 



use of the calculation of the order of the subgroup of the sandpile group of T„ generated by 



the root, Proposition 4.1.7 



A function H on the vertices of a directed graph G is harmonic if 
H{x) = ^—r^ y H{y) 



outdeg(x) 



for all vertices x. 



Lemma 4.2.6. Let H he a harmonic function on the vertices of G. Suppose chips on G 
can be routed, starting with u(x) chips at each vertex x and ending with v{x) chips at each 
vertex x, in such a way that the initial and final rotor configurations are the same. Then 

H{x)u{x) = H{x)v{x). 

xeV{G) a:eV(G) 

Proof. Let u = uo,ui, . . . ,Uk = v be the intermediate configurations. If Uj+i is obtained 
from Ui by routing a chip from Xj to yi, then 

^ H{x){u{x)-v{x)) = Y,H{x,)-H{yi). (4.12) 

xeV{G) i 
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If the initial and final rotor configurations are the same, then each rotor makes an integer 



number of full turns, so the sum in (4.12) can be written 



i xev{G) 

where N{x) is the number of full turns made by the rotor at x. Since H is harmonic, the 
inner sum on the right vanishes. □ 

Let Tn be the regular tree of degree d and height n, with an edge added from the 
root r to the sink o. Denote by {Xt)t>Q the simple random walk on r„, and let r be the 
first hitting time of the set consisting of the leaves and the sink. Fix a leaf z of T^, and let 

H{x) = P,(X, = z) (4.13) 

be the probability that random walk started at x and stopped at time r stops at z. 

The quantity H{r) can be computed by a standard martingale argument. Recall 
that the process 

Mt = a" 1^*1 

is a martingale, where a = d — 1 and |x| denotes the distance from x to the sink. Since Mt 
has bounded increments and E^r < oo, we obtain from optional stopping 

= E^Mo = W^rMr =p + {l- p)a~"- 

where p = Pr(^r = o). Solving for p we obtain 

P. {Xr = o) = (4.14) 

In the event that the walk stops at a leaf, by symmetry it is equally likely to stop at any 
leaf. Since there are a'^~^ leaves, we obtain from (4.14) 



a"- ^ a"- — 1 

Lemma 4.2.7. Let a = d — 1. If the initial rotor configuration on Tn is acyclic, then 
starting with ^^-5f chips at the root, and stopping each chip when it reaches a leaf or the 
sink, exactly one chip stops at each leaf, and the remaining " chips stop at the sink. 

Moreover the starting and ending rotor configurations are identical. 



116 



Proof. By Theorem 



4.2.5 



and Proposition 



4.1.7 



the element Cr S RR{Tn) has order 



so the starting and ending rotor configurations are identical. Fix a leaf z and let H be 



the harmonic function given by (4.13). By Lemma 4.2.6 and (4.15), the number of chips 
stopping at z is 



1 



a-l 



H{r) = 1. 



Since there are aP' ^ leaves, the remaining — a 



n-1 



a-l 



chips stop at the sink. □ 



Let T be the infinite d-regular tree. The ball of radius p centered at the origin in 



o G T is 



Br, 



{x € T : |x| < p} 



where is the length of the shortest path from o to x. Write 



#Bp = 1 + (a + 1) 



1 



a-l 

As the following result shows, provided we start with an acyclic configuration of rotors, the 
rotor-router aggregation cluster An is a perfect ball at those times when an appropriate 
number of chips have aggregated. It follows that at all other times, the cluster is as close 
as possible to a ball: \i hp < n < 6p+i then Bp C An C -Bp+i- 

Theorem 4.2.8. Let An he the region formed hy rotor-router aggregation on the infinite 
d-regular tree, starting from n chips at the origin. If the initial rotor configuration is acyclic, 
then Ahp = Bp for all p >0. 

Proof. Define a modified aggregation process A'n as follows. Stop the n-th chip when it 
either exits the occupied cluster or returns to o, and let 

^n = ^n-l U {^n} 

where x'n is the point where the n-th chip stops. By relabeling the chips, this yields a time 
change of the original process, i.e. A'^ = ^/(n) for some sequence /(I), /(2), . . .. Thus it 
suffices to show A'^^ = Bp for some sequence ci, C2, . . .. We will show by induction on p that 
this is the case for 



Cp = 1 + (a + 1) ^ 



t=i 



a-l ' 



and that after Cp chips have stopped the rotors are in their initial state. For the base case 
p = 1, we have ci = a + 2 = d+l. The first chip stops at o, and the next d stop at each 
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of the neighbors of o, so A'^^^ = Bi. Since the rotor at o has performed one full turn, it is 
back in its initial state. 

Assume now that ^Cp_i = ^p-i a-^d that the rotors are in their initial acyclic 



state. Starting with Cp — Cp_i chips at o and letting each take one step, there are chips 
at each neighbor of o. Since the rotor at o has performed an integer number of full turns. 



the rotors remain in their initial acyclic state. By Lemma 4.2.7 exactly one chip will stop 
at each leaf z £ Bp — -Bp_i, and the remainder will stop at o. Thus A'^ = Bp. Moreover, by 



Lemma 4.2.7 once all chips have stopped, the rotors are once again in their initial state, 
completing the inductive step. □ 
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4.2.3 Recurrence and Transience 

In this section we explore questions of recurrence and transience for the rotor- 
router walk on regular trees. We aim to study to what extent the rotor-router walk behaves 
as a deterministic analogue of random walk. We find that the behavior depends quite 
dramatically on the initial configuration of rotors. 

Let Tn be the d-regular tree of height n. We collapse the leaves to a sink labeled b 
for boundary, label the origin o and treat it as both the source and a sink. We will examine 
the hitting rates of the two sinks for various configurations of the model, and compare to 
the expected result for random walk. 

Since the origin is a sink, the recurrence and transience of chips on one principal 
branch of the tree are independent of the rotor configurations on the other branches. There- 
fore, as shown in Figure [43] we focus on a single branch Yn of T„. By the root r of Yn we 
will mean the unique node with \r\ = 1. To each rotor direction we associate an index from 
{1, . . . with direction d corresponding to a rotor pointing to the parent vertex. Rotors 
cycle through the indices in order. 

Lemma 4.2.9. Let Yn be a principal branch of Tn, the regular ternary tree of height n. If 
all rotors initially point in direction 1, then the first 2" — 1 chips started at r alternate, the 
first stopping at b, the next stopping at a, the next at b, and so on until all the rotors again 
point in direction 1. 

Proof. Induct on n. For n = 2 the result is obvious as there is only one rotor, which sends 
the first chip in direction 2 to b, the next chip up in direction 3 to o, and the third chip in 
direction 1 to b, at which point the rotor is again in its initial state. 

Now suppose that the lemma holds for Yn~i. Let L and R be the two principal 
branches of Yn- We think of L and R as each having a rotor that points either to b or back 
up to r. The initial state of these rotors is pointing to r. The first chip is sent from the 
origin to R, which by induction sends it to b. Note that the root rotor is now pointing 



towards R, the i?-rotor is pointing to b, and the L-rotor is pointing to r (Figure 4.6 1). We 
now observe that the next four chips form a pattern that will be repeated. The second 
chip is sent directly to o (Figure |4.6[3) and the third chip is sent to L which sends it to b 



(Figure 4.6;). The fourth chip is sent to R, but by induction this chip is returned and then 



it is sent to a (Figure 4.61). Finally, the fifth chip is sent to L, returned, sent to R, and 
through to b (Figure [4.6^). Note that the root rotor is now again pointing towards R, the 
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(c) (d) (e) 




Figure 4.6: The four-chip cycle, which begins after the first chip has been routed to b. 



i?-rotor is again pointing to b, and the L-rotor is again pointing to r. In this cycle of four 
chips, the two branches R and L see two chips apiece. This cycle repeats 2""^ — 1 times, 
and each subtree sees 2^~^ — 2 chips. 

Recall that we sent the first chip to R, so that it has seen a total of 2"^^^ — 1 chips. 
By induction, all the rotors in R are in their initial configuration. We have sent a total of 
2" — 3 chips. The next chip is sent to o, and the last to L which sends it to b. Now L has 
seen 2"~^ — 1 chips so by induction all of its rotors are in their initial configuration. The 
root rotor is pointing towards L, its initial configuration. We have sent a total of 2" — 1 
chips, alternating between b and o, and all of the rotors of Yn are in the initial configuration, 
so the inductive step is complete. □ 



We remark that the obvious generalization of Lemma 4.2.9 to trees of degree d > 3 
fails; indeed, we do not know of a starting rotor configuration on trees of higher degree which 
results in a single chip stopping at o alternating with a string of d — 1 chips stopping at b. 

Consider now the case of the infinite ternary tree T. A chip performing rotor- 
router walk on T must either return to o or escape to infinity visiting each vertex only 
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finitely many times. Thus the state of the rotors after a chip has escaped to infinity is 
well-defined. Therefore, we can properly define R{m) as the number of chips that stop at 
o after m chips have executed the rotor-router walk beginning at r. The following result 
shows that there is a configuration of the tree for which the rotor-router walk behaves as an 
exact quasirandom analogue to the random walk, in which chips exactly alternate returning 
to the origin with escaping to infinity. 

Proposition 4.2.10. Let T he the infinite ternary tree, with principal branches labeled 
y(2)^ Qjj^ y(3) correspondence with the direction indexing of the rotor at the origin. 
Set the rotors along the rightmost branch of Y^^'> initially pointing in direction 2, and all 
remaining rotors initially pointing in direction 1. Let E{m) be the expected number of chips 
that return to the origin if m chips perform independent random walks on T. Let R{m) be 
the number of chips that return to the origin if m chips perform rotor-router walks on T. 
Then \E{m) — R{ra)\ < ^ for all m. 



Proof. Lemma 4.2.9 implies that for the branches Y^'^'^ and Y^'^\ the chips sent to a given 
branch alternate indefinitely with the first escaping to infinity, the next returning to o, 
and so on. Likewise, chips sent to y^^^ will alternate indefinitely with the first returning 
to o, the next escaping to infinity, and so on. Since chips on the full tree T are routed 
cyclically through the branches beginning with Y^'^\ we see that the chips too will alternate 
indefinitely between escaping to infinity and returning to the origin, with the first escaping 



to infinity. Thus R{m) = |_yJ • Taking n ^ oo in (4.14) we obtain E{m) = y, and the 



result follows. □ 

Lemma 4.2.11. Let Y^ be a principal branch of Tn, the d-regular tree of height n. If all 
rotors initially point in direction d — 1, then the first n — 1 chips return to a before hitting 
height n ofY^. 

Proof. Induct on n. For n = 2, there is a single rotor which sends the first chip from r to 
o. Now suppose the lemma holds for l^-i. Let Zi, . . . , Z^-i be the principal branches of 
Yn. The first chip placed at r is sent directly to o. By the inductive hypothesis, the first 
n — 2 chips that are sent to each branch Zi are returned to r before hitting height n of Y^. 
Thus each of the next n — 2 chips is sent to Zi, returned to r, sent to Z2, and so on until 
it is sent to Zii_i, returned to r and then routed to o. □ 



121 



Our next result shows that, perhaps surprisingly, the rotors can be set up so as to 
make rotor-router walk on the d-regular tree recurrent. 

Proposition 4.2.12. Let T he the infinite d-regular tree. If all rotors initially point in 
direction d — 1, then every chip in an infinite succession of chips started at the origin 
eventually returns to the origin. 



Proof. By Lemma 4.2.11 for each n, the n-th chip sent to each principal branch Y returns 



to the origin before hitting height n + 1 of T. □ 

We continue our exploration of recurrence and transience on the infinite ternary 
tree T, but now we allow for arbitrary rotor configurations. We again focus on a single 
branch Y . For a given rotor configuration we define the escape sequence for the first n chips 
to be the binary word a = ai . . . a„, where for each j, 

0, if the j*^ chip returns to the origin; 

1, if the j*^ chip escapes to infinity. 

As noted above, a chip cannot stay within a finite height indefinitely without returning to 
the origin, so a is well-defined. 



We define a map 'ip on an escape sequence a = ai . . . a„. First we rewrite a as the 
concatenation of subwords bi - ■ - bm where each bj G {0, 10, 110}. Since at least one of any 
three consecutive chips entering Y is routed directly upwards by the root rotor, at most 
two of any three consecutive letters in an escape sequence o can be ones. Therefore, any 
escape sequence can be factored in this way up to the possible concatenation of an extra 0. 
Now we define ip{a) = (c, d) by 



{cj,dj) = 



(0,0), 


if 6, 


= 




(1,1), 


if&, 


= 110 




(0,1), 


iib, 


= 10 and #{i < j\bi 


= 10} is odd 


.(1,0), 


ifb, 


= 10 and #{i < j\bi 


= 10} is even 



Now we define the map cp on a pair of binary words c and d each of length m by 



122 



(f){c, d) = bi - ■ ■ bm, where 



= < 10, 







110, 




(0,0) 



(1,0) or (0, 1) 



(1,1) 



Note that (f>is a left inverse of i/j, i.e. (j) o ij;(a) = a, up to possible concatenation of an extra 



Lemma 4.2.13. Let Y be a principal branch of the infinite ternary tree. Fix a rotor 
configuration on Y with the root rotor pointing up. Let c and d be the escape sequences for 
the configurations on the left and right sub-branches of Y , respectively. Then (p{c, d) is the 
escape sequence for the full branch Y . 

Proof. We claim that each word bj is the escape sequence for the full rotation of the root 
rotor. Suppose first that {cj,dj) = (0, 0) so that bj = 0. In this case, since each sub-branch 
will return the next chip it sees back to r, the next chip that enters the full tree will be 
sent up to the origin, and the root rotor will once again be pointing up to the origin. If 
{cj,dj) = (1,0), the next chip entering the full tree will be routed to the left sub-branch 
where it escapes to infinity. The following chip will be routed to the right sub-branch and 
then back up to the origin, and the root rotor will once again be pointing up. In this case 
we have bj = 10. If {cj,dj) = (0, 1), the next chip entering the full tree will be routed to the 
left sub- branch, back up to r, and then to the right sub-branch where it escapes to infinity. 
The following chip will be routed directly up to the origin leaving the root rotor pointed 
up once again. Again, in this case bj = 10. Finally, if {cj,dj) = (1, 1), the next two chips 
entering the full branch will escape to infinity, one through each sub-branch. The following 
chip will be routed directly up to the origin, once again leaving the root rotor pointing up. 
In this case we have bj = 110. □ 

We observe that if the root rotor is not pointing up and c and d are the escape 
sequences of the left and right sub-branches respectively, then a = ^(c, d) is not the escape 
sequence of the full branch. In order to calculate the escape sequence of the full branch 
we define extended escape sequences on the sub- branches, c' and d' . Suppose the root rotor 
initially points to the left sub-branch. Then c' = Oc and d' = d. Suppose the root rotor 
initially points to the right sub-branch. Then d = Oc and d' = Od. Now a = 0(c', d') is the 
escape sequence of the full branch. 



0. 
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We now introduce the condition that is central to characterizing which words can 
be escape sequences: 

any subword of length 2^ — 1 contains at most 2^~^ ones (Pjt) 

We next show that the map V' preserves this requirement. 

Lemma 4.2.14. Let a he a binary word satisfying (P^) and let ip{a) = (c, d) as defined 
above. Then c and d each satisfy (Pfe_i). 

Proof. Let c/ be a subword of c of length 2^~^ — 1 and let d' be the corresponding subword 
of d. Let a' = (t){c', d'), which is a subword of aO. The formula for <p guarantees that a' has 
one zero for each letter of c' , so a' has 2*^-^-1 zeros. Since the last letter of a' is zero, and 
a satisfies (Pfe), it follows that a' has at most 2^~^ ones (else after truncating the final zero, 
the suffix of a' of length 2^ — 1 has at most 2^^^ — 2 zeros, hence at least 2^^^ + 1 ones). 

Let m be the number of ones in c'. Since the instances of (0, 1) and (1,0) alternate 
in the formula for ip^a) = (c, d), it follows that d' must have at least m — 1 ones. Since 
the number of ones in c' and d' combined equals the number of ones in a', we obtain 
2m — 1 < 2^~^, hence m < 2^~'^. The same argument with the roles of c and d reversed 
shows that d has at most 2^^"^ ones. □ 

Lemma 4.2.15. Let a = ai . . .an be a binary word of length n. Then a is an escape 
sequence for some rotor configuration on the infinite branch Y if and only if a satisfies (Pjt) 
for all k. 

Proof. Suppose a is an escape sequence. We prove that a satisfies (P^) for each k by 
induction on k. That a satisfies (Pi) is trivial. Now suppose that every escape sequence 
satisfies (Pfe-i) and let c and d be the extended escape sequences of the left and right sub- 
branches respectively. As we saw above, a = (j){c, d) up to the possible concatenation of an 
extra zero. Let a' be a subword of a of length 2^ — 1, and let il){a') = (c',d'). Then there 
are words c" and d" each of which is a subword of c or d, and which are equal to d and d', 
respectively, except possibly in the first letter; moreover the first letters satisfy c'-^ < c'{ and 
4 < d'l 

By the formula for ip, the number of ones in a' is the sum of the number of ones 
in d and d' . If c' has length at most 2^"^ — 1, then since c and d satisfy (Pjt_i), each of d 
and d' has at most 2^'"^ ones, and therefore a' has at most 2^~^ ones. On the other hand. 
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if d has length at least 2^~^, then the number of zeros in a! is at least 2^^"^ — 1. Thus a' 
has at most 2^~^ ones, so a satisfies (.Pfc). 

The proof of the converse is by induction on n. For n = 1 the statement is trivial. 
Suppose that every binary word of length n — 1 satisfying (P^) for each k is an escape 



sequence. Then by Lemma 4.2.14 ip^a) = (c, d) gives a pair of binary words each satisfying 
(i-fc) for all k. If c and d have length n — 1 or less, then they are escape sequences by 
induction, hence a is an escape sequence by Lemma 4.2. 13[ If c and d are of length n then 
the definition of ^p implies that aj = for all j, in which case a is an escape sequence by 



Proposition |4. 2. 12[ □ 

We can now establish our main result characterizing all possible escape sequences 
on the infinite ternary tree. 

Theorem 4.2.16. Let a = ai . . . a„ be a binary word. For j G {1,2,3} write a^^^ = 
ajttj+sajj^Q .... Then a is an escape sequence for some rotor configuration on the infinite 
ternary tree T if and only if each a^^^ satisfies (P^) for all k. 

Proof. Let Y'^^\ Y^'^\ and y^^^ be the three principal branches of T assigned so that the 
rotor at the origin initially points to Y'^^\ Then a is the escape sequence for T if and 
only if a^^'^ = ajaj^^aj^Q ... is the escape sequence for y(j). The result now follows from 
Lemma 14.2.151 □ 
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Chapter 5 



Conjectures and Open Problems 



5.1 Rotor-Router Aggregation 

A number of intriguing questions remain unanswered about rotor-router aggre- 



gation in Z . While we have shown in Theorem 1.2.1 that the asymptotic shape of the 



rotor-router model with a single point source of particles is a ball, the near perfect cir- 



cularity found in Figure 1.3 remains a mystery. In particular, we do not know whether 



an analogue of Theorem [1.2. 2 holds for the rotor-router model, with constant error in the 



radius as the number of particles grows. 



Equally mysterious are the patterns in the rotor directions evident in Figure 1.3 
The rotor directions can be viewed as values of the odometer function mod 2d, but our 
control of the odometer is not fine enough to provide useful information about the patterns. 
If the rescaled occupied region y^vr /n A„ is viewed as a subset of the complex plane, it 



appears that the monochromatic regions visible in Figure 1.3 in which all rotors point in 
the same direction, occur near points of the form (1-1-22;)""'^/^, where z = a+bi is a Gaussian 
integer (i.e. a,b E Z). We do not even have a heuristic explanation for this phenomenon. 



Figure 5.1 shows the image of ^1,000,000 under the map z l/-z^; the monochromatic 
patches in the transformed region occur at lattice points. 

Laszlo Lovasz (personal communication) has asked whether the occupied region 



An is simply connected, i.e. whether its complement is connected. While Theorem 1.2.1 
shows that An cannot have any holes far from the boundary, we cannot answer his question 
at present. 

Another question is whether our methods could be adapted to internal DLA to 
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Figure 5.1: Image of the rotor-router aggregate of one million particles under the map 
z 1-^ 1/z^. The colors represent the rotor directions. The white disc in the center is the 
image of the complement of the occupied region. 



show that if n = Wrfr*^, then with high probability Br--clogr C In, where In is the internal 
DLA cluster of n particles. The current best bound is due to Lawler |30J, who proves that 
with high probabiUty -Br_ri/3(iogr)2 C In- 

Finally, let us mention two simple variants of aggregation models which appear 
to have interesting limit shapes, but about which we cannot prove anything at present. 
The following conjecture addresses a variant of the rotor-router model. Analogous vari- 
ants can also be defined for internal DLA and the divisible sandpile, and by analogy with 



Theorem 1.1.2 we expect all three models to yield the same limiting shape. 



Conjecture. Let An be the region formed from rotor-router aggregation starting with 
n particles at the origin in Z^, if rotors on the positive x-axis reflect particles back in the 
direction they came from; all other rotors operate normally. Then as n — > oo, the asymptotic 
shape of the rescaled region n"^!"^ An is the cardioid in whose boundary is given by the 
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Figure 5.2: Rotor-router aggregation started from a point source at the origin in Z^, modi- 
fied so that sites on the positive j;-axis (a) reflect particles back in the direction they came 
from; or, (b) send all particles in the downward direction. 
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This apparent cardioid is pictured in Figure 5.2 on the left. If instead we alter 
the rotors on the positive x-axis so that they always point downward, we obtain the spiral- 
shaped region on the right. Our simulations indicate that a limiting shape exists for this 
variant as well, but we do not have a conjectured formula for the boundary curve. 



5.2 Sandpile Aggregation 



While Theorem 1.2.3 improves on the best-known bounds of [ST and [18] for the 
shape of sandpile aggregation in Z'^, it does not settle the question of existence of a limiting 
shape. There is also an intriguing pattern in the shapes of the sandpile aggregates Sn,H of 
n particles in which every site starts with a hole of depth H. These shapes are pictured in 



Figure 5.3 for n = 250,000 and three different values of H. 



Question. In dimension two, is the limiting shape n ^^'^Sn,H a regular polygon with 4H+12 
sides? 
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Figure 5.3: Abelian sandpile shape Sn,H started from a point source of n = 250, 000 particles 
in 7?, for three different values of the hole depth H. 

Simulations indicate a regular polygon with some rounding at the corners; it re- 
mains unclear if the rounded portions of the boundary become negligible in the limit. If 
the limiting shape is not a polygon, it would still be very interesting to establish the weaker 
statement that it has the dihedral symmetry -D4H+12- 

The only case that is even partly solved is H = —2: Fey and Redig [18] prove that 
Sn,-2 is a cube in Z'^, but the growth rate of the cube is still not known. It would be of 
interest to show that the volume of the cube Sn-2 is of order n. 
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